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Introduction 



Quantum Chromodynamics (QCD) is a quantum field theory, based on an SU{Nc) non 
abelian gauge group, born in order to describe strong interactions. Presently it is a well 
defined theory as well as the best candidate nowadays available. 

The success of QCD originated from the fact that it provided precisely, for Nc — the 
observed symmetries of strong interaction (such as the statistic of the baryons) and no 
more. Prom a dynamical point of view, its outstanding property of asymptotic freedom 
— due to the non abelian nature of the gauge group — was able to account for the scaling 
properties of cross sections experimentally observed. Purthermore, even if not rigorously 
proven, QCD gives strong indications of color confinement, e.g. in lattice simulations. 

Asymptotic freedom means that the effective coupling, as defined by the renormaliza- 
tion group, becomes vanishingly small when large space-like momenta (with respect to 
the QCD scale Aqcd — 200 MeV) are transferred. This kind of processes, called hard 
processes, can be investigated by means of perturbative methods. 

On the opposite side, soft scattering, hadronization and all long distance effects, un- 
avoidable in any strong reaction, involve strong coupling features most of which are far 
from present computational possibility. 

At the basis of several applications of perturbative QCD is the factorization. Por 
some measurable quantities, factorization theorems exist which allow the separation of 
short distance (perturbative) physics from long distance (non perturbative) physics of 
observable hadrons. It should be noted that what can be actually evaluated are not 
absolute values of observablcs corresponding to a given choice of variables, but rather the 
evolution of observables in the variables space. 

A crucial role in the factorization methods is played by the number and the rela- 
tive values of the hard scales involved in the process. In the traditional deeply inelastic 
scattering (DIS) and in the old accelerators physics, the virtuality of the transferred 
momentum is the only hard scale (the center of mass energy being of the same order 
of magnitude). Por this single-scale processes, perturbative QCD predicts the evolution 
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in of the relevant quantities as a power series in aJ^Q^^. The natural framework for 
studying this class of phenomena is the coUinear factorization in which only the longitu- 
dinal (with respect to the incoming hadron) degrees of freedom of the on-shell partons 
are present. The transverse degrees of freedom, peculiar of the interacting theory, give 
rise to logarithms of Q^. These large logarithms, which need to be taken into account 
to all orders in perturbation theory, are resummed by the DGLAP equations and are 
responsible for the scaling violations, i.e., for the deviation from the pure scaling behavior 
one would obtain by neglecting the parton-parton interaction. 

The coming of high-energy (for the time being) colliders has entered a new era in which 
the energy is a scale much harder than the transferred momentum. The electron- 
proton collider HERA at DESY is of particular importance, since high energy DIS has 
opened the road to new and interesting physics. The large available energy in the center 
of mass of the coUiding particles s ~ (300 GeV)^ allows to investigate a wide kinematic 
region. It is possible to reach values of larger than 10^ GeV^ and very small values 
of the Bjorken variable x ~ 10~^. The structure functions have shown to undergo large 
scaling violations towards high Q^, especially at low values of x. Besides, a steep rise of 
the structure functions has been observed stimulating the interest of a considerable part 
of both the experimental and the theoretical community. 

In this context, the perturbative QCD description by means of the DGLAP approach 
reveals itself very successful, even beyond the expectations, in the sense that starting with 
reasonable parametrizations of the parton distribution functions at rather low values of 
Q\ ~ 0.35 GeV^ — definitely outside the perturbative domain — the evolution of the 
structure functions is very well described by a next-to-leading order DGLAP fit. 

The basic issue remains to justify or motivate the particular shape of the input parton 
densities entering the DGLAP evolution equation. The most pessimistic approach is that 
the problem is a non perturbative matter, since the perturbative evolution works even 
with initial conditions at very low Q\. And actually it could be so. 

However one can also argue that we are not allowed to use DGLAP equations outside 
the perturbative domain and that it would be preferable to start the evolution at some 
higher point of Q\ of the order of some GeV^, where the perturbative theory is presumably 
trustworthy. In this case the right input functions which are needed to describe the data 
present a (small) power-like rise in \jx. For instance, the gluon distribution — playing 
a leading role in high energy processes — has the form f^^\x,Q^) oc a;^'^^*^^-' where 
A(gg ~4 GeV2) -0.17. 

The question then arise: can one justify such a power-like shape for the partonic 
densities? The question reminds us of Regge theory, where it is expected a power-like 
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growth of the cross sections with s (the latter being proportional to l/x). Even if Regge 
theory is not based on perturbative physics, nevertheless they should be someway related. 
This hypothesis is confirmed by the fact that, at very low values of Q^, the x-growth 
exponent of the structure functions reaches the value A ~ 0.1, and this suggests a smooth 
junction with the soft pomeron value cjp ~ 0.08, i.e., the universal exponent governing 
the s-growth of total cross sections. 

It remains to see whether perturbative QCD can explain a small- a; rise of the structure 
functions with a power-like law and with a correct exponent in the intermediate ~ 
1 10 GeV^ regime. 

In connection with that point there are the so-called two-scale processes, such as 7*7* 
scattering and forward jets, where at the endpoint of the QCD evolution two hard scales 
are present. Here the use of perturbative theory should be more suitable in order to 
describe the and energy behaviour. Prehminary results seem to indicate a growth in 
energy compatible with s^'^'^^\ A(Q^) ~ 0.3 for of order of 3 ^ 40 GeV^. 

Those high cnergy-not very large regimes arc referred to as semi-hardiegiTaes. Here 
the perturbative series can be slowly converging or even unstable, because the higher order 
terms are accompanied by large Ins and may be as important as the first ones. In this 
case, finite order calculations are no longer reliable and rcsummation techniques must be 
devised in order to take into account all the important terms. 

Different approaches have been adopted in order to study gauge theories at asymp- 
totic energies, e.g., Regge theory. This "old" theoretical problem was for the first time 
investigated at a perturbative level in the 70's by the russian school where the large loga- 
rithms of the energy are resummed by means of the BFKL equation, predicting at low x a 
power-like growth of the structure functions, but with a too large exponent X{as) — 0.55 
for values of ctg ^ 0.2 as in the HERA range. 

In order to obtain a quantitative agreement with the experimental data, a huge theo- 
retical effort ('89 -'98) has been devoted to analyse high-energy QCD beyond the leading- 
logarithmic approximation. Last year the next-to-leading logarithmic (NLx) BFKL ker- 
nel was found. It immediately appeared that the NLx corrections to the kernel are quite 
large and of particular relevance for their negative sign. In particular, the "pomeron 
singularity" , which should provide the s-growth exponent for high energy cross sections, 
undergoes a drastic reduction reaching its maximum value (~ 0.11) for quite a small value 
of as — 0.08 and in the HERA range it becomes even negative. 

A pathological consequence of the large size of the NLx corrections is that, for larger 
values of the coupling, they may even provide negative cross sections in the very large 
s- limit! The fact is that, for values of ctg > 0.1, such corrections are so large that they 



4 



cannot be taken literally, suggesting an intrinsic instability of the perturbative series in 
the effective parameter aglnl/x. 

A direct analysis shows that responsibles for this instability are the large coUinear 
contributions to the kernel. The latter are single or double logarithms of the transverse 
scales and determining the process. Single logarithms (as In fe^/feg) are strictly 
related to the well known logarithmic corrections to scaling in the naive parton model, 
providing the scaling violation predicted by the renormalization group. Double logarithms 

In^fcVfc^) are due to the mismatch between the factorization scale Sq = kk^ entering 
the high energy factorization formula and the Bjorken scale sb = Max{k'^,kl) which is 
the relevant scale in the coUinear limits k^ ^ k^ and fc^ <^ k^. 

Since coUinear contributions are determined by the renormalization group, it is manda- 
tory to develop a unified picture where both renormalization group constraints and small-x 
features are taken into account. Single and double Infe^/fep are known to all orders in 
In Therefore, the first step of the improved small-x formulation consists in the resum- 
mation of the coUinear contributions to the BFKL kernel, consistently with the full Lx 
and NLx expressions. The second step concerns the new method for solving the improved 
small- a; equation. Starting from the factorization property of the solution in a pertur- 
bative times a non perturbative factor (up to higher twist corrections) the perturbative 
part is determined by an expansion with respect to a new parameter u — the moment 
index with respect to x. The use of u instead of ctg turns out to be more convenient in 
the small-x region where the whole Q^-dependence of as(Q^) is important and has to be 
taken into account. 

The outline of the thesis is the following: the first chapter introduces the basic objects 
of our study and recalls the coUinear properties of QCD in the special context of DIS. 
Scaling violations are explained both with the operator product expansion analysis and 
in terms of the renormalization group improved parton model, in order to identify the 
results of the latter with the field-theoretical quantities (e.g. the anomalous dimension) 
of the former. 

In Chap. we concentrate on high energy physics by starting with an overview of 
Regge theory in connection with the phenomenological analysis of processes at asymptotic 
energies. As far as the perturbative treatment of small-x processes is concerned, we recall 
the framework of high energy factorization, its connection to the coUinear one and we 
show how the resummation of the leading \nl/x leads to the BFKL equation. 

The subsequent chapters contain the author's original contributions, obtained in col- 
laboration with M. Ciafaloni and G.P. Salam and published in Refs. P7|, ^|5^ ,^ . 



In Chap. 1^ we present the generalization of the high energy factorization formula in 
NLx approximation. In this context we perform the NLx analysis for both the impact 
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factors 1^ and the BFKL kernel. We present the main results as well as the problems 
concerned with the large NLx corrections. 

With Chap. ^ we get to the heart of the improving procedure of the NLx approxima- 
tion. The starting point is the identification of the collinearly-enhanced physical contri- 
butions as the main responsibles of the instability of the BFKL hierarchy. The collinear 
singularity resummation procedure devised in Refs. [^,^ leads to the improved small-x 
equation which incorporates exact leading and NLx BFKL kernels on one hand and renor- 
malization group constraints in the relevant collinear limits on the other. The basic idea 
consists in using a new expansion parameter such that, in the corresponding series, the 
terms of order higher than the second are indeed small corrections [Q. 

Besides exhibiting the good qualities of both DGLAP and BFKL formulations, this 
improved small-x equation shows several tricky aspects mainly due to the running of the 
coupling, the need of regularizing the Landau pole, the diffusion in the infrared region 
and the existence of two perturbative exponents for the high-energy behaviour. For in- 
vestigating such questions, a simple model [|^, built up with the main ingredients of 
the RG-improved formulation, has been introduced and studied in Chap. H. Within this 
model, many quantities (among which the full solution of the ensuing equation) can be 
computed exactly — at least numerically — and the strong coupling physics can be taken 
into account so as to see its reflection on the relevant outputs and a clear interpretation of 
the small-x growth exponents is at hand. In two-scale processes, this method can find a 
direct application, at least in an intermediate small-x, large-Q^ regime where perturbative 
physics appears to dominate. For increasing values of the energy, there are signals of a 
new transition mechanism in consequence of which the perturbative features are lost and 
"the pomeron", i.e. the energy growth exponent, belongs to the realm of non perturbative 
physics. 

The main outcome of the work presented in this thesis is the development of a novel 
small-x expansion which turns out to be much more stable than the BFKL method and 
which is affected by small next-to-next-to-leading corrections. In particular, the improved 
method provides a "hard pomeron" index ujs{as) — 0.27 ^ 0.32 for ctg — 0.2 ^ 0.3 which 
should be observable in two-scale hard processes. 

We are confident that, in two-scale processes, this method can serve to make reliable 
predictions, at least in the perturbative regime previously mentioned. Whether or not it 
applies to single-scale processes like DIS, remains a question to be investigated in more 
detail. 



Chapter 1 



Perturbative Treatment Of Hard 
Processes 

1.1 What's asymptotic freedom? 

The systematic study of strong interactions and the successful description of deeply inelas- 
tic phenomena in terms of free partons inside the hadrons demanded an asymptotically 
free theory, i.e., a theory with vanishingly small coupling at short distances, to be the 
right and fundamental one. 

Non abelian gauge theories are unique among renormalizable theories in 4 space-time 
dimensions in having this peculiar feature. The quantum charge giving rise to strong 
interaction, called color, is carried by both the quarks, the fermionic spin 1/2 matter 
fields, and the gluons, the bosonic spin 1 gauge fields, mediators of the color interaction. 

Hadronic interactions are effectively described by quantum chromodynamics (QCD) 
which is a gauge theory based on a SU{3) group (for a review, see e.g. ^ and also 
0). QCD is a conformally invariant theory. Nevertheless, the renormalization procedure, 
required in order to remove the divergences of the perturbative expansion, breaks the 
conformal invariance and causes the dimensionless coupling g to depend on dimensional 
kinematic variables. However, contrary to abelian theories, the self interaction of gluons 
"spread out" the color charge and an anti-shielding effect is produced. This is responsible 
of weakening the color charge at small distances and presumably to provide the strong 
coupling regime at large distances leading to confinement. 

In practice, for a process characterized by the scale Q^, the color strength := g'^/4:TT 
among partons must depend on in a way constrained by the RG equations, according 
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to which 

PMQ')) ■■= = -boaliQ') - haim + ■ ■ ■ , (1.1) 

where the 6j's are pure numbers completely determined by the structure of the gauge 
group (and, for i ^ 2, from the renormalization scheme). 

If 6o > asymptotic freedom holds. In QCD 

bo = ^ , 1.2 

" 127r ' ^ ' 

where A^^c = 3 is the number of colors and = 3 -i- 6 is the number of flavors with masses 
less than Q, and hence Bq ~ 0.72 ^ 0.88. 

The solution of Eq. ( |1 . 1| ) is known and at lowest order reads 

«s(Q') = ^ /q. , (1.3) 
bo In -ri — 

where Aqcd (A from now on) is a new mass parameter — not present in the original 
Lagrangian — introduced by the renormalization procedure. It is apparent that, for 
^ A^, as <^ 1 so perturbation theory and Eq. ( |1.3| ) are reliable, while for < A^ we 
are in strong coupling regime. The QCD scale A has to be determined experimentally. 
Its value of about 200 MeV sets the lower bound of the perturbative domain. 



1.2 Hard and semi- hard processes 

According to the above considerations, we define hard processes those hadronic pro- 
cesses involving a large momentum transfer ^ A^, i.e., greater than or equal to 
some GeV^. In this case, perturbative QCD provides quantitative predictions. 

Among hard processes we must distinguish semi-hard ones where various (hard) 
momentum scales of different orders are present, because the coefficient of the perturbative 
series may become large and the latter has to be improved with resummation techniques 
of classes of Feynman diagrams. 

In this chapter we will consider only hard processes. Semi-hard processes will be 
discussed from Chap. 0. 
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1.3 Inclusive means perturbative 



As soon as we start a perturbative calculation, we are confronted with divergent integrals. 
Besides the ultra-violet (UV) divergences, causing the running of the coupling constant 
and, as we shall see, the scaling violations of structure functions (SF), the perturbative 
theory suffers also infra-red (IR) divergences, due to the presence of massless particlesQ. 

IR singularities appear in the calculation of diagrams representing 



• emission of gluons with vanishing 4-momentum {soft singularities;) 

• emission of massless partons coUinearly to the incoming one (co//mear singularities). 



How to recover finite estimates for cross sections? We have to realize that, from a physical 
point of view, it is not possible to distinguish a single quark state from one in which the 
quark is accompanied by an arbitrary number of softQ or collinear gluons. These degen- 
erate states belong to the same physical state. When evaluating transition rates, both 
the initial and final states should be physical ones. Then, the Kinoshita-Lee-Nauenberg 
theorem p| — stating that in a theory with massless fields, transition rates are free of the 
IR (soft and collinear) divergence if the summation over the initial and final degenerate 
states is carried out — ensures to obtain finite results in completely inclusive processes, 
i.e., processes in which nor quarks neither gluons are registered in the initial or final states. 
We can quote, e.g., the total e~ annihilation cross section, the jet cross section at fixed 
angular and energy resolution, the energy flow, etc.. These quantities are calculable in 
perturbative QCD by assuming that the sum over parton states equals the corresponding 
hadronic sum, in the same spirit of the parton model. In short, hadronization and all long 
distance effects are unimportant when neglecting the structure of the hadronic initial and 
final states. 

When one or several partons are identified by measuring their momentum, the hy- 
potheses of the KLN theorem are no longer fulfilled and the perturbative calculation 
exhibits large logarithms due to collinear singularities. However, thanks to the factor- 



ization theorem of collinear singularities (cfr. Sec. |1.7|) , this logarithms can be resummed 
{improved perturbation theory) and they give rise to non trivial anomalous dimensions 
determining scaling violations. We refer to these processes as inclusive processes with 
registered partons. We note that the soft singularities of virtual corrections and real soft 



^In hard processes, light quarks can be considered massless in comparison with the much larger hard 
scale Q^. 

^I.e., in the limit of vanishing energy. 
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emissions cancel each other because of the inclusive nature of the process. A few ex- 
amples: DIS (the registered quark in the initial state being the active parton), inclusive 
single particle distribution in e"*" e~ annihilation (the registered quark being in the final 
state), jet cross section at fixed invariant mass of observed particles, etc.. 

Finally, when one or more hadrons are identified, strong coupling effects become im- 
portant. These exclusive or semi-inclusive processes are therefore sensitive to IR physics 
and require a particular perturbative treatment (see, e.g., 0) in addition to hadronization 
models. 

The aim of this thesis is to test and to improve high-energy perturbative QCD. We 
will be concerned mainly with high-energy DIS, double DIS, forward jets and in general 
with the class of inclusive processes with registered parton described above. 



1.4 Deeply inelastic scattering 

In order to introduce one of the main physical processes studied in the field of strong 
interactions as well as to present some basic mathematical apparatus relating physical 
observables in terms of field-theoretical quantities, we start by specifying DIS. 

An effective way to study the structure of the hadrons is to probe them by means of 
point-like particles, such as electrons or photons. In DIS of electrons by proton^, e.g. in 
the HERA configuration, a beam of electrons of momentum pi collides on a target beam 
of protons of momentum p2- By recording only the scattered electron p[ without care of 
the hadronic final state, one obtains a highly inclusive measure. The conventional Lorentz 



invariant kinematic variables are defined, in the notations of Fig. |1 . 1| , as follows: 



q : = 


Pi - p'l 


momentum of the virtual photon probing the proton 


s : = 


{Pl+P2)^ 


center of mass energy squared; 


: = 




virtuality of the photon; 


W': = 


{P2 + q? 


invariant mass of the hadronic shower X; 


u : = 


P2 ■ q/M 


electron transferred energy in the LAB frame; 


X : = 


QV^P2 ■ q 


Bjorken variable; 


y ■= 


P2 ■ q/pi ■P2 


transferred energy fraction in the LAB frame; 



^Hereafter we will consider a proton as the target, but any hadron can be considered as well. 
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X 



Figure 1.1: DIS: pi (p[) represents the incoming (scattered) electron, p2 the incoming proton, 
q the exchanged virtual photon and X the hadronic shower. 

If the target proton has no spin polarization, the differential cross section can be 



What can we say about the structure functions? Even if the final state is inclusive as 
far as the hadrons are concerned and we limit ourselves to sufficiently high transferred 
momenta Q^, long distance effects are important in this kind of reaction. This is due to 
the fact that, in the initial state, a particular hadron (the proton) is recorded by measuring 
its momentum and therefore the process is not completely inclusive — it belongs to the 
second-type ones described in the previous section. 

Suppose anyhow to undertake the perturbative calculation by assuming the proton 
composed of some quarks and gluons. At tree level, the calculation is lead back to the 
computation of the electron-quark Born cross section, but as soon as one consider 1-loop 
corrections, gluons set in and IR singularities (besides the UV ones) appear giving infinite 
cross sections. 

Let's stop for a while and try to understand what's happening by looking at scattering 
in the center of mass (CM) frame. At high energy ^/s ^ M, the fast motion causes the 
proton to be Lorentz contracted in the direction of the collision and its internal interaction 
are time-dilated. Thus, the lifetime of any virtual partonic state is lengthened and the 

^Here we are concerned only in the EM part of the neutral current, neglecting exchanges. 
^We have neglected a term of relative order NP / s in the coefficient of i^2- 
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time it takes the electron to traverse the proton is shortened. When the latter is much 
shorter than the former, the proton will be in a single virtual state characterized by a 
definite number of partons during the entire time the electron takes to cross it. Since 
the partons in practice do not interact during this time, each one may be thought of as 
carrying a definite fraction z of the proton's momentum^. 

Furthermore, if we consider large momentum transfers 3> M^, then the resolving 
power of the virtual photon is much smaller than the transverse size of the proton and 
the electron will be able to interact with only a single parton (provided the density of 
partons is not too high). 

Under these conditions, the electron-parton amplitudes contribute incoherently, hence 
the scattering probability on the proton is given by the sum of the probabilities of scat- 
tering on the partons. The latter, in turn, can be thought to be inside the proton with 
certain probability distributions. 

We have finally the following simple probabilistic interpretation: the electron-proton 
cross section is given by the sum over partons of the probabilities for electron-parton 
scattering, each of them weighted with the density distribution of that parton to be 
found inside the proton. This is the essence of the factorization theorem of collinear 
singularities which, for the structure function, states 

Ff\xM') = Y^\\zf^l^\zM')FfC-.alQ^)) (z = L, 2) . (1.5) 

Here f^^^^\z, Q'^)dz is interpreted as the probabilityQto find a parton of type a (= gluon, 
u, u,d,d,s, . . .) in a hadron of type H carrying a fraction z to z + dz of the hadron's 
momentum {p2 = p in the following), while Fj^^\x , aslQ"^)) is the partonic structure 
function (for the process ea eX) with respect to the rescaled Bjorken variable^ 

2{zp) ■ q z 

Note that the partonic distribution functions (PDF) have an explicit Q^-dependence. In 
our intuitive picture this is due to the fact that additional partons — generated by virtual 
processes — can be resolved by increasing the value of Q^. 

^The CM frame is, in this case, a very good approximation of the so called "infinite momentum 
frame" in which the proton is highly boosted so to have only a (very large) light-cone component of the 
momentum. 

^In the literature an alternative notation is also used, differing for a z factor from ours. For instance, 
the gluon density reads f^^\z, Q^) — zg{z, Q^) etc.. 

®The hat* indicates partonic quantities, e.g., p — zp is the struck parton momentum. 
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1.5 The parton model 



If we completely neglect interactions among partons, the parton content of the proton is 
fixed, i.e., independent of Q^, since virtual processes are absent. Furthermore, the virtual 
photon can be absorbed only elastically by quarks, as depicted in Fig. |1.2|. Since partons 




Figure 1.2: Elastic scattering of a quark by absorption of a virtual photon. 

are considered massless objects in DIS, it holds 

= {zp + qY = {zpY + + 2zpq = z 2pq — 

or, equivalently 

, = 9L = j: =^ x = - = l 
2pq z 

This means that the Bjorken variable x coincides with the momentum fraction of the active 

quark and that the partonic SF are proportional to 6{1 — x). The explicit calculation of 

the partonic SF yields (omitting the hadron suffix H hereafter) 

Fi^\x,a, = 0) = elS{l- x) , (1.6a) 
Fi'\x,a, = 0) = 0, (1.6b) 

Ca being the EM charge (in unit of the electron charge) of type a partons. 

The last equation (Callan-Gross relation) follows from the fact that a longitudinally 
polarized vector boson cannot be absorbed by a spin-1/2 quark (cfr. App. However 
the QCD corrections provide a non-zero coupling between longitudinal virtual photons 
and partons, so that Fl = 0{as). The SF measurements show that Fl ^ F2 (see Fig.p,.3[), 
confirming the spin 1/2 property of quarks. 

Because of the Q^-independence of both the PDF and the partonic SF, the SF depend 
only on the dimensionless variable x, as one would expect from a simple dimensional 
analysis. This is the so called scaling behavior ||^ . The Q"^ dependence of F2 for various 
values of x is reported in Fig. |1.4| . Approximate scaling is visible, especially around 
X ~ 0.2. Anyway, scaling violations are also evident, and their origin is the subject of the 
next sections. 
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Figure 1.3: The longitudinal proton structure function Fl data compared with the F2 structure 
function fit (the black line). 
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Figure 1.4: The proton structure function F2 versus for various values of x. 
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1.6 Scaling Violations 

Historically, the first approach to next-to-leading order^ (NLO) corrections for DIS made 
use of the OPE. Even if a clear interpretation of the underlying physics will be available 
in the parton model picture, this formal approach will be useful in order (i) to classify the 
main contributions to scaling violations (by means of the so called twist index); (ii) to show 
explicitly the connection between anomalous dimensions of certain composite operators 
and evolution kernels for partonic distributions; (iii) to provide a first explanation of 
factorization of perturbative and non perturbative effects in strong interactions. 

1.6.1 Operator product expansion (OPE) 

Starting from the expression ( [A.4b| ) for the hadronic tensor 



W,.{p,q) = d^^e'^^i J2 (P,A|[J^(^), J.(0)]b,A) , (1.7) 

A=spin 

we note that the commutator of the EM currents vanishes for z'^ < 0. Furthermore, in 
the Bjorken limit oo, x fixed, the dominant domain of integration turns out to be 

^ < const/Q^, i.e., a thin sheet upon the forward light cone. 

Notice that the product of two local operators at light-like distances is singular. Never- 
theless, these singularities can be embodied in c-number functions by mean of the OPE ^ 
that near the light-cone takes the form|^ 

oo 

J{z)J{0) = J2Y. 0%..^^ , (1.8) 

0^=0 j 

where the c-number coefficient functions C^\z'^) are eventually singular at = and 
the composite operators 0^\^...^„ are well-defined objects, apart from the renormalization 
infinities. Anyhow, the singular behaviour in coordinate space is entirely carried by the 
coefficient functions. 

In the Bjorken limit, all the masses can be neglected and the only dimensional variables 
in configuration space are the coordinates. A naive dimensional analysis of Eq. (|1.8|) shows 
that, near the light-cone, the behaviour of the coefficient functions should be 

C'^\z^) ~ (^s^-s+^CrfL--) ^ ^2^0, (1.9) 



^Next-to- leading in as. 

^''We omit for simplicity the Lorentz and flavour indices. 
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where 3 = [J] and c?!^ = are the natural dimensions of the EM current and of 

the composite operators respectively. 

In the non interacting theory the conformal invariance would require such a simple 
power-like behaviour. However in the interacting theory the above simple argument does 
not hold and in general logarithmic corrections to Eq. (|1.9|) develop after the renormal- 
ization procedure. 



We see from Eq. (L9) that, the smaller the twist value d'^ — uj is, the more singular 
c!^^ will be, and then the dominant terms as (Q^ —>■ oo) are controlled by the 

lowest twist operators, higher twists being suppressed by powers of Q^. 

There are two kinds of leading twist operators (twist = 2) occurring in the OPE ( |1.8| ) 
for DIS: 



O^h^ = F,,xD,, ■ ■ ■ 1),.-,F,^' , (1.10a) 
OS-M. =V^^V^^..---^^..V^^'^ , (f=l,...,iV^), (1.10b) 



Nf being the number of (massless) flavours and D ("D) the covariant derivative for gluons 
(quarks). The former are gluonic operators and the latter are quark operators. 

It is to be stressed that the composite operators of Eqs. ( |1.10|) appearing in the OPE 



do not depend on the hard probe (photon, weak boson, gluon, etc.) coupled to the 
hadronic tensor, and their hadronic matrix elements A^^ (cfr. Eq. ( |1.13|) ) describe only 
the partonic content of the initial particle H. In this sense, they are universal, i.e., 
the same for any hard reaction involving the same hadron H. However, the composite 
operators describe the long distance properties of the particles and therefore they cannot 
be calculable within perturbative QCD and their ignorance have to be compensated by 
some modellistic assumption or, better, by experimental measurements. 

The coupling to the probe is accounted for by the coefficient functions, which are 
process dependent. They contains the informations of the short distance interactions 
between the exchanged boson and the constituents of the hadron as well as of high transfer 
momentum processes among partons. The following analysis will show how the former 
part of these short distance processes corresponds to the partonic SF and the latter one 
can be included in the PDF which acquire a (logarithmic) dependence on the virtuality 

of the probe. 

What is needed to derive the Q^-dependence of the SF is 

• to relate the objects of the OPE to the SF; 

• to determine the Q^-dependence of the coefficient functions (the composite operators 
being independent of z and hence of Q^). 
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The relations between OPE and SF have been derived in |p and express the moments 
(i.e. the Melhn transforms) of the SF|^ 

dr 

FM)--= / —x-F{x,Q^) (1.11) 

in terms of the Fourier transforms^ of the coefficient functions 

j d^ze^'^'z'"' ■■■z''-C2\z^) =:q'''---q^'-c!^l^{Q^) (1.12) 
and of the matrix elements 

^ MO^il-.f^Jp, A) = A^LiP^i ■ ■ -Pf,^ + terms containing (1.13) 

A 

in the famous moment sum rules 

^u,(Q^) = E^^^C'«(Q2). (1.14) 



J 
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The determination of the coefficient function dependence on Q can be done by using 
the RG technique [10|. The renormalization procedure, necessary to handle the infinities 
stemming from the perturbative expansion, requires the introduction of a mass parameter 
/i (the renormalization scale) not present in the original lagrangian. As fi is an artificial 
parameter, the physical quantities — in our case the SF — do not depend on it. This con- 
straint determine the /^-dependence of the coefficient functions, of the composite operators 
and of the coupling as well. 

In addition, the dimensionless coefficient functions and the coupling can depend on 
the massive parameter and fi only through the ratio Q^/yU^- Therefore, their Q^- 
dependence is completely calculable within perturbative QCD and, by means of Eq. ( |1.14D , 
we obtain the Q'^ dependence of the SF. 

In particular, the RG equation for the coefficient function is 

(^^l. + 2(3ia^)^^5-^' + {jJ' CW(QV/^^a,(/.2))=0 Vj , (1.15) 

where the indices j and k label the different types of operators in Eqs. (|1.10| ) and is 
(the opposite of) their anomalous dimension matrix: 

7. := -ZZ' (-"1;^-) > Og:^-: = [Zj' o«.;— 1^^^' , (1.16) 

(1) 2 

Tc. =: "s Tc. + Tc. + ■ ■ ■ ■ (1-17) 



^^You can think F to be any of the SF F2 or F^. Only the coefhcient functions carry the "2" or "L" 
dependence. 

^^Up to a numerical factor 2™7r'^i'*. 
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The solution of Eq. ( p..l5| ) is found by introducing the running couphng function of 
Eq. (|1.1|) and can be written as a path-ordered exponential 



C(j)(QV/^',«s(/i')) = Texp 



jk 



da CW(l,as(g^)) 



:i.l8) 



showing that the coefficient functions depend on only through the running coupling 



1.6.2 The "running" parton distribution functions 



The formal results ( |1.14| ) and ( |1.18| ) can be given a simple and elegant interpretation in 
terms of PDF which depend upon Q^: we can write at LO in InQ^ 



Ifui/bo 



Cl^ihO) X [l + 0(l/lng^ 



jk 



:i.l9) 



First of all, switch off for a minute the strong interaction, so that QsiQ"^) = «s(At^) = 



and 



[1.20) 



If we take the Mellin transform of Eq. ( |1.5D by defining 

FS{c^m)--= ['-x-Ft\x,a.{Q')) 
Jo X 

fL'\Q')--= rdxx-/(^)(x,g2), 

Jo 

the moments of the SF turns out to be simply 



(1.21) 
(1.22) 

(1.23) 



In this situation it is natural to identify the short distance and long distance factors as 



:i.24) 
:i.25) 



Equivalently, one can obtain the parton model picture ( |1.5| ) by performing the inverse 
Mellin transform of Eq. ( 1.2CI|) by suitably identifying PDF and partonic SF. 
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We can go further by applying the same reasoning in the interacting case, in which we 
use again Eq. ( |1.24| ) — because at LO the partonic SF are Q^-independent — together 
with the new identification 



a) 



(1) 



'~/u:/bo 



^e\Q')- (1-26) 



ab 



in such a way that the Q^-dependence is incorporated in the PDF. 

At this point we differentiate the above relation with respect to In Q"^ and, taking into 
account Eqs. (|1.1|) and (|1.3|) , we end up with the evolution equation for the moments of 
the PDF 

d 



dlng^ 



fL'\Q') = J2<Q'^(^-^^'f^'^(Q'^ ' ^ = f'"f'g- (1-27) 

b 

This set of 2Nf + 1 coupled differential equations can be diagonalized by analysing the 
properties of the PDF under flavour symmetry — which is an exact symmetry if quark 
masses are neglected. 

Clearly, the gluon density f^^^ is invariant under flavour transformations, as well as 
the "modulus" of the quark vector 

Q') := J2 f^'\^^ Q") + f''\^^ Q") (1-28) 
f 

called quark singlet density. 

The non-singlet (NS) components of the quark densities, transforming according to 
the adjoint representation of SU{Nf), are constructed by subtracting from the various 
fermionic PDF the common singlet component: 

/(NS,0(^,g2):=/(f)(^,g2)__l_^(E)(^^g2) _ (^ 29) 

Since the NS quark densities carries different quantum numbers (flavours) they don't mix 
and renormalize independently. Furthermore, the flavour group commutes with the colour 
group, so all NS densities evolve with the same quark-to-quark anomalous dimension 7'''' 

^ fL''''HQ') = a.miiyL''''HQ')- (1-30) 



dlnQ2^" 

Only quark singlet and gluon operators mix in the renormalization and require a 2 x 2 
anomalous dimension matrix such that 
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Using Eq. (|1.22|) and the splitting functions implicitly defined by 



=: dz P^\z) , (1.32) 



27r 

we can invert Eqs. ( |1.30|) and (|1.31| ) in x space and reproduce the famous Dokshitzer 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations 

d 



dlnQ2 
d 

dlng2 



/(Ns.f)(^,g2^ = ^£M rd.P^^(^)/(^s-f)(.,Q^), (1.33a) 
27r ^z' 

f^^\x,Q') rd.p-(^)/(^)(.,g^) + 2iv,p^«(^)/fe)(.,g^), 

ZTT J ^ z z 



1.33b) 



f\z%\-) f''\z,Q') +P^^{-) f^\z,Q') . (1.33c) 



dlnQ^' ' ' 27r 

It is also possible to diagonalize the anomalous dimension matrix of the singlet sector by 
determining the eigenvalues 



7^ = 1 (tS^ + 7f ± V^(7S^ - 75^)^ + 8Nn'.V^ (1.34) 



and the eigenvectors 



e^ = fL'^+^§^e\Q'), (1-35) 



so that 

d 



dlng^ 



fL''\Q') = MQ'hifL''\Q')- (1-36) 



To summarize the results of the previous sections, we have shown that a certain 
number of observables of hard inclusive processes can be expressed as a convolution of 
process dependent, perturbative calculable, partonic cross sections and universal, non 
perturbative, partonic distribution functions inside hadrons whose dependence on the 
virtuality of the hard probe is however computable within the framework of pertur- 
bative QCD. This means that we are not in such a position to say what the value of 
F2{x = 0.1,Q'^ = 10 GeV^) is, but if we measure with enough accuracy the x-dependence 
of a sufficient number of structure functions at a certain value of ^ A^, we can predict 
their behaviour at higher values. Furthermore, the consequent knowledge of the PDF 
allows us to predict absolute values of cross sections for completely different hadronic pro- 
cesses. This is very important as a test of the theory as well as experimentally, in order to 
check normalization on the cross sections, to estimate the occurrence of yet unobserved 
phenomena, etc.. 
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Ta[z' + {l-zf] 


+ 3u; + 4 
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2Ca 


+n ^ + ^(1-2) +S{i-z) —^-—L 
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Ca = N, .,Cp = -j-^ ■ ^« = 2 ■ / ■= 1 [/(^) - /(l)l.'K^) - /(I) giz) dz 





Table 1.1: The 1-loop partonic splitting functions and anomalous dimensions. 



1.7 The improved part on model 

In the following, we shall show an alternative approach for dealing with scaling violations 
in DIS. This method presents several advantages with respect to the OPE shown before. 

On one hand, there arc very few reactions where one can justify the use of the OPE. It 
is therefore of great importance to rephrase the physics in the language of a very general 
QCD-improved parton model by which we can study a larger variety of phenomena. 

On the other hand the improved parton model gives the LO results in a much more 
intuitive and economical way and permits, with reasonable efforts, to perform NLO cal- 
culations (e.g., the NLO anomalous dimensions). 

In addition, — and this is the main reason of this lengthy discussion — the partonic 
framework clarifies the connection between general field-theoretical predictions, such as 
factorization and OPE, and the resummation of certain classes of Feynman diagrams. A 
very similar formalism will be very useful when dealing with high energy reactions, where 
Regge theory, BFKL dynamics, high energy factorization and RG constraints meet in 
describing semi-hard processes. Also in this case the partonic picture turns out to be of 
fundamental importance to extend the calculation to ncxt-to-lcading level and to take 
into account a whole series of sub-leading contributions to all orders, which is the central 
purpose of the present thesis. 
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1.7.1 Next-to-leading order parton model 



Having convinced ourselves of the factorization properties splitting strong coupling from 
perturbative physics, we are going to attack the problem from another point of view. 
We assume the hadron undergoing DIS to be composed of several partons a = f, f, g 
with distribution densities f^^\^), ^ being the longitudinal momentum fraction of the 
parton. We assume also that parton-photon and parton-parton interactions can be treated 
perturbatively. The validity conditions and consequences of these hypotheses are what 
we wish to present. 

Since among partons only quarks carry EM charge, the partonic interaction at Born 
level is represented uniquely by the diagram in Fig. |L^, where p = is the incoming 
quark momentum. The ensuing partonic tensor reads 



Wl^f\p,q) = el5il-x)T^,iP,q) 

r \ ^ ( , 1 pq \ ( ^ pq 

^Ap^ ^) = - 2 y^^'' '^)^PqY^~^'^^) " 

whose longitudinal and transverse projections give the partonic SF of Eqs. 



(1.37a) 
;i.37b) 



61). 





(a,0) 



Figure 1.5: The diagram representing the partonic tensor W 
incoming quark flavour an p its momentum. 



at Born level, a labels the 



At NLO, we should take into account all the insertions of a gluonic line in Fig. Eo 



such as Fig. Ola and also gluon-photon coupling via a quark-loop like in Fig. iLBb. 



Let's start computing the real emission diagram of Fig. |1.6| a in the kinematic regime 
^ and x,l — x not too small, where we expect scaling violations. In this regime, 
quarks can be considered massless, but in so doing the integration over the internal mo- 
menta k diverges when the transverse component k vanishes, i.e., when the internal quark 
and the outgoing gluon are emitted coUinearly (coUinear singularity). In our notations, all 
boldface transverse vectors have to be considered D — 2 dimensional vectors with euclidean 
components, such that fc^ > 0. 
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By means of a suitable regularization, by introducing for instance a (small) cutoff 
A, we end up with a finite result which has the interesting feature of containing a term 
proportional to InQ^/A^. The same behaviour is shown by the diagram in Fig. |1.6| b. 

This logarithmic factor comes from the ^-integration 

/ V^l^^ ' fcmin = A^ ; fcL. = Q^—:^ = (small const)Q^ (1.38) 

It is crucial to note that to the IR logarithmic coUinear singularity corresponds a InQ^ 
term arising from the hard k component of the emitted quark. This is obvious from a 
simple dimensional analysis, being the only dimensional scale which can accompany 
the cutoff A^ in the logarithm. Conversely, to each InQ^ is in fact associated a collinear 
logarithmic singularity so that one can investigate the structure of the latter for deter- 
mining the former. Here resides the intimate connection between collinear singularities 
and scaling violations. 

The presence of a large InQ^ at 0{as) should warn us as far as the convergence of the 
perturbative series is concerned. In fact, taking for as its running value as{Q'^) ~ 1/ InQ^, 
the NLO correction ~ as{Q'^)lnQ'^ ~ 1 may be large even for small coupling. This is 
exactly what happens at higher orders, where in the 0(a") contribution we find terms 
proportional to q;"(Q^) ln"((5^) ~ 1. Therefore the relative importance of all the sub- 
leading corrections may be the same, and in order to provide quantitative estimates, all 
of them ought to be taken into account through a suitable resummation procedure. 

The first step of a resummation procedure consists in rearranging the perturbative 
series according to a new hierarchy, the first term of which is formed by the "largest" 
contributions ~ [aglnQ^]" : n = 0, 1, ■ ■ ■ (leading logarithmic (LL) approximation), the 
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second one containing one power of ttg niore than the first one, i.e., ~ as[tts InQ^]" (next- 
to-leading logarithmic (NLL) approximation), the third one (NNLL) ~ ^^[as InQ^]" and 
so on. 

The second step requires the determination of all the LL terms in the perturbative 
series by identifying an iterating common structure — typically an integral kernel — 
representing the building block of the corresponding Feynman diagrams. 

In the third step one expresses the unknown resummed object by means of an equation 
— typically an integro-differential one — whose basic ingredient is the iterating structure 
previously mentioned. 

Finally, one has to solve the resumming equation as best as he or she (or the computer) 
can, to check its consistency, to try to fit the data and, possibly, restarting from the 
beginning in the next logarithmic approximation. 

In this section we are dealing with the real emission, 1-loop correction to SF at leading 
In level. For both |1.6| a and b processes, the logarithmic term stems from the integration 
over the internal quark transverse momentum (Eq. (|1.38|) ). The link with the coUinear 
singularity tells us that the coefficient of InQ^ is given by the residue of the pole. 

What about other real emission diagrams? Following the resummation program, we 
have to identify the diagrams furnishing the LL terms we are looking for. It turns out 
that, adopting axial gauges, i.e., including only physical transverse polarizations in the 
gluon propagators, only the diagrams L6a and b contributes in the LL approximation. 

It is convenient to introduce the Sudakov parametrization 

k = zp + zq + k , q := q + xp , x:=— — , (1.39) 

2p ■ q 

which are nothing but the components of k with respect to the light-like vectors p and 
q' plus the transverse component. In our notations, fc^ = 2zzp ■ q' — and k^ > 0. 
For vanishingly small values of k, the mass-shell constraints of the outgoing partons force 
^ ~ and z ~ X so that the only relevant variable (upon which the residue depends) is 
the momentum fraction x of k with respect to p. 



The residue of the coUinear singularities of Figs. |1.6| .a and b turns out to be respec- 
tively0 



^P'\z) = ^6'^P^\z) , (1.40a) 
^pbg(^) = ^P^^{z) , (1.40b) 

ZTT ZTT 



^^From now on, both for the sphtting functions and for the anomalous dimension coefficients we shall 
omit the 1-loop suffix (1) which will be always understood, unless explicitly specified. 
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where q denotes a generic quark or antiquark and the quark-to-quark and gluon-to-quark 
sphtting functions readQ 



1 + 2^ 

P^^{z) = Tn[z' + (1 - zf] 



;i.4i) 

:i.42) 



The partonic tensor W^/^u at first order in ag assumes then the form {p = ^p) 



k ^0 ^ 



;i.43) 



It is tempting to suppose that the hadronic tensor W^^{p, q) could be obtained by summing 
up the partonic ones q) weighted with the parton densities f^^\C)- Including the 

1/^ flux factor for compensate the fractional momentum of the incoming parton with 
respect to the proton, we try writing^ 



/ 



(a) 



:i.44) 

;i.45) 



where in the last equality we have changed the order of integration ^ u and introduced 
the shorthand notation for convolutions 



:i.46) 



in which a sum over partonic indices is eventually understood. 
A couple of remarks are in order: 



we are still in presence of the IR singularity A — which prevents us from inter- 
preting Eq. (|1.45|) clS cl reliable expression for PV^i,; 



^^The qq splitting function is not in its full form, having not yet considered the contribution of virtual 
corrections to diagram 1.5. 

^^The quark absorbing the photon carries a fraction of momentum x with respect to the proton and 
thus x/S,{^ l!) with respect to the incoming parton. 
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• we have used perturbative theory in the (strong couphng) IR region fc^ ~ —k^ 
without being allowed to. is not yet running, but can we expect to have obtained 
the right answer? 



The way to handle and to eliminate these drawbacks is easily explained. Let fip ^ A 
be a sufficiently hard mass scale, called factorization scale. We believe perturbation 
theory for |/c^| ~ fc^ ^ while we cannot say anything for < fj^p. However, the bare 
parton densities are completely unknown, so if we insert the non perturbative contribution 

^ fc^ < in / redefining the partonic distributions, we don't introduce any new 
uncertainty in Eq. (|1.45| ). Furthermore, in so doing the new partonic distributions absorb 
the coUinear singularities which do not affect W^^ anymore. In practice, we write, at 
0{as) accuracy. 



27r A2 



l + ^^ln^i^)®/=(l + gln^P)®(l + ?^ln^P 



27r A2 



:i.47) 



and define 



(a) 



(0 



;i.48) 



to be considered as an input parameter. 

By writing the W tensors in terms of (partonic) structure functions, Eq. ( |1.45| ) yields 



{i = 2,L). (1.49) 



Obviously, the RHS of the above equation cannot depend on the artificial parameter ^p. 
It is natural to define the Q^-dependent PDF asQ 

d^ 

b 

ending up with the desired factorization formula 



f'\z,Q') := f'\z,^.l) + ^In^ 5: r f P^'{^)f'\i,^^l) (1.50) 

^TT lip J z <, S 



:i.5i) 



16 



To be fair, up to now we can define the "dressed" PDF of Eq. ( 1.50 ) as well as the "bare" ones in 
Eq. ( 1.48| ) only for quarks (a = f, f). In the next section we will justify the extension to the gluonic case 
we are already using here. 
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Our complete ignorance of f^^\z,jj?p) prevents us from determining the PDF f^^\z,Q'^) 
at a given value of Q^. Nevertheless, their full Q^-dependence is contained into the 
logarithmic factor of Eq. ( |1.50| ). Taking the derivative with respect to InQ^ yields 



d 



ding- 



2^ 



The identification of the splitting functions P^^ in Eqs. ( |1.33|) and (|1.52|) is straightforward. 
We are not far from completing the bridge connecting OPE and improved parton model. 

Before embarking in virtual corrections and higher order diagrams, let's derive — for 
future reference — the solutions of the AP equations in this fixed-coupling situation. In 
moment space Eq. (|1.52|) reads 



ding- 



f,f,g). 



1.53) 



This set of 2Nf + 1 coupled differential equations can be diagonalized as usual by in- 
troducing the "non-singlet" , "plus" and "minus" components for PDF and anomalous 
dimensions 7^^. The solutions are then readily obtained: 



fL''\Q') = {QT''ci^^ a = NS,+,-. 



:i.54) 



The 1-loop virtual correction diagrams contributing with Ing^ (Fig. L7) are both UV 
and IR divergent. 







(a) 



(b) 



Figure 1.7: 1-loop virtual diagrams: (a) self-energy and (b) vertex correction. 



The UV divergences are regulated with the introduction of the running coupling con- 
stant whose argument has to be of the order of the virtuality of the probing particle (the 
photon). 
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The IR divergences, on the other hand, originate when the incoming quark emits a 
gluon with vanishing 4-momentum (soft singularity), i.e., z = 1, z = and A; = in the 
quark Sudakov variables, or when the quark emits a collinear gluon. A corresponding 
term proportional to 5(1 — z) has to be added to the quark-to-quark splitting function in 
Eq. dm]) 



y 



1.55) 



which reproduces the quark-to-quark splitting function in Tab. 1.1 



1.7.2 All order resummation of leading logarithms 

To complete the resummation program, we should determine the LL contributions to all 
order in ag. An accurate analysis of higher order diagrams shows that: 

• in axial gauges, the real emission diagrams contributing in LL approximation are 



the ladder type ones like in Fig. 1.8 



the phase space region generating large In corresponds, in the Sudakov variables 
of the internal particles 

ki = Zip + Ziq + ki , (1.56) 

to ordered transverse momenta 

> kl> kl> ■■■kl> X\ (1.57) 

in fact, the j-th internal parton plays the role of the virtual photon in regard to 
the lower part of the diagram, so that the arguments of Sec. |1.7.1| can be repeated 
showing that the integration for k'j^^ > is suppressed; 

the virtual contributions in LL approximation corresponds to vertex and self-energy 
corrections to the ladder diagrams; 

in addition to g ^ q and q ^ q splitting, we find g ^ g and q g processes 
involving collinear singularities; 

the mass-shell constraints of the outgoing partons forces the Zj variables to be 
ordered along the ladder: 

X ~ 2:1 < 2:2 < ■ ■ ■ < 2;„ < 1 , (1.58) 

while the Zi components in the ordered transverse momenta region ( p.. 57] ) are very 
small. 
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Figure 1.8: Ladder-type diagram contributing to DIS in the leading logarithmic approximation. 



The n-loop contribution to the partonic tensor is 



dk 



A2 fci Jx^ k 



^-as{kl) 



Zj 2lT Zj 



(1.59a) 
(1.59b) 



27r 



(1.59c) 



By summing over all n G N and integrating over ^ with the partonic densities just like in 
Eq. ( pri^ ) we obtain, for the SF, 



F,(x,Q2) = / de 5^AW(g^)5(™K7)/^^Ho 

a,ai n=0 ^ 



(« = 2,L), (1.60) 



where, from Eqs. and (f]^), F^^^^ = e\ and fJ'^^ = 0. 

The integral over transverse momenta can be performed by introducing the variable 



L(Q^ 



dk^ 



A2 



as (A; ) = — In 



1 «.(A2) 



(1.61) 



and gives 



n! 



(1.62) 
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By performing a Laplace transformation with respect to L, we obtain 

/•OO 1 

The integral over Zj variables is easily evaluated in Mellin space by defining 

^(^aia,n)_ j z""^ E^^^^'") (z) 

Jo 

— ST^ ^aia2 . . a„_ia„ a„a _ f-^nl ^l^ 

/ J lUJ lUJ lU) [ IU)\ 1 

a2---an 

SO that 



p. _ D(aia,n) r(a) 

a,ai n=0 

a,ai n=0 V / aia 



:i.63) 



:i.64) 



(1.65) 



In (cu, p)-space, the above results can be summarized in the following diagrammatic rules: 



^ = t 



,qq 



= y 



1 
P 



gq 




= T 



gg 



By inverting the Laplace transform we obtain 

.c+ioo 



2m 



(a) 



a,ai 



By splitting the exponential factor 
exp{7L} 



r«s(A2)i 






-y/bo 


r«s(A2)i 


■y/bo 


L«s(g2)j 








.as(PF). 





and redefining the bare parton densities 



■y/bo 



f(b) 



ab 



;i.66) 



(1.67) 



:i.68) 
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we get an expression for the SF free of singularities which is the cu-space analogue of 
Eq. ( [1.49| ). The moments of the Q^-dependent PDF are defined by 



L"s(Q2)J 



f/bo 



;i.69) 



ab 



in terms of which we finally obtain 



;i.7o) 



i.e., in x-space, Eq. (|1.51| ). 

The evolution equation for the PDF, obtained by differentiating Eq. (|1.69| ) with respect 
to InQ^ reads 



ding- 



(a = f,f,g). 



;i.7i) 



The main differences of the resummed formula (|1.71| ) with respect to the 1-loop one 
( |1.52| ) are that (i) it has been obtained without approximations apart from the LL re- 
strictions, (ii) the coupling is running and (iii) the evolution equation is extended to the 
gluon densities where 7^^ and 7^^ anomalous dimensions are present. 



Chapter 2 



Small-x hard processes 



The coming of high energy colhders is making it possible to investigate strong interaction 
processes in very pecuhar kinematic regimes, where the CM energy is much larger both 
than the hadronic masses and of the transferred momenta. 

From a theoretical point of view, this new class of phenomena is of great importance 
by testing the asymptotic behaviour of cross sections in the high energy limit on one side 
and to check QCD theory in a much wider kinematic domain on the other side. 

From the phenomenological point of view, the discovery of a marked rise of structure 
functions at HERA |[T^ has increased the interest of physicists for understanding the high 
energy features of QCD. What is actually observed (Fig. |2.1| ) is a surprising growth of 
the F2 structure function towards low values of x. Note that 1/x is nearly proportional 
to the CM energy squared of the virtual photon-proton system 

w^c^^n-x) ^ ^ (x<i). 

X X 

The growth is particularly marked for large values of Q^, but it is also evident down to 
very low values (Fig. |2.2| ). 

In this chapter, after a brief introduction to Regge theory, we expose the general 
framework of the perturbative treatment of high energy inclusive strong interactions. 

2.1 Regge behaviour 

Before a fundamental theory for strong interactions — such as QCD — was established, 
the study of scattering of hadronic particles relied basically on rather general assump- 
tions on the scattering matrix, such as Lorentz invariance, crossing symmetry, unitarity. 
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Figure 2.1: x-dependence of the F2 structure function at different values of Q^, showing the 
steep rise towards small-x for all values of . 



causality, analyticity, asymptotic states etc.. The subsequent development led essentially 
to what is known as Regge theory |jl3[, which determines the asymptotic behaviour of 



cross sections in the high energy limit regardless the strength of the coupling, i.e., inde- 
pendently of perturbation theory. 

In this section we outline some important results of Regge theory which will be useful 
for our study on semi-hard processes. Keeping in mind our interest in high-energy inclusive 
reactions, unitarity allows us to relate the total cross section for the scattering of two 
initial particles "1" and "2" to the elastic scattering amplitude of the same particles. 

a a' 

More precisely, by denoting with A4^^J{s,t) the scattering amplitude for the s-channel 
(s > 0, — s ^ t ^ 0) process 

(pi,ai) + (p2,a2) (p'i,a'i) + (^2,32) ; s := {pi + P2) , t := {pi - p\) , 
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Figure 2.2: Small-x rise of F2 for very low values of Q^. 

(aj labelling the discrete quantum numbers) the optical theorem states thatQ 

<r{pi,P2){s) = i '^mMllll{s,t = 0) . (2.1) 

Consider on the contrary the t-channel (t > 0, — t ^ s ^ 0) process 

{Pi, ai) + {-p[, a[) {-P2, 32) + {P2, a's) , 

where the antiparticle 82 of "2" is now outgoing and the antiparticle a[ of "1' " is incoming. 
By exploiting the crossing symmetry and the analyticity properties of the amplitude, we 
can say that the corresponding amplitude is given by the same function t), or 

better, by its analytic continuation from the s-channel to the t-channel region. 

In the t-channel CM frame [pi — pi = 0), the scattering angle 6 between pi and —p2 
is given by 

2s 

cos^ = l + — ; -l^cos^^l. (2.2) 

t 

By expressing the scattering amplitude by means of {cos9,t) instead of (s,t), it is easy 
^AU masses are supposed to be much smaller than the CM energy y/s and are thus neglected. 
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to perform the expansion in series of Legendre functions 



oo 



M{s,t) = M{cose,t) = 16n > (2/ + l)Mi{t)Pi{cose) (2.3) 



1=0 



which is called the partial wave series for Ai. 



We know from quantum mechanics that P/(cos^) is the emission pattern of a state 
whose total angular momentum (in the CM system) is /. Hence we interpret the partial 
wave amplitude Mi(t) as the contribution to the total amplitude corresponding to the l^^ 
angular momentum component. 

Coming back to the s-channel domain, | cos6\ > 1 loses its physical meaning. Nonethe- 
less, thanks to analyticity, the contribution from the l^^ partial wave amplitude is still 

167r(2/+ l)MKt)P«(l + y) . 
In the high energy limit s ^ |)f:| it holds 



t' 2'r2(/ + i)Vt, 

which indicates a power-like behaviour in s for the elastic scattering amplitude whose 
power equals the angular momentum of the state mediating the interaction. Moreover, 
the leading contribution to the amplitude will be given by the highest spin intermediate 
state allowed in the t-channel process. 

As far as the differential cross section is concerned, we have (Q^ = —t) 

da 1 |^(,,-g2)|2 



dg2 167r s 



2 



s''-' . (2.5) 



At this point we can anticipate that, in QCD, the leading contributions to cross sections 
will be given by gluon exchanges in the t-channel. In fact, s and u channel exchanges 
are suppressed by the large denominators 1/s and 1/m of the virtual particle propagators 
with respect to t-channel ones which involve 1/t {s ~ \u\ ^ |t|); the highest spin field in 
the QCD Lagrangian is the spin 1 gluon which, according to Eq. ( p. 51) is responsible of 
nearly constant differential cross sections. 

What we have sketched are nothing but rough estimates, which give us an idea about 
the "order of magnitude" of the phenomena we are investigating. It is possible to push 
further the partial wave analysis in a beautiful mathematical description where the am- 
plitudes Mi{t) are analytically continued to complex values of the angular momentum 
variable /. 
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The complex functions / t— > Mi{t) may present poles or branch cuts whose positions 
ai{t) are smooth functions of t. Having in mind the t-channel process (t > 0), we expect 
the amplitude to have a pole corresponding to the exchange of a physical particle of 
spin J and mass m, thus 3i : Q;j(m^) = J. It turns out that the functions t ^ C(i{t), 
which are called Regge trajectories, assume physical (semi)integer values when t equals 
the squared mass of a particle or of a resonance. For each set of quantum numbers (except 
spin!) there is a definite Regge trajectory along which hadronic particles carrying those 
quantum numbers lie. 

However also non-physical value of ai{t) are important: each Regge trajectory con- 
tributes to the scattering amplitude in the s-channel. The kind of contribution depends 
on the particular nature of the singularity associated to the trajectory, e.g., 

Miit) ~ _ — Mis,t) ~ .-W[ln.]^W-i . (2.6) 

as if a particle of "effective spin" a{t) were exchanged. We are not surprised to encounter 
logarithmic corrections to the naive power-like estimates: these are normal features in 
interacting theories. 

According to the optical theorem, we argue the total cross section to be expressed — 
apart from logarithmic corrections — as a sum of powers of s involving the intercepts of 
the various Regge trajectories: 



^c,s"'W"i ™ CMs"-(°)-\ (2.7) 



where aAf(O) is the largest among the intercepts. 

A careful analysis carried out from Froissart, based on general grounds and which 
assumes a small range force (as the strong interaction is) shows that unitarity limits the 
asymptotic behaviour of the scattering amplitude, setting an upper bound on the largest 
s-growth exponent: Q;Af(0) — 1 ^ 0, i.e., the total cross section of any reaction cannot 
grow faster than some power of In s. 

In this connection the phenomenological analysis seems to present a sort of puzzle: a 
Regge inspired fit like Eq. (|2.7|) with two terms {i = 1,2), as performed by Donnachie 
and Landshoff |]14|, describe very well the experimental data, with the same Reggeon 
intercepts a;R(0) — 1 = —0.45 and ap{0) — 1 = 0.08 for different processes (see Fig. P73| ). 
However, the latter clearly violates the Froissart bound. The answer to this apparent 
paradox has to be found in the kinematic domain so far investigated which is evidently 
outside the regime where unitarity effects become important. 

Anyway, there is evidence of a dominance of Regge-particles in the high energy region 
so far investigated. The lower Regge trajectory corresponds to the family of particles 
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5 



Figure 2.3: The Donnachie-Landshoff fit to (a) proton-proton and antiproton-proton and (b) 
to photon-proton total cross section. 



including the p, u, etc. mesons, and has odd parity and charge conjugation. This is the 
reason of the different coefficient ck in the pp and pp cross sections. 

The upper trajectory does not distinguish charge conjugated particles — it has the 
quantum number of the vacuum — and does not correspond to any particle till now 
observed. This elusive particle has been given the name of pomeron and is the responsible 
of the high energy behaviour of cross sections, at least before the unitarity regime, where 
pomeron self-interactions are no longer negligible. 

Is perturbative QCD able to describe the pomeron? and, in general, the high energy 
processes? These are the questions we wish to give an answer. 



2.2 Perturbative analysis of DIS in leading Inl/cc 
approximation 

We are ready to start the analysis of DIS in the high energy limit s ^ ^ A^, i.e., 
from Eq. ( A.13| ), xy <^ 1. In particular, we are interested in the kind of process where 



X <^ 1 and y <^1. 

The large virtuality ^ of the virtual photon should justify the use of pertur- 
bation theory, regardless of the value of s. However, in this regime where there are two 
very different hard scales, the QCD perturbative expansion is affected by large coefficients 
which, to order a" and for inclusive observables, are (9[(ln I^^/Q^)™^] ~ In™' 1/x : m ^ n. 
We have to face again a resummation problem, this time with respect to the effective 
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expansion parameter ctg In 1/a; which can be of the order of unity for small values of x 
even if ^ 1. The leading In 1/x (Lx) approximation consists in taking into account all 
the perturbative terms of order [a^ In the next-to-leading In 1/x (NLx) considers all 

contributions like as[cts In l/x]" and so on. 

In order to resum the leading logarithms of x a technique is usually adopted which is 
in certain aspects similar to the one employed in the resummation of the leading \nQ^. 
First of all, the coUinear factorization formula for SF has to be replaced with a cor- 
responding high energy (fe-dependent) factorization |]T5|]. The latter should take into 



account the larger phase space available and should generalize the former. The second 
step is to calculate the evolution (in x-space) of the "parton density" factor which obey an 
equation obtained by Balitskii, Fadin, Kuraev and Lipatov (BFKL) By performing 
the transverse-space integration of the parton densities with the partonic cross section, 
which can be explicitly evaluated, one obtains a factorized expression for the high energy 
(small-x) SF where both partonic SF and PDF are resummed. 



2.2.1 High energy factorization 



As we noticed in Sec. ^.l| and we shall show in more detail in Sec. the terms contribut- 
ing to the Lx approximation arise from gluon exchanges in the t-channel. Since gluons 
can couple to the (virtual) photon only via quarks, we are led to consider diagrams where 
single gluon (Fig. |2^a) and multiple gluon (Fig. |2]lb) exchanges are possible. 





(a) 



{b) 



Figure 2.4: Single (a) and multiple (b) gluon exchanges in the high energy limit of photon- 
proton total cross section. 



If we work in axial gauges, the multiple gluon exchange diagrams are suppressed 
by powers of as(Q^), because the gluons' endpoints have a very large relative velocity. 
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the upper one being nearly coUinear to pi and the lower one nearly coUinear to p2 (see 
Eq. O)). 

Single t-channel gluon exchanges in Lx approximation are then factorizable as we will 



briefly explain in the following. According to Eq. ( [A.4|) , the contribution to the hadronic 
tensor from Fig. |2.4|a can be written 



H^i- ab ab 

AZ{q,k)G»^{p2,k) 



(2.8) 



47r J (27r)4 

(summation over colour indices a, b is understood) where 

ab P a b 

Ki3{(l^k):= / d$(p3,P4) M'^{q,k,p3,p4)M''^{q,k,p3,p4y 
denotes the lowest order qq contribution to the 7*g* 7*g* absorptive part (Fig. p75|) 

ab 

integrated over the final particle phase space and G'"^(p2, k) is the full pg* pg absorptive 
part, including the gluon propagator. 



(2.9) 




+ 






Figure 2.5: Feynman diagrams for photon-gluon coupling at lowest order in ag- 



in the high energy, small-y regime, the differential cross section is essentially given 
by the F2 structure function, and can be obtained by the eikonal approximation at the 
electron-photon vertex. The eikonal approximation can be applied when the components 
of the transferred momentum q are much smaller than that of the external momenta pi and 
p[; it consists in the replacement of the helicity conserving vertex Ucr{p'i)'y'^Ucr{pi) — 2pi 
and neglecting the non conserving one, i.e., to approximate the leptonic tensor ([A.3|) 



{pi, q) ~ 2pi^2piv . 



After contraction with the hadronic tensor we obtain 



2pi^2pi, W^' = -Q'Fi + 2s 
and hence, using Eq. ( |A.13| ), 
F2 ~ 



y/2 2 ^ 2s ^ 
1L^F2 ~ —F2 

y y 



{y « 1) 



A'^h ab 1 ab 

AZ^{q,k)-G-^ip2,k) 



(2.10) 



f2.111 



(2.12) 



39 



By adopting the Sudakov parametrization 

q = ypi + yp2 + q (2.i3a) 

k = zpi + zp2 + k (2.13b) 

the high energy hmit constrains y <^ y <^ 1 and it turns out that the dominant integration 
region in the k variable is the one with fixed k^^ z <^ z and |A;^| ^ A;^ = 0{yzs) = 

The next step is to show that only a single gluon polarization contributes in the 



diagram of Fig. p^ a. Since the amplitude A^^ =: Ai^^T involves a colourless object (the 

ab 

photon), the tensor A'^^ =: A^^^^d""^ satisfies the Ward identities of electrodynamics 

= fc^A^^ = (2.14) 
and thus the following decomposition holds: 

Ao.p: = {^Pi,){^Pi.)A>:^^^{q,k) (2.15a) 

=^1 9ap] - A2JIT ka -Qa ] [ kfS -Qp (2.15b) 



P V " qk^"J V ^ qk 

where the dimensionless Lorentz invariant functions Ai depends on Q^, fe^, 2qk ~ yzs, 
2{ypi)k = yzs or, equivalently, on the rescaled variables x/ z and Q"^ /k"^. 

In the kinematic regime mentioned above, the amplitudes Ai and A2 satisfy 



• z' 

for fixed values of k? j zs (they are exactly equal when A;^ = in order to cancel the 
spurious pole at A;^ = in Eq. ( |2.15|) ). This is due to the fact that the corresponding 
diagrams Fig. ^.5| involve quark (spin 1/2) exchanges, vanishing like ~ x for increasing 
energy (see Eq. (|2.5|) ). Therefore, the dominant contribution in Eq. ( p.l2|) is obtained for 



Ai~A2~0(-) 



fixed values oi xj z = 0{\) and hence z ^ 0. On the other hand, when z — 0, A;^ ~ — fc^ 
and k'^ /zs is fixed, the polarization tensor in front of A2 can be replaced with 



^ fe^ PiaPW 1 



zs 



l + 0{z}^) 



s 



(2.16) 



which clearly shows an 1/z enhancement factor with respect to the Ai one and a "selec- 
tion" of the gluon polarization along p2- 

By using Eqs. (|2.15|) and (|2.16|) , Eq. (|2.12|) can be rewritten as 

ah 
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By defining the off-shell partonic cross section 



^^i-.S) ■= -M-S) - (2.18) 



-z z ^z k^ k 

and the unintegrated gluon density 

dk^ k^ 



T[z, k):= y ^ -^PicPif^ E ^ k) (2.19) 



the high energy, fc-dependent, factorization formula is estabhshed 

'Uz - 



F2{x,Q')= [ - [ d'ka,{-,%)n^,k). (2.20) 

Jo ^ J Z K 



Eq. (|2.18| ) shows that the hard partonic cross section (3"2 is given by the diagram in 
Fig. ( |2.5| ) where the soft {z 0) off-shell gluon is coupled to an external fast quark (or 
gluon) with an eikonal vertex It is important to note that (T2 is a gauge invariant 
object, despite the off-shellness of the "incoming" virtual gluon k. The reason for this is 
that the eikonal coupling induces physical polarization for the incoming gluon. 

The formula ( p.20| ) can be written in diagonal form in Mellin space. With the defini- 
tions of Eqs. (|C.19|) we get 

= I d'k a^{^)T^{k) (2.21) 

or, equivalently, 

^2,^(7) =^-(7) ^-(7) • (2.22) 
2.2.2 Relation with the coUinear factorization 



The high energy factorization formula (|2.20| ) is a generalization of the collinear one ( p..5| 



which holds not only for s ~ ^ but also in the semi-hard regime s ^ ^ A^. 
This can be easily seen if we remember that, for s ~ Q^, the leading contributions to 
the cross section arise from emission of partons whose transverse momenta are strongly 
ordered. In particular, the integration over k is important only for k^ <^ Q^, where the 
hard cross section assume the typical form of collinear emission 

^(^'fJ)- E^a^l^f^^^^'(^) ik'«Q')- (2.23) 

a£{partons} 

In tu-space this means that 

= (E^-) f «s7S^ln|J-^.(fc) . (2.24) 
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Since 2^ is just the partonic structure function for the quark singlet density, we 

identify 

fi'^KQ') = f d'k a^2NpJ\n^r^{k) . (2.25) 

Taking the derivative with respect to In yields 

d 



ding^ 



HQ') =o^.2NpJ j d'k r^{k) (2.26) 



once we have made the further identification 

flf\Q') = f d'kJ^Uk). (2.27) 

In conclusion, J?-L(fc) represents the unintegrated gluon density in fc-space which is related 
to the usual gluon PDF by Eq. ( |2.27| ). The former gives the most important contribution 
to the quark singlet density in the framework of high energy factorization. 

However, when s ^ we cannot exploit the collinear picture to describe the 
evolution of the singlet density, since a does not necessarily contain a leading InQ^. 
Rather, a has to be considered the 0{as) coefficient of the gluon density contributing to 
F2. The fc-convolution of a with the unintegrated gluon density provides the resummation 
of the large In 1/x in the coefficient function, as we shall see in Sec. p. 3. 3 . 



2.3 Resummation of Inl/cc in the leading approxi- 
mation 

The high energy factorization formula (|2.2(]|) shows that the leading high energy contri- 
bution to semi-hard DIS is governed by the unintegrated gluon density for small values of 
the gluon momentum fraction z. As one can easily realize, it is not possible to determine 
perturbatively the gluon density, since long distance effects are unavoidable when dealing 
with hadrons. Nevertheless, in a particular regime (to be specified soon), the perturbative 
analysis allows us to extract information about the dependence of the unintegrated gluon 
density on the "hard" variables z and A; involved in the hard gluon-photon vertex. This 
is done by means of an evolution equation in z-space obtained more than twenty years 
ago by the russian school. 

In this section we will go over the road that led to the BFKL equation |T^. The 



latter resums all the leading logarithmic (Ins) coefficients of the ag perturbative series 
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for inclusive semi-hard processes. The following presentation will be useful in order to 
identify the basic physical quantities (in connection with the DGLAP ones) and to fix the 
framework for the NLx formulation of the BFKL equation (Cap. ^) and for its subsequent 
improvement (Chap. 

Following the idea of factorization of non perturbative effects, we can guess that the 
probability to find a gluon with momentum fraction z with respect to the parent proton 
and transverse momentum h should be given by a convolution of bare partonic densities 
f^^^ with a perturbative calculable function Q^^^ which embodies all kinds of intermediate 
processes between the bare parton a and the hard gluon. At high energies, only gluon 
exchanges in the t-channel contribute in the Lx approximation. Therefore we assume a 
sort of factorized expression for the unintegrated gluon density of the form 

T^(k) ^ J d^ko fi^Kko) g^k, ko) (2.28) 

which stresses again the relevance of the transverse degrees of freedom. Here f!f\kQ) 
denotes the Mellin transform of the probability density of finding a gluon with transverse 
momentum fcg inside the proton; it is to be considered as an unknown input function to 
be eventually specified in particular phenomenological applications. 

The function Quj{k, ko) is the basic object we wish to investigate; it will be extensively 
studied throughout this thesis. 



2.3.1 Leading \nl/x BFKL equation 

The function Qu){k,kQ) represents the set of gluon exchanges contributing in Lx approx- 
imation to high energy scattering of two strong interacting objects. Also high energy 
scattering of colourless particles is governed by gluon exchanges: in this case the coupling 
to gluons occurs via quark or EW boson intermediate states. 

In order to determine Q^)-, the simplest situation we can consider is high energy scat- 
tering of two partons a and b. At Born level, the cross section is dominated by one gluon 
exchange in the t-channel (Fig. |2.6| ) corresponding to the amplitude 

Mf_^^ = -i2sT:,,-^Vl,^ , s = {p,+p2f , t={p^-p\f (2.29) 

where summation over repeated colour indices is understood and the parton-gluon vertices 
are 

K'. ■= 9 5a,,A3 r''" (a = g = gluon) , (2.30a) 
r:,, := -ig 6x,x, Ta'a (a = q = quark) . (2.30b) 
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Figure 2.6: The lowest order Feynman diagram contributing to leading Lx order. 

Since gluon exchanges do not affect the ffavour of the quarks and the polarization/hehcity 
Aa' = -^a is conserved, only colour degrees of freedom and kinematics have to be taken 
into account. Note also that the longitudinal variables corresponding to the Sudakov 
decomposition k = zpi + zp2 + k are small in the high energy regime: z = \z\ = \t\/ s 1. 




Figure 2.7: The two particle final states virtual correction diagrams contributing at leading Lx 
order. Iteration of gluon exchange in the s-channel (a) and in the u-channel (b). 

The 1-loop virtual corrections to the elastic amplitude give rise to a Ins coefficient. 
In covariant gauges, only the "box" and "crossed" diagrams of Fig. contributes with 
logarithms of s. These amplitudes can be projected in the subspaces of the irreducible 
representations of 8 8 where the couple of exchanged gluon lives. In other words, a 
given amplitude can be considered as a sum of terms representing the exchange in the 
t-channel of colour structures belonging to the various 8 (E) 8 irreducible representations. 
When summing the amphtudes ^^2-^2 A^2^2 of Fig- ll]^^ ^^"^ ^ respectively, which 
kinematically corresponds to the replacement s ^ u, only the octet representation of the 
antisymmetric part (8 ® 8) a = 1 © 8 © 10 © 10 interferes constructively. For this reason 
the colour factors can be well represented by the parton-gluon vertices of Eqs. (|2.30|) . 

The octet part of reads 

MS::l = g's'N,r:,j{s,t)ri,, (2.31) 
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where the loop integral / yields^ 

The logarithmic factor stems from integration over the longitudinal variables of k. In the 
massless theory the scale of the energy sq is of the order of the transverse momentum 
squared of the outgoing particles. In the Lx approximation, however, the choice of Sq is 
irrelevant, because with an other choice Sq the difference In s/sq — In s/sg = Insg/so is 
independent on s and therefore subleading. 

The coefficient Q{—t) of the logarithmic term is an IR divergent two-dimensional 
integral. The IR divergence will cancel when considering real emission correction to the 
Born amplitude (see App. |C.1|) . 



We note that among the subleading 1-loop corrections we have neglected there are the 
self-energy and vertex-correction diagrams which determine the running of the coupling 
constant. Accordingly, in a Lx treatment of the radiative corrections, as must be regarded 
as fixed. 

By collecting Eqs. ( |2.29| ), ( |2.31| ) and ( |2.32| ) the full 1-loop expression for the elastic 
amphtude results 

[M^'^+M^%^, = _,2sT:,, ^ + ^M-t)Hs/so) ^.^^ ^2.33) 

t 

Eq. (|2.33| ) is intriguing, since the logarithmic correction looks like the ag truncation of 
{s / So)"^^^~'^^ which would correspond to Regge behaviour. The 2-loop correction has been 
computed ||T6| and confirms this assumption. In conclusion, we make the ansatz that, in 
Lx approximation, the amplitude for elastic scattering with the gluon quantum numbers 
in the t-channel is 

M2^2 = -i 2s T:,, ^-1^^ r;,, (2.34) 

which states the reggeization of the gluon. 

Upon colour averaging and by using the phase space measure for two particle final 
states (including the fiux factor (2s) 

d0(^) = ^ldfc, (2.35) 



^The 1-loop Regge-gluon trajectory cisfi^*'-' is mostly known in the literature as o;^^-'. We prefer to 
adopt a different notation in order to avoid confusion with the uj variable and because of its connection 
with the leading BFKL kernel K'-^\ In this chapter the suffix will be understood. 
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Eq. (|2.34| ) yields the two-particle final state contribution to the differential cross section 

da3b^2 = h<^^\k) (-) h'^^\k) dk (2.36) 



in terms of the Born impact factors 




Figure 2.8: The Lipatov effective vertex for the emission of a real gluon: the blob represents 
gluon emission both from the exchanged gluon and from the external particles. 



Let's now consider the 1-loop order real emission corrections. At Lx level, the high- 
energy kinematics favours the emission of an additional gluon (Fig. |2.8| ) of momentum qi 
and polarization e in the central region qt ^ Qi ^ \qi\ ■ The corresponding amplitude is 

A<(°Z3 = -i2.r- ir^^^ir-,, (2.38) 

T^,.,(fci,A;2):= ig T'^'' Uli)J'ik„ h) , (2.39) 
where = — k2 , U = kf ~ —k] and 

J"(fc,. fc.) = + + K (5:* + 4) + pS + 4) (2.40) 

is the Lipatov effective vertex. 

The 1-loop virtual corrections contributing to the three-particle final states in Lx 
approximation can be computed, and their final effect — according to gluon reggeization 
— consists simply in the replacement of the propagators 

- - , Si = {q, + q,^,y {^ = l,2). (2.41) 

U U \so/ 

in the amplitude ( [2.38| ). 

After polarization and colour averaging, by using the three-body phase space measure 
{zi being the momentum fraction of ki with respect to pi) 

1 d..d<,.d<=.. 
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and the properties of the J vertex 



q'Uh, h)=0 =^ *'x'W^ = = 4^ (2-43) 

A ^ 

we obtain the three-particle final state contribution to the differential cross section 

da3b^3 = /^a(fci) - ^ - h,ik2)—dk^dk2 (2.44) 

whose leading logarithmic character follows from the zi integration with infrared boundary 
zi > sq/s (pure phase space would yield ^ < 2:1 < 1 — j-). 

To summarize the results till now obtained 

• the huge part of the cross section arise from t-channel small virtualities gluon ex- 
changes; 

• the longitudinal phase space grows logarithmically with s and the differential cross 
section is almost independent of the longitudinal variables (apart in the proximity 
of their upper boundary]^) whose integration yields just the In s term; 

• most of the longitudinal phase space corresponds to small values of the Zj Sudakov 
variables; 

• the mass-shell conditions of the outgoing particles causes the transverse momenta 
of the exchanged (reggeized) gluons to be the only relevant degrees of freedom for 
the amplitudes, furthermore kf ~ — fc^. 

Going to higher order in as is then straightforward. 

The scattering amplitude with 2 + n particles in the final state receives its Lx con- 
tribution in the phase space region where the outgoing particles are strongly ordered in 
rapidity. At tree level, this corresponds to (half) a gluon ladder (Fig. |2.9|) in multi-Regge 
kinematics (MRK): 

ki = ZiPi + ZiP2 + ki {i = 1, . . . ,n + 1) (2.45a) 

1 = 20 > ^1 > ■ ■ ■ > Zn+i = (2.45b) 

s 

— = ^1 < ^2 < ■ • • < Zn+l < 1 (2.45c) 
s 

where the squared transferred momenta = kf ^ —kf are of the same order and much 
smaller than s. Because of the mass-shell constraints and total momentum conservation, 



^I.e., for zi —> 1 where the longitudmal components of qi get large. 
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Figure 2.9: The n + 2 final particle gluon ladder contributing at Lx to high energy scattering. 
The thick wavy line represents reggeized gluon exchanges. 



we can take {zi : i = 1, . . . ,n} and {kj : j = 1, ■ ■ ■ , n + 1} as independent variables; 
accordingly the phase space measure reads 

n J n+1 

d</,(2+") = ^^—^ TT — TT dfe. • (2.46) 



The real emission amplitude is nothing but the generalization of Eq. (|2.38|) , where 
outgoing gluons are coupled to internal lines by the effective vertex ( 2.401 ). Finally, the 
virtual corrections provides the reggeization of the exchanged gluons with the replacement 
of the propagators as in Eq. ( |2.41|) , yielding [p!6tp!7[ 



M2^2+n = -i2.r- iElll^T^^ 



[S2/S0) 



with the shorthand notation Vti = tj). 

Having so obtained the 2 + n particle differential cross section 



T-^^^J-^I^ (2.47) 



n+l 



TT 



as\" (si/so)^=2^^ ■ ■ ■ (s„+i/so)^=2^-+^ fr dz,. 



n+l 



1=1 ' j=i 



(2.48) 



there remains to sum up over all values of n from to infinity. This is easily done in 
Mellin space with respect to the uj variable conjugated to s by defining 



cr 



ab •" 



-d. (s\-^ 



(2.49) 
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(2.51) 



By using the multi-Regge kinematic relations 

s.^'-^ql (z = l,...,n + l) , , . = ^i^^^, (2.50) 

we can perform the change of variables {s,Zi} —>■ {si} such that 

ds dzi dzn dsi dsn+i 

S Z\ Zn Si Sn+l 

and, after integration over the sub-energies Si G [sq, oo[, we obtain 

— ^ = /i3(/c)/ib(fco) (2.52) 

E°°/"j, J, 1 as 1 tts 1 
/ dfcs ■ • • dkn — 2 2 ' 

where k = —ki and fco = kn+i are the transverse momenta of the scattered partons a' 
and b'. The sum on the RHS is just the series expansion of the operatoif] 

^ ^ +-^^^^^-^ + --- (2-53) 



^ _ /C{v) _ /c(R) cu-/C(v) c^-/C(v) cu-ZCW 
where the real and virtual BFKL kernels are 

/C^^) (fc, fc') := «s2^](fc2)(52(fc - k') . (2.54b) 

It is evident that the distribution Quj{k, fco), determining the cross section apart from 
the impact factors and h^, represents all exchanges and emissions stemming from the 
two gluons that couples directly to the incoming partons. Therefore, we identify Qi^{k, fco) 
as the function introduced in Eq. ( p.28| ) and representing the evolution of the unintegrated 
gluon density at high energy. 

Finally, we can invert the Mellin transform ( p.49| ) to s-space and represent the Lx 
partonic differential cross section as 



d(Tab(s 



dfc dfc 



- = fP: (-] K{k)g^{k,k,)K{ko) . (2.55) 
J 2ni \soJ 



The parton-parton cross section just analized, due to the chromatic nature of the scatter- 
ing particles, suffers severe Coulomb singularities, since the impact factors h(k) ~ (fc^)~^ 
prevent the definition of a total cross section for coloured objects. 



We ask the reader to forgive us for the "naive" notation! 
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On the other hand, by high energy factorization, Eq. ( p.55| ) is vahd as well for scattering 



of colourless particles, provided the corresponding impact factors are employed. For 
instance, in the case of virtual photons, the ensuing impact factors — depending on the 
virtuality of the photon — assume the form 

= h^mk) = (2.56) 

which is non singular for fe^ — and hence provide, upon k and feo integration in 
Eq. (|2.55|) , a finite total cross section (see Sec. p.2.3|) . 



The identity 

[a;-(/C(^) + /C(^))]e;^ = l (2.57) 
suggests that Q^^ is the solution of the Green function equation 

ujQu.{k, fco) = 5\k -ko)+ J dk' IC{k, k')g^{k\ ko) (2.58) 

where IC = /C(^) + /C^^) is the Lx BFKL kernel. 

The gluon Green's function (GGF) Qi^ can be determined once a complete set of eigen- 
functions for the integral operator K, is found. In App. y it is shown that the functions 



Um{k) = {ky-^ —j= , 7 = i + iz/ , z/GR , mGZ (2.59) 



form a complete set of eigenfunctions for /C whose eigenvalue is asX(7? '"^) where 

X(7, m) = 2V'(1) - ^(7 + ^) - ^(1 - 7 + ^) • (2.60) 

The finite expression of ^(7, m) confirms the cancellation of the IR divergencies of the 
real and virtual kernels. 

The GGF can be expressed by the spectral representation 

rfi.i.\ ^ Z"^^'"" d7 UUk)f;,m{ko) 

g^{k,ko)=y — 2.61 

which, by deforming the contour of integration to the left for k^ > k^, can be evaluated 
by the residue method. 



50 



2.3.2 Resummation of the anomalous dimension 

The full determination of the Lx GGF permits the high energy resummation of both the 



gluon anomalous dimension |]T6[ and of the coefficient function [|T5| (at Lx level). 

First of all, let's note that the unintegrated gluon density ( |2.28| ) satisfies the BFKL 
equation 

uJ'^ik) = uj:F^^\k) + jdk /C(fe,feo)^i°nfco) (2.62) 

where an uj has been factorized in the inhomogeneous term ujj-'!;j'\k) = /i,^-*(fc) for later 
convenience. 

If we adopt an azimuthally symmetric input function J-'!;j^\ then only the m = 
subset of eigenfunctions f^^m contributes to the expansion of JF^,'^'' and J-'uj. In this case 
the spectral representation corresponds to the Mellin transformation with respect to k^ 
in 7-space (App. p.2|) , where Eq. ( p.62|) reduces to the simple algebraic equation 

By inverting Eq. ( p^.22| ) in favour of the variable, we get 

Since we are interested to ^ A^, we displace the contour of integration to the left 
so as to enclose all the singularities whose real part is less than 1/2. Provided a smooth 
input function t!^^ has been adopted, the singularities in the integrand of Eq. ( p. 641) can 
arise either from the hard cross section ^"^(7) or from a vanishing denominator. 

(7(^(7) has poles both at 7 = — corresponding to its collinear behaviour for Q"^ ^ k^ 
— and 7 = 1 (Q^ <^ fc^) |15|. We explicitly show the former by writing 



^.(7) =: -hM ■ (2.65) 

7 



In the denominator of Eq. ( 2.64|) we have zeroes corresponding to spectral points of 



Qu}, i-e., when lj = asxil)- The symmetry, hermiticity and boundedness properties of 
/C causes xil) to assume real values only for 7 — 1/2 G I (the imaginary axis) where 
Xil) ^ Xm ■■= x(l/2) = 41n2 and for 7 G M. 

liuj > otsXm then the rightmost singularity 7l in the half plane Ke 7 < 1/2 — governing 
the high-Q^ behaviour of F2 — is a simple pole determined by the conditions 

^ = ^sX{li) , 0<7l<^ , 7l = 7l(— )• (2.66) 
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The other singularities to the left of 7l give terms which are suppressed as (almost integer) 
powers of (higher twist) and will not be considered in the following^. Anyhow, they 
may be important in moderate and intermediate-Q^ regimes. 

By applying the residue theorem to Eq. ( p.64| ) yields 
In order to extract the anomalous dimension, we note that in the as — > limit 

7l-|x'(7l)I-i , :f^S\i.) ^ :f^S\Q) , /^.(tO - /^.(o) (2.68) 

UJ 

and both F!f^{Q) and /ia;(0) are finite^. 

By comparing Eq. (|2.67|) with Eq. ( p..54| ), which gives — up to a constant coefficient 
function factor — the dependence of the SF in the fixed-Os formulation, we identify 
7L(as/i^) as the Lx gluon anomalous dimension. 7l resums all the leading as/oj terms of 
the gluon anomalous dimension 

^ U) ^-^^ \ijJ ) ijJ \ijJ J \ijJ ) 

n=0 

or, in x-space, the as In 1/x terms of the gluon splitting function 

^P^^(x) = ^Vc„ asln- . (2.70) 

n=0 ^ ^ 

Outside the anomalous dimension regime, i.e. for uj < asXm, the function 7l develop 
a branch-cut singularity at as/^ = 1/Xm corresponding to local non-invertibility of the 
eigenvalue function: x'{1l = 1/2) = 0. 



2.3.3 Resummation of the coefficient function 

By applying the same complex integration techniques to the integrated gluon density 
whose Mellin transform is 

clO' fQ^V^ ^«'d,^„(,).^.i^52) . (2.71) 



Q- UV J ' 7 ^ i_ 



^The eigenvalue functions xil, corresponding to conformal spin m ^ 0, which has been neglected 
in the present analysis, would contribute with subleading poles. In fact, the position of the rightmost 
pole of x(7j™) in the half-plane 3?e7 < 1/2 is 7 = —\m\/2 

^Even if h^ij) is singular for 7 — > in massless quark emission processes, nevertheless it is 0(q;s), 
and in the ckg ^ limit the singularity arising from the implicit as dependence of 7l is cancelled by the 
explicit as dependence of hi^. 
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it is easy to see that, in the anomalous dimension regime uj > asXm, 

A"i^'^-{%yS^r ''''' 

The expression ( ^.671 ) is then consistent with the standard QCD factorization theorem 
in which the remaining factor /i(^(7L(iSs/^)) represents the Lx gluon coefficient func- 
tion [|l^ for F2. The resummed effect is incorporated through the ots/oJ dependence of 7l 
of Eq. ( p^.69| ) and the 7-dependence of h^j. 



2.4 High energy behaviour of the structure functions 



The inverse Melhn transformation in cj-space of Eq. ( |2.67| ) allows us to determine the 



small-x behaviour of F2 according to the BFKL resummation of the leading Inx. 

Regarded as a function of a;, -F2,a;(Q^) in Eq- ( P-67|) presents its rightmost singularity 



at 

u = u;p(as) := "sX™ = 4 In 2 (2.73) 

both for the presence of the branch-cut of 7l(q;s/i^) and for the vanishing of x'(7l) the 
denominator. By assuming F^^^ and h not too strong dependent on cu and 7 around cUf 
and 7l = 1/2 respectively, we have 

"^^'^^^ ^ - lioo ^ (2.74) 

where A = lu - lu^, x'L = x"(V2) = 28C(3), = 4^^? (l/2)/i^,(l/2) and 



A ' 7L(^)-i-aAV2 , t = ln%. (2.75) 
The integration in the complex A-plane can be performed explicitly and gives 

^,„(g'^).^,^__^^j e.p|-^-^|, (2,76) 

Several facts can be learned from the last expression: 

• at high energy the leading anomalous dimension 7l saturates at the value 1/2; 

• correspondingly the coefficient function h^{'~fi,) undergoes a large enhancement with 
respect to its "coUinear" value where 7l ^ since hi^{l/2) is usually much larger 
than h^{0) because of the singularity of h^^ at ■y = 1; 
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apart from small logarithmic terms, the F2 structure function — and hence the total 
cross section — presents a power-like rise towards small values of x. 



As explained in Sec. the last point is in disagreement with the general results of 
Regge theory: the Froissart bound cxtot ^ const(lns)^ is violated. The reason is that the 
s-channel unitarity constraints for scattering amplitudes are not fulfilled in Lx approxi- 
mation. Nevertheless, a power-like growth in x, or equivalently, in s, has been observed 
for the structure functions. The agreement of the Lx prediction ( |2.76| ) is only qualita- 
tively (power-like growth), the small-x exponent Uf = as41n2 being too large to fit the 
data, but the order of magnitude is OK. as is not yet running, NLa; corrections have to 
be taken into account... 



Does that mean that we are on the right track? 



Chapter 3 



High energy processes in 

next-to- leading Inl/x approximation 



The analysis of high energy scattering in QCD by means of the high energy factorization 
and of the BFKL approach has clarified the mechanisms of logarithmic enhancement of 
the cross sections and has provided a powerful tool for the resummation of the large 
coefficients of the perturbative series in the Lx approximation. Nevertheless, we have 
remarked that the BFKL equation, at Lx accuracy, is not satisfactory both from the 
theoretical and from the phenomenological point of view. 



Concerning the last point, the small-x rise of F2 shown in Figs. |2.1| and is compatible 
with a power-like growth 

F2(x,Q2)~C(Q2)x-"(«'). (3.1) 
Figs. 3J- show the measured value of A((5^). Note that X{Q^) is larger than the pomeron 



intercept up ~ 0.08 but the former joins smootly the latter for < 1 GeV^. 

The small- a; growth exponent appear to be much smaller than the Lx BFKL result 
'^p(tts(Q^)) — 2.77ds{Q'^) (where the running coupling replaces the fixed one). In fact, 
for a typical value of ~ 30 GeV^ so that as — 0.2, the experimental exponent is 
A(Q2) ~ 0.3, definitely smaller than ujp{0.2) ~ 0.55. 

Furthermore, A(Q^) increases with Q^, while ci;p(q:s(Q^)) decreases. The DGLAP 
equations, on the contrary, provide an accurate description of the rise of A(Q^). 

The rise of A((5^) is a consequence of scaling violations, and one can roughly under- 
stand the underlying mechanism by observing that 

7^^ + 75^ bofor^^l 
2iV/75s + [ ^ for ^ 1 ' 
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Figure 3.1: Small-x growth exponent for the F2 structure function versus Q^. 

so that the Q^-evolution enrich the low-a: (a; — > 0) content of the PDF at the expenses of 
their large-x (tu +00) one. 

Indeed the HERA small-x F2 data are very well described by NLL DGLAP fits (see 
Figs. p.l|J2.2| ), provided suitable parton distributions are taken at the beginning of the Q^- 
evolution. In particular, nearly constant (in x) PDF at very small Qq ^ 0.35 GeV^ are 



able to take into account the observed behaviour for the SF [|18| . As an example, consider 
the singlet density — governing the high energy regime — f^~^\x, t), t := InQ^/A^, whose 
small-x behaviour is controlled by the singular part of the gluon anomalous dimension 

+ ^ Os _ J_ 

uj bout 

The DGLAP equation, for a flat starting density 



(3.3) 



<•{+) 



'{x, to) = const , 



has the simple solution 



J Ll. 



'j^~^\t) = const X exp 
The x-distribution is then given by 



/W(x,t) 



27ri 



buj to 



(3.4a) 
(3.4b) 

(3.5) 
(3.6) 
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which, by a saddle point evaluation, gives the asymptotic solution 



/«(x,*)-o,.tI!;^exp{^ll„llni}. (3.7) 

Eq. ( |3.7| ) shows essentially the t-evolution obeyed by the SF. The x-dependence is however 
incomplete, because of the crude approximation (|3.4b|) . 



Since the DGLAP approach is based on a purely perturbative ground, we may wonder 

2 




whether we are allowed to start the evolution at so a small value of as a fraction of 



GeV^. The same evolution for F2 is obtained with power- like behaved input PDF [|19| 

f{x, Q^) oc x-^-^^ for ~ 4 GeV^ . (3.8) 

This is a reasonable starting point for apply perturbation theory. 

Can perturbative QCD explain the power-like shape ( p.8|) for the PDF? In principle 
it should, since as(4 GeV^) ~ 0.3 should be enough small. Of course, the high energy 
regime involved at small-x requires the resummation of the logarithms of 1/x so that 
BFKL appears to be the most natural approach. On the other hand, we should not 
expect BFKL to give a reasonable answer at higher values of , since it ignores the In 
resummation. 

The natural field of applicability of BFKL evolution is the class of processes with 
two hard scales. We can mention 7*7* scattering ||2^, which occurs, e.g., in double DIS 
(Fig. |3.2| a), the two hard scales being the (large) virtualities of the photons. Another 
two-scales process being studied concerns the so called forward jet (Fig. ^]^b), where the 
two hard scales are the photon virtuality and the jet transverse energy. By selecting large 
rapidity differences between the two observed objects we demand a large BFKL evolution, 
since we allow strong ordering in longitudinal space along the ladder. On the contrary, 
we require the two extreme transverse scales to be as close as possible, so as to suppress 
DGLAP evolution. 

While 7*7* measurements are not entering the small-x domain yet [^, the forward 



jet analysis is affected by rather large uncertainties mainly concerning hadronization phe- 
nomena and limited statistics, giving consistent but not very significant agreement with 
BFKL resummation 



We will not deal with unitarity corrections, which are hopefully not too important for 
the perturbative indices that we shall determine. On the other hand, their study for the 
full "pomeron" behaviour is outside the scope of the present thesis, and we refer to the 
available literature on the subject pc 
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Figure 3.2: Examples of two-scale processes: (a) double DIS: two highly virtual photons couples 
to a gluon ladder by quark loops; (b) a high Et jet is emitted in the forward region with respect 
to the incoming proton. 



On the contrary, we are going to extensively study the next-to-leading Inl/x (NLx) 
corrections to the BFKL equation and to the resummed anomalous dimensions in semi- 
hard processes with one or more hard scales. 

In this chapter we present a formula which generalizes at NLx level the factorization 
properties already determined at Lx accuracy by Eq. ( p.55|) and we compute on one hand 
the NLx BFKL kernel ||2^,^ , which is supposed to be independent of the external probes, 
and on the other hand the impact factors pB|[27[], which characterize the probe. 



The separation of impact factors and kernel is made on the basis of high-energy k- 
dependent factorization (Sec. p.2.1| ) which is therefore to be extended at NLx level . 

Particular care has to be taken in the use of a proper factorization scheme, because 
the separation of the cross section in impact factor contributions and gluon Green's func- 
tion suffers of some ambiguity, analogous to the one between coefficient functions and 
anomalous dimensions in coUinear factorization. In fact, the subtraction of the leading 
In s terms involves a prescription which is not unique, not only for choosing the scale of s, 
which is undetermined in Lx approximation, but also for the form of at finite energies. 
Therefore, some probe-independent NLx terms can be attributed to either the impact 
factors or to the Green's function, depending on the scheme being adopted. 

First of all, we shall introduce a definition of impact factors on the basis of a proper 
scale choice in such a way to prevent the presence of spurious infrared divergences. The 
definition, carried out for partonic particles in a first time, will be generalized to the case 
of colourless probes, which is the most important one for applications to DIS and to heavy 
quark processes. 
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Secondly, we shall proceed to the determination of the NLx part of the BFKL kernel, 
discussing at length the relevance of the choice of the energy-scale sq. 

Finally, we shall extract the physical quantity, i.e., the resummed anomalous dimension 
and the hard pomeron intercept stemming from the NL calculation, by discussing also 
the running coupling features. By pointing out the large size of the NLa; corrections, we 
will prepare the ground for the next chapter, where a further improvement of the BFKL 
equation will cure a lot of pathologies of the Inl/x hierarchy. 



3.1 Next-to- leading high energy factorization 

In the calculation of the Lx cross section, two basic properties have led us to the factorized 
form ('E.55I): 



• the dominance, in multi Regge kinematics (MRK), of the amplitude with the gluon 
quantum numbers in the t-channel (only for them there is no cancellation between 
s and u channel contributions); 

• the gluon reggeization, which determines in each t-channel gluon exchange the cor- 
rection to the propagator (|2.41| ) coming from all order leading virtual contributions. 

The above properties holds also in the NLx approximation, as the explicit calculations 
show. 



In conclusion, we will assume the factorization formula ( p. 55 ) to be valid also at NLx 



level, with appropriate NLx impact factors and GGF. This assumption has to be checked 
"a posteriori" , and we anticipate that this is actually the case, even though the GGF has 
not a simple resolvent structure as in Eq. ( p.53|) . 



3.2 Impact factors 

Among the ambiguities concerning the definition of the NLx impact factors and kernel, 
the most important is due to the dependence of some NLx features on the determination 
of the energy-scale Sq in the In s dependence of the cross section. 

In hard processes like DIS, the scale of s is taken to be Q^, the virtuality of the photon 
or EW boson involved. Thus the SF are basically dependent on the Bjorken scaling 
variable x ~ Q'^/s, with scaling violations induced by asiQ"^). Similar considerations can 



be made for double DIS [20| of quarkonium production, where two hard scale are present. 
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On the other hand, the computation of cross sections in process where the probe 
doesn't couple directly to the gluons involves diagrams with at least one loop more than 
the ones with partonic probes. For this reason, the high energy cluster expansion (see 
Sec. 3.5. 1|) needed for the definition of the NLx kernel has been mostly investigated in the 
case of parton-parton scattering, in which no physical hard scale is present. 

While the total cross section of this process has severe power-like Coulomb singulari- 
ties, it is hoped that by fixing the transverse momenta k and of the fragmentation jets 
(corresponding to the virtualities of the exchanged gluons. Fig. one is able to define 
a two-scale hard process. A similar but not identical procedure was devised by Mueller 
and Navelet in hadron-hadron scattering (cfr. Sec. p.3|) . 

We then consider high energy scattering of two partons a, b = quark (q) or gluon (g) 
with momenta pi and p2 and we factorize it in impact factors and GGF contribution as 
in Fig. O. 




Figure 3.3: High energy factorization diagram for parton-parton scattering in dijet production. 



Formally, we can write the colour and polarization averaged differential cross section 



do- 



ab 



K{k) Guj{k, ko) hb{ko 



(3.9) 



dfcdfco J 27ri \so{k,ko) 

By assuming the GGF to be the resolvent of a NLx kernel G^j = [^^ ^ /C]~^ and by 
expanding in ctg the impact factors and the kernel 



h{k) = h^'^^k) + a,h^^\k) + ■ ■ ■ 
K,{k, k') = a.K^"^^ (fc, k') + alK^^^ {k,k') + --- 

we obtain, after a;-integration, the 1-loop expansion of 



dcr 



ab 



dfc dfco 



So 



'■a "-b 



+ 



(3.10) 

(3.11) 



(3.12) 
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It is easy to recognize, in the second term of the above expansion, the leading In s term 
whose coefficient provides the Lx BFKL kernel, whilst the constant (in s) ones yields 
the impact factor corrections. However, the ambiguity in the energy scale Sq affects the 
definition of the impact factor correction, and a "scale fixing" prescription has to be 
introduced. 



We define the impact factors by the following factorization procedure p6i[27[] : 



(i) subtract the leading term with a reference scale, consistently with the IR and 
coUinear properties of the process; 

(ii) interpret the remaining constant at single-fc factorization level as the 1-loop correc- 
tion h^^\ 

The definition of the impact factor correction h^^^ can be accomplished by the 1- 
loop calculation of the differential cross section, which involves the known particle- 
particle-reggeon (PPR) vertex at virtual level (the reggeon trajectory correction giving 
a Lx contribution) and require an accurate treatment of the squared matrix element for 
one particle production in the fragmentation regionsQof the incoming particles p6|,|27|. 



Because of the UV singularities stemming from the loop integrals in the virtual correc- 
tion diagrams, we adopt dimensional regularization in D = A + 26 space-time dimensions. 
In this way we regularize at the same time the IR singularities when considering real and 
virtual emission separately. Definitions and details are deferred to App. O. 



3.2.1 Quark impact factor 

Let's start by considering quark-quark scattering in which an additional gluon g of mo- 
mentum q is produced. According to the Sudakov parametrization 

ki = zipi + Z1P2 + fci (fci = -k) Zl = ^— , (3.13a) 

(1 - Zi)S 

k2 = Z2P1 + Z2P2 + k2 (fc2 = fco) Z2 = -- ^— , (3.13b) 

(1 - Z2)S 

we define the "rapidity" of the gluonQ y- = ln{zl^/s/\q\) which can assume the values 
y e [-Y,Y] where F:= \n{^/\q\). 

^The fragmentation region of the incoming particle a is the phase space region corresponding to 
outgoing particles being almost coUinear to a. In the CM frame, this corresponds to large and positive 
rapidities with respect to a. 

^The usual rapidity variable would he y = In ((zi + Z2)y/s/\q\) which reduces to the the definition of 
the text for Z2 ^ zi, i.e., in the central region and in the fragmentation region of a. 
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In the central region, corresponding to \y\ ^ Y, the squared matrix element is correctly 
obtained by the amplitude 

A ^^^^ 1 ^ 



,2£ 



dzi dfcidfc2 



(3.14) 



which is nothing but Eq. (p.44|) without reggeization of the gluon. 

Outside the central region {\y\ ^ 1), Eq. ( p.l4| ) does not hold at NLx accuracy. The 
real emission amplitude, assuming the gluon in the a quark fragmentation region, can be 
computed by means of high energy factorization 



\M 



qq^qgq I 



u 4 cd 1 2 



(3.15) 



c,d,X 



i.e., by coupling the fragmentation tensor A^f^^ of Fig. |3.4| to the b quark with an eikonal 
vertex, because of the strong ordering between the sub-energies si := {p[ + qY <^ ^2 : = 







Figure 3.4: Quark-to-gluon fragmentation tensor A'^j^. 



By using the three-body phase space at NLa; accuracy 



we get 



, (NLx) 
UOqq^qgq 

dzi dki dk2 



h(^\k,) h(^\k. 



X 



C^cts (1 - 2;i)qr-(qr - 2;ifc2) , CpOL^ 



zfki 



{q - ziki f 



vr q^{q- Zik^ 



where 



7,,{z,e) = -[l + {l-zf+ez^] 



(3.16) 



(3.17) 



is related to the quark-to-gluon splitting function. 

We can easily check that dcqq^^qsq matches dcjqgq ( ^.14| ) in the central region 2:1 <^ 1 
— the Cp term being suppressed by a factor z1 — so that we can take it as the right 
differential cross section in the whole positive q'^ rapidity range 



2/>0 



(3.18) 
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In the remaining half phase-space y < the cross section has the same expression provided 
we exchange ki ^ k2 , z\ z^. 



It is worth noting that the splitting function ( |3.17| ) is factored out in Eq. ( p.l6| ) even 



outside the collinear regions -C fe^ ~ ^2 and fc^ <^ ~ thus suggesting a smooth 



extrapolation between collinear and Regge regions pO 



The cross section in Eq. ( |3.16|) contains two colour factors which have a simple interpre- 



tation, depending on the collinear singularities involved. The Cp term, with singularities 
at = ((g — zik'i)^ = 0), comes from the Sudakov jet region, in which the emitted 
gluon is collinear to the incoming (outgoing) quark. 

On the other hand, the Ca term is not really singular at either = or (qr — 2:1^2)^ = 
0, except for zi ~ 0, which correspond to the central region. It comes from the "coherent" 
region in which the gluon is emitted at angles which are large with respect to the piPi 
angle, and is thus sensitive to the total qq' charge Ca- This is the region we are interested 
in, which is relevant for the energy scale, because it tells us how the leading matrix 
element, valid in the central region, is cut-off in the fragmentation region. 

We can simply realize the angular ordering phenomenon in the Ca part of the real 
emission contribution to the cross section ( p.l6|) upon azimuthal averaging in ki at fixed 
q, obtaining 



dCTqq^qgq 



dzi dki dk2 1 ^ TT^g^ TT 



The latter is the coherence effect prescription of Ref. [|30| in which the polar angles of the 
emitted gluon and recoiling quark obey 

^>y^ ^ Q,>Q,^,. (3.20) 

Z\ i- — Zl 



Roughly speaking, for > k^ the half phase space \/q^Js < zi < 1 is available in full, 
whilst, for kl <C q^ it is dynamically restricted to \/q^Js < zi < |qr|/(|fci| + |qr|) ~ |qr|/|fci|. 
With the parallel consideration in the other half phase space, the longitudinal integration 
of ( p.l9| ) provides the logarithmic factor 

£^cb, %,Memi - - + {1 „ 2} . '" Max(,.fc.)Max(,.fc,) f''") 

suggesting the natural scale of s to be Sq = sm'-= Max(qr, ki) Max(qr, ^2) and hence the 
subtraction to be adopted in order to isolate the 1-loop impact factors. 



63 



More precisely, we define as leading real emission term of the differential cross section 
in the half phase space y > the quantity 



Q{\q\ - Zi\ki\) o^Ca 1 




(3.22) 



Zi TT TTeQ V S 

which we subtract from Eq. (|3.16|) and interpret as the real emission contribution to 



the ash^a\ki)h^^' {k2) term in the expansion 



.(0) 



In this way we have fulfilled the 



first prescription for the definition of the impact factors, because the {Ca part of the) 
subtracted term is finite in the q ^ limit even for z = and therefore, according to the 
second prescription, after integration over all but one transverse variable, it defines an 
impact factor correction with the right coUinear singularities. Explicitly, after subtracting 
( p.22|) from (|3.16|) , dropping the "lower" impact factor hq\k2), and integrating over 
zi > a/ q'^/ s and fci at fixed ^2 we get 



real 



Ca 



dfci /iW(fci) 





3 



dzi Tgq(2i, e 
2^'(1) - 



[I- Zi)q-{q - zik2) Q{\q\ - Zi\ki\) 



^2 



zikiY 



(3.23a) 



for the Ca part of real emission and 



% [ dfcl h^^\k,)^ !\z, Tgq(^l,£) 
t^A J T^eQ Jo 



zfk^ 



[q - zik2y 



c, 



as 



C, 



+ 4-^'(l) 



^2 



(3.23b) 



1 3 

for the Cp part. Note that the lower bound of the integrations has been pushed down 
to zero since, in the s oo limit, the error is negligible. The real emission correction to 
the lower impact factor are provided by the left half phase space contribution y < 0. 

The 1/e pole in Eq. (|3.23b|) comes as no surprise, because it corresponds to the coUinear 
divergence at fc^ = due to q ^ g transition of the initial massless quark. This is what 
we intended with "consistent with the coUinear properties" . A different subtraction could 
have produced a double pole in e which is not consistent with the coUinear singularity. 

As far as the Cp term is concerned, we do have a double pole in e. This is not a 
problem, because, as we shall see soon, the whole Cp contribution is cancelled by virtual 
corrections. 

The correction to the cross section due to virtual emission, including subleading effects, 
for general parton-parton scattering, can be extracted from the amplitude of Ref. [B^] 



l + l](°)(-t) In 



-t 



1 + asn 



{+) 

bb 



J- bb' 



(3.24) 
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where F^^, is the particle-gluon vertex of Eqs. ( |2.3CI| ). The Ins term provides the 1-loop 
Regge gluon trajectory (Eq. ( p.32b|) ) and the finite terms n*^'^^(A;^) define the PPR vertex 
correction in the hehcity conserving channel. 



In the quark case 



n(+) = - 

2 



11 



I2e 



85 7r2\ Nf f I 5 \ Cf f I 3 ^ 



36 



26 



2\ £ 



k 

(3.25) 



The virtual contribution to the 1-loop cross section is given by the interference between 
(|3.24| ) and the Born amplitude ( |2.29|) together with the two-body phase space measure 
( |B.3D yielding 



do", 



(1) 

qq^qq 



dfci dfc2 



2r]W(fc?)ln^ + 4nW(fc?) 



S{ki-k2) 



(3.26) 



The first term in the above equation is the virtual part of the leading BFKL kernel. The 
second one provides the virtual contribution to the last two terms in Eq. ( p.9|) , so that 



h^^\k, 
h^^\k. 



virt 

Ca 
virt 



11 



85 vr^ 



Nf 



12e ^ V 36 ^ 4 / ^ A/", V 6£ 1^ 

.2\ £ 



k'^V 



(3.27a) 
(3.27b) 



When adding the RHS of Eqs. (|3.23|) and ( |3.27|) , one can immediately see the cancellation 
of the Cp terms and obtain the complete 1-loop impact factor correction [E^] 



/.«(fc) 



1 fllN^ - 2Nf 



12Nr 



1 ^3 _ 1 \ 1 /67 _ TT^ _ SA^ 
2e V 2 2 2 V 18 6 9Nr 



2\ £ 



(3.28) 



3.2.2 Gluon impact factor 

We consider now gluon-gluon scattering. In this case we have to distinguish two different 
final states for three-particle production. In addition to the gluon g' there may be: 



i) two gluons; 



a) a quark-antiquark pair. 
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gg ggg 



Assume that the additional gluon is emitted in the gluon a fragmentation region. By 
using again high energy factorization, the scattering amphtude in the Feynman gauge is 
given by 



gg^ggg 1 3 5 ^Mi/i3M5i^ 1,2 

rvry 



1 b 

g 2^2 Tb 6x^X2 

2 



(3.29) 



where the amphtude A, corresponding to the diagrams of Fig. was found in Ref. ^1 
to be 



+ 



'P3 ^ P5 
+ I /is ^ At5 

, as ^ 05 , 



-Pi ^ P5 
/il <-> /i5 

oi ^ as. 



where the current 



D^(x,i/,z) = — 
x-y 



y-z - P2-P4 ) P2 + ^^^{y-z - P2-Pa)x^ + {x-p2)y^ 

P2-Z J X-Z 



(3.30) 



(3.31) 



is a function of the momenta x", y", z", with x-y = x"ya- 



The squared hehcity amphtudes corresponding to Eq. ( p.30|) were given in Ref. 



Since we work in D = 4 + 2e dimensions, we perform exphcitly the polarization sum on 
Eq. ( |3.29| ). The averaged squared matrix element that we obtain is actually independent 
of the space-time dimensionality (i.e. ^-independent) so that we finally get 



dcr, 



H Hf7ur = ^i°Hfci)4°^(fc2)^^ 

azi dki ak2 & " 



X 



CaOs zfk^ + (1 - zi)^q'^ + 2:1(1 - zi)ki-q 



7T 



{ziki + (1 - Zi)q) 



(3.32) 



This expression explicitly exhibits the symmetry in the exchange — fci ^ q ; zi ^ 1 — zi 
due to the identity of the gluons emitted in the fragmentation region of gluon a, but will 
be used for the softer gluon only {zi < i), in order to avoid double counting. 
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The function 



4 + 



2^i(l-^i) 



(3.33) 



is related to the gluon to gluon sphtting function, and is factored out in Eq. ( p.32| ) in the 
whole fragmentation region, as for quark scattering. In this case also the cross section 



matches the leading expression ( 3.14 ) in the central region 2:1 <^ 1, and thus is taken 
to be valid in the half phase space ( p.l8| ) in which two of the three gluons have positive 
rapidity. 



gg ^ qqg 



The main contribution to the cross section from this kind of final states is reached when 



the fermion pair belongs to the same fragmentation vertex, as in Fig. |3.6| a. Other graphs 
like Fig. P]^b are suppressed by a factor \t\/s because of fermionic exchange. 





{b) 



Figure 3.6: Quark-antiquark emission amplitudes: (a) The pair is emitted in the fragmentation 
region of the upper incoming gluon; (b) This kind of diagram is suppressed by a factor \t\/s. 



Assuming qq being emitted in the fragmentation region of gluon a, and labelling q 
with "3" and q with "5" (ps = g), the corresponding amplitude is 



M 



qqg 



K2 



where Bifl is the sum of the diagrams depicted in Fig. |37 

This calculation is very similar to the qq qgq cross section, and yields 
ddn 



^qqg Ar ^(0)/i, ^^ Qis-^e '^ggj^l-,^) 



X 



2vr, 

Caois -^i(1 - zi)ki-q ^ CfOs J_ 
TT g2 (^ziki + (1 — zi)qy n 



(3.34) 



(3.35) 



67 



H , a 




V , b 






Figure 3.7: Gluon to quark fragmentation tensor i?^^. 



in which the function 



'^cig{,zi,e) — 1 



2zi{l-zi] 
1+e 



(3.36) 



related to the gluon to quark splitting function, is regular in Zi G [0, 1]. Hence these final 
states do not produce any In s term, and need no subtraction. 

The virtual correction to the amplitude, including NLx effects, is given by an expres- 
sion analogue to Eq. ( |3.24| ) where Ili^\k'^) is replaced by the gluon- gluon- reggeon vertex 
correction in the helicity conserving channel 



— K + 



11 

12e 



67 vr2 \ Nf 

h 4^ (1) 

36 4 ^ ^ 7 N, 



l---^'(l) 
6e 18 ' 



(3.37) 



Now, with the result of experience of Sec. |3.2.1| , we can obtain the gluon impact factor in 
a more direct way: let's introduce the effective fragmentation vertex Fa(zi,fei,qr) of the 
incoming particle a defined in such a way that the 1-loop differential high energy cross 
section for ab scattering with possible emission of a new particle c with positive rapidity 
with respect to a is given by Fa times the exchanged Regge-gluon propagator and the full 
impact factor of particle b. In practice the effective fragmentation vertex represents the 
cross section of a with a Regge-gluon, including PPR virtual corrections and one-particle 



emission in NLx approximation, as in Fig. 3.8 





Figure 3.8: Effective parton Regge-gluon fragmentation vertex F^. In the RHS the circular 



blobs denotes the fragmentation tensors of Figs. 3.4, 3.t and 3.7; the elliptic blobs corresponds 
to the PPR vertices n(+) . 



In order to extract the 1-loop correction to the impact factor, we have to eliminate 
the Ins term by subtracting the leading fragmentation vertex — i.e., Eq. ( p. 22 ) without 
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the second impact factor — containing the singularity coming from the sphtting 

e(g - ziki] 



functions ^Pga: 



a. 



1 



'0 T^eQ" 

From Eqs. ( |3.32| ), ( p.35| ) and ( |3.37| ), we get the gluon fragmentation vertex 



(3.38) 



(3.39) 



+ Nf — 



"s Tqg(2l,£:) 



-zi{l - zi)ki-q 



Cf 1 



Lg2(^ifci + (l-^i)g)' Ca 

zlk\ + (1 - z^fq^ + 2;i(l - 2;i)fci-g 



g2 (z^fci + {l-z^))q^ 



(3.40) 



Having identified zi as the momentum fraction of the softer gluon in gg ggg scat- 
tering, to avoid double counting due to final gluons identity, we restrict the naive half 
phase space to zi G [■^, |]. In so doing, the 1-loop correction to the gluon impact factor 
is given by 

I (llNc-2Nf\ 1 1/ llN^+{2 + e)Nf 



h^\k) 



e V 12iVe 



2Nre 



6 



(3.41) 



1 CF 2NfTR /2 1 \ 1 /67 5iV/ 
^2iV^a4 e U ^ 3 V ^ 2 lis ~ T ~ 9iV^ 



2\ £ 



/.W(fe) . 



3.2.3 Colourless impact factors 



It is possible to define impact factors for colourless sources — as in double-DIS processes 
— in a way analogous to the one adopted for partons. In this case, there is an additional 
dependence on the u variable and on the external hard scale(s) so that, at leading level, 
the impact factors 



(3.42) 



yield a nontrivial fci-dependence, which has been explicitly computed in several proces- 
ses [Q. In Eq. ( p.42| ) the function /a has the role of setting ki of order Qi in the total 
cross section, and has no fe^ = singularity at all, because there is no initial colour charge. 

Therefore, by translating angular ordering into the cu-space threshold factor (q/ki)'^, 
the analogue of Eq. ( p.38|) becomes 



c^sh^P^iQi, fcs) = dki dzi Fa,^(Qi, zi,ki,q) - a^hfl^Qi, ki] 



LOTleq^ \ki 



a.e(|fci|-|q|) 



(3.43) 
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In this case, no g = nor fci = singularities are expected in F^, so that it is again 
important that the leading term be subtracted out with the angular ordering constraint. 



We thus conclude that the definition of the impact factors in Eqs. ( |3.38| ) and ( |3.43| ) 



defines a self-consistent fc-factorization scheme of Eq. (|3.9| ) for both coloured and colourless 
sources. 



3.3 Collinear factorization and finite parts 

We want now to analyze the structure of the quark and gluon impact factors derived in 
Sees. ( |3.2.1| ) and ( |3.2.2| ). The explicit expressions reported in Eqs. ( |3.28| ) and ( |3.41| ) show 



the presence of 



a /9-function coefficient — h 



71 , llNr-2Ni 



the finite part Pba{uJ = 0) in the u expansion of the Mellin transform of the a ^ b = 
q,g splitting functions for the incoming parton a under consideration, as depicted 
in Tab. fO); 



, , . , /C /67 5Nf 
a common and e- finite factor — : = 



Nr V 18 6 9Nr 



c^a 




J^z^P,,{z,s) dz 


q q 


Cf 




O + Oioj) 


g q 


Cf 


z J 




q g 




^ 2z(l-z)l 
1 + e J 




g g 


2Nc 




2N, nN, + i2 + e)Nj 


i . 6 '^^-^^ 


6 ^^^""^ 



Table 3.1: Partonic splitting functions and Mellin transforms in uj-space. 
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This suggests writing the complete one- loop impact factors ha{k) = h^a\k) + ash^a\k) 
for both a = q, g in the following manner: 



a. 



_^^/,(o)(/,) + J- y h^^){k)^ + -^h^^\k) ic (^Y 



(3.44) 



bo fk 

1 a. 



2\ £■ 



2tx 



/.f)(fc)/c+5:/.(°)(fc)^(^ 



2\ e 



The factor in front of the last expression provides the renormalization of the coupling 
constant, and by introducing the one- loop running coupling 



as{k^) = as 



1 - &o — — 



(3.45) 



it can be incorporated in hf^ itself. 

Therefore, the impact factors in Eq. ( |3.44| ) assume the form 



/ia(fc') = /ii°)(fc;«s(fc'))(i + ^/c) + 
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/,(o):^ + /,(o):^ 

^ e ^ e 



2\ e 



(3.46) 



and thus satisfy the DGLAP equations with splitting functions Pcaioj = 0). 

After factorization of the coUinear singularities contained in the latter, the only finite 
renormalization is the one implied by the factor (^^ + ^ ^) > which is universal, i.e. 
independent of the parton type. 

Several comments are in order. First, the collinear singularities in do not contain 
the ~ l/cij terms of the splitting functions, which have been subtracted out in the leading 
term. In this context, we remark the difference between the double- A; cross section defined 
here through fc-factorization, and the double-mini-jet inclusive cross section defined by 
Mueller and Navelet 



In the latter case the cross section is inclusive over all fragmentation products of the 
incoming partons which are not identified, and measures the gluonic fe's only because of 
collinear strong ordering. Therefore, it factorizes the gluon distribution density in the 
parton, rather than the impact factor, with all its collinear singularities included and, 
furthermore, can be applied only at leading In s level because of strong ordering. 

On the other hand, in our case the definition of the fc-dependent cross section is more 
precise, but can be done theoretically rather than experimentally, because one needs not 
only to identify all fragmentation products, but also to subtract out the central region 
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tail — which appears to be hard to do experimentally. Because of this difference, the 
coUinear singularities to be factored out are different in the two cases. 

As a second point we remark that a universal renormalization, with the same /C coef- 



ficient, holds also for the soft part of the one-loop time-like splitting functions [0. This 
part has a next-to- leading 1/(1 — zi) singularity, the leading one having a logarithmic 
factor, and is therefore analogous to the impact factor, which corresponds to the NLa; 
constant piece in the high-energy limit. 



The above analogy is perhaps a hint |35] to explain the universality found here. How- 
ever our result does not seem to follow in a clearcut way, because of the difference between 
coUinear and high-energy factorization pointed out before. 



3.4 Which of the energy-scales? 



The problem of choosing the energy scale so{k, ko) is not concerned only with the coUinear 
behaviour of the impact factors, but also with the validity of the NLx factorization formula 
( p.9|) to all orders in and with the definition of the GGF. 

One could ask whether Eq. ( p.9|) really holds to all orders when sq = sm is adopted 
as energy-scale. If not, is it possible to properly define the impact factors and the GGF 
in such a way that, for a different choice of Sq the factorization formula is valid? 

These questions naturally arise when noticing that the scale sm of Eq. ( p.21| ) has the 
defect of not being factorized in its ki,k2 dependence. This fact is argued to contra- 
dict multi-Regge factorization of production amplitudes which implies short-range 
correlations of the fc's. Furthermore, the use of 0(2, 1) variables [^] implies that Regge 
behaviour in the energy s should involve the boost (in the following fcj := |A;j|, q:= \q\) 

s 



cosh(C) 



2kik2 



(3.47) 



thus suggesting the choice Sq = kik2. 

We are going to show that, if the factorized scale sq = kik2 is chosen, then at 1- 
loop level the definition of the GGF has to be modified with respect to Eq. (|2.53|) by 
introducing two cu-independent additional kernel. 

The 1-loop cross section restricted to half the phase space zi < q/y/s can be repre- 
sented, apart from the b impact factor and the Regge-gluon trajectory, by the effective 
fragmentation vertex 



dzi F^{zi,ki,q) 



In 



In 



Max(fci, q) 



+ V^{ki,q) 



ki 



- Qgk^ In-^ + r]a{ki,q) 
fci 



(3.48) 
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(©qfci being a short hand notation for Q{q — ki)) where the constant (in s) function r^a 
vanishes for q = and is related to the 1-loop impact factor by 

/ dfci /^(o)(fci)^^a(fci,g) = ash^^Hk^) . (3.49) 

fixed fc2 

By adding the second half of the phase space, which is simply obtained by interchanging 
a ^ b and 1 2, we reproduce the leading logarithmic term \ns/kik2 with the desired 
energy scale, the ?7's provides the impact factor corrections and two new terms appear. 



contributing to the 1-loop cross section as [EO] 



/i(°)(fci) a,HL{k,,k2) + asHn{ki,k2) 



H:=HR{ki,k,) = HL{k2,k^ 



a. 



h':\k,) 



0gfc2 In 



(3.50) 
(3.51) 



-nsq" ' ' k2 

According to Eq. ( p.l2|) , the impact kernels and cannot enter the leading kernel 
/^(o) ^ furnishing a subleading contribution and, at the same time, they cannot be absorbed 
in the impact factors, because of their unacceptable collinear behaviour 



j dk, hf\k,)H{k,,k2) = h^^\k, 



2£2r(i + e; 



+ 0{e) 



(3.52) 



as the double e-pole points out. 

Therefore, the kernels Hl and have to be included in the GGF as subleading 
contributions in the impact points: the equation [EB] 



to be interpreted at operatorial level, replaces Eq. (p.53|) . 



(3.53) 



To summarize, if we adopt a collinear safe energy-scale, meaning that such a choice 
provides automatically IR finite impact factors, the impact kernels can be set equal to 
zero, but in this case the energy term in Eq. (|3.9|) is not factorized in its fc, fco dependence; 
on the other hand, if we want to use a factorized energy-scale, we have to introduce the 
impact kernels ( |3.51|) , in order to retain IR finite impact factors. Of course, the final 
result is the same. 



3.5 Two- loop analysis 

The calculation of NLx correction to the BFKL equation has required a large theoretical 
effort, the main part of which has been done by the Russian school and the Florence group. 
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Besides the difficulties of the technical calculation for evaluating 2-loops amplitudes, the 
determination of the NLx part of the BFKL kernel involves several conceptual problems 
mainly concerned with the identification and the subtraction of the leading term in a 
consistent way. 

The previous section had the purpose of cleaning up the subtraction procedure, but 
it is to be said that, chronologically, such an analysis was the last step but one in order 
to derive the NLx kernel. For its determination, in the "simplest" case of parton-parton 
scattering, a 2-loops calculation is needed. In fact, by expanding Eq. ( |3.9| ) to 2-loops in 
NLx approximation by using Eqs. ( p.53| , p.l0| , |3.11|) , we obtain, in addition to Eq. ( p.l2| ) 



plus the contribution of the impact kernels (|3.50|) , 



-c^^jln^A 
+ ln 



^^(0)^(0)2^(0) 



dfcdfeo I Sq 



(3.54) 



up to NNLx terms. It is evident that K^^'^ can be unambiguously extracted once the 
energy-scale Sq, the 1-loop impact factor h^^^ and the impact kernel H^, Hji are known. 



3.5.1 The cluster expansion 

To obtain the 2-loops differential cross section in NLx approximation, we have to consider, 
besides real emission in MRK and Lx virtual corrections, also processes in which the 
additional power of ctg is not compensated by a In s. This is accomplished by admitting 
that any (but only one) pair of the outgoing particles carry a small — i.e., not increasing 
with energy — invariant mass. This is known as quasi-multi-Regge kinematics (QMRK), 
in distinction to MRK where all the sub-energies of all pairs are large compared to the 
transverse momenta. At virtual level, we have to consider corrections not including one 
In s. 

After these premises, we are led to consider three kinds of processes: (i) elastic diffu- 
sion, (ii) one particle and (Hi) two particle production. In (ii) and (Hi) there can be a 
cluster of two particles with small invariant mass. 



Elastic diffusion 



The simplest process whose NLx part of the kernel is identified is the elastic one (i). As 
Fig. shows, there are three different 2-loops structures contributing at NLx level: the 
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(V 



V. 



1-loop 



'Lx 



NLx 



Lx 



2-loop 



2^2 



2—4 



Figure 3.9: Cluster expansion. The various class of cut diagrams contributing to the 1-loop 
and 2-loops Lx + NLx cross section are represented. Symmetric virtual correction diagrams are 
understood. The dashed/dotted lines joins classes of diagrams that are combined in order to 
form effective fragmentation vertices in the differential cross section. 



leading gluon trajectory iteration, the product of gluon trajectory and PPR vertex cor- 
rection and the pure NLx gluon trajectory. The former, given in Eq. ( p.32b|) , belongs to 
the iteration of the Lx kernel in Eq. ( |3.54| ). Also the PPR corrections, defined as the non 
logarithmic contribution to the 1-loop elastic amplitude has already been derived p9 



Therefore, the latter contribution, the NLx gluon trajectory, can be unambiguously de- 
fined, by observing that in the elastic process the only physical scale is provided by the 
transferred momentum —t^k^. In dimensional regularization OS^^ reads 
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Three-particle final states 



When considering one-particle production, the full amplitude with 1-loop corrections can 
be written 



M 



(1) 



2^3 



M 



(0) 



2^3 



1 + 



+ 1.(0) ^,2^ ^ + m 4) + 1.(0) (,2) (m p + m 4) ] } + m(H3 (3.56) 



where Mgl^^g is given in Eq. (g), si and S2 are the sub-energies variables defined in 
Eq. ( p.41| ), n is the vertex corrections of Eqs. ( p.24| , p.25| , p.37| ) and M2-I3 is the reggeon- 
reggeon-gluon (RRG) vertex correction amplitude. 

The RRG vertex corrections didn't take place in the 1-loop evaluation of the cross 
section and their interference with the Born amplitude M2^^ give a pure NLx contribution 
to 2-loops. On the contrary, the trajectory and PPR corrections are important in the 1- 
loop expansion ( p.l2| ) contributing to both the Lx and NLx term, as explained in Sec. p.4| . 
Therefore, M2H3 provides a so called "irreducible" contribution to the NLx kernel, whereas 
the remaining terms of M2H3 contribute to the NLx kernel after the subtractions of the 
leading term, of the impact factor corrections and of the impact kernels. 

The three ensuing structure contributing to the 2-loops cross section (n("''\ .(o) and 



M) are drawn in Fig. ^]9|. The one-gluon irreducible contribution is given by p9| 



da. 



{2,irr) 
2^3 



dfci dfc2 



a^/^Wirff) In 
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(3.57) 



2n {q^f^ 
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In 



kt 



2\ 1 



where, for simplicity, an azimuthal average has been performed. 

We now consider the one-particle emission with trajectory corrections diagrams to- 
gether with the elastic ones carrying one PPR and one trajectory correction, as the dashed 



line in Fig. indicates. 

Let's restrict to particle emission in the positive rapidity half space y > and let's 



take the elastic diagrams with PPR correction only in the upper line, as in Fig. |3.10| , the 
lower impact factor having been omitted. The first two diagrams are nothing but the 
components of the effective fragmentation vertex (cfr. Fig. p.8|) times the leading virtual 
term il^^^ In s in Eq. ( |3.56D and yield 



da. 



(2) 
2-+3 



dfci dfc2 



y>0 
a+b 



/i(°Hfc2) r d^iF3(^i,fci,g)l.(o)(fc 



hf{k^)^^'\k 



In ^ + In 

zs 



dzi Fa{zi,ki,q)\n- 



(3.58) 
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Figure 3.10: Reducible contribution to the NLx kernel from three-particle final states: PPR 
corrections in the upper line (the symmetric virtual correction is understood) and real emission 
in the positive rapidity half space y > with 1-loop reggeon correction in the lower and upper 
propagator. 



since S2 — Zis. It is easy to realize that we get Ins terms from both the integration 
boundaries if you replace the NLx fragmentation vertex with the leading one 



1 a« 



Z TTsq 



(3.59) 



and perform the 2;-integration. Therefore, the effective vertex Fa provides an important 
NLx correction with respect to its Lx analogue. 

On the other hand, the third diagram in Fig. |3.10| yields 



da. 



(2) 
2^3 



dfci dfc2 



y>0 



'q/Vs 



dz,F,{z,,k,,q)ln^^\kl) 



1 *i , 1 ^1 



dzi F^{zi,ki,q)\n- 



(3.60) 



since si ~ / zi. Here is another effective fragmentation vertex which, just like Fa, 
reduces to the RHS of Eq. ( |3.59| ) in the small 2;-region and for all values oi z > qj \fs in 
\jX approximation. 



However, Eq. ( 3.60 ) provides a logarithmic contribution only in the vicinity of the 
lower integration boundary, and by making the substitution (|3.59|) for Fa, the error is 
only NNLx. 

In this way, by taking the derivative of Eqs. ( ^.58| , |3.60D with respect to Ins and by 



using Eq. ( |3.48|) , we can perform the 2;i-integration exactly obtaining 



d do^"^ 



2-»3 



51ns dfci dk2 



a+b+c 



(3.61) 
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In-^ + ^2 In ^ + 2(^2^/1 
kik2 ki 



a. 



+ 2n2HL \ /if^ 



where we have set fij := fi*^°''(fc^) and ?7j := ?7a, (fc«, q)- Including the other half phase space 
y < Q ({1 2}, {L ^ R\)i inverting the Ins derivative and carrying out the transverse 
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integration of the //-terms in order to reproduce the impact factor corrections, we end up 



with the "reducible" part of the one-particle emission cross section 

a. 



do". 



(2,red) 



2^3 



dfci dfc2 



\ (fii + ^2 



In^ 
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^]2)^ In ^ + Hi:m2 + i^iHr 
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(3.62) 



Four-particle final states 

There remains to consider two-particle production at Born level (and the left over class 
of diagrams of the previous section). The contribution coming from a qq pair in the final 
state does not require any particular identification procedure, since it doesn't appear at 
Lx level, giving therefore an irreducible contribution 



do". 
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2^4 
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(3.63) 



(3.64) 



The two-gluon production, on the contrary, gives Lx and NLx contributions, and a 



treatment similar to that of Sec. |3.5.1| has to be done. First of all, we identify another 
irreducible contribution to the kernel, which consists in the production of a gluon pair in 
the central region with a small invariant mass. 




Figure 3.11: Two-gluon emission in the central region in QMRK. 



To be explicit, according to Fig. 3.11, we parametrize the exchanged momenta by 



1,2) 



(3.65a) 
(3.65b) 
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so that the interesting kinematic region corresponds to the QMRK 



l>;Zl~-2>-22^ — . 

s 



(3.66) 



The irreducible contribution to the 2 — > 4 cross section is then defined by subtracting, in 
the region (|3.66| ), the leading term: 
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(2,irr) 
2^4 
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The irreducible contribution to the cross section is [p£, 42 



rln-(2'i") 



dfci dfcj 



s"'a -'^2g 



h 



(0) 



(3.68) 



in terms of the two-gluon irreducible kernel [|39|,|42 
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where we have introduced the "collinear" kernel 
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in which an azimuthal average has been performed, and the dilogarithmic one 
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(3.71) 



where := Max(fei, fc^^ 

At this point, to complete the calculation, it remains to evaluate the subtracted term 
in the RHS of Eq. ( p.67|) to NLa; accuracy. To do this, we add to it the one-particle 
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emission diagrams with PPR corrections (see the dotted hne in Fig. |3.9| ) and perform the 
following analysis: let qi be the gluon with the largest rapidity and define 

y,:=lnz,^ , G [-F,, F,] , r,:=ln^. (3.72) 
Qi Qi 

Working at fixed transverse momenta, we split the phase space into three regions: 

A: < ?/2 < Z/i < , 

B : -Fa < 1/2 < < < Fi , 

C : -Fa < 2/2 < yi < . 



Then we add the one-particle emission diagrams with PPR corrections that we distinguish 



y <0 



C 



y>0 



y<0 



B 

Figure 3.12: How 2^3 and 2 — > 4 diagrams are combined in order to reconstruct the effective 
fragmentation vertices. Symmetric virtual corrections diagrams are understood. 



into four types — according to the sign of the rapidity of the emitted particle and to the 
(upper or lower) leg which the PPR correction is attached to — and we associate them 
to the three regions A, B and C as in Fig. p.l2| . Repeating the analysis of Sec. p. 5.1 



we reproduce the effective fragmentation vertex in the relevant phase space regions, and 
remembering that, when yi ~ ?/2 the subtracted term is exactly given by the last one in 
Eq. (|3.67|) , we are allowed to write, at NLa; accuracy. 



J (2,red) _ _ 

^^^-^ - h^^\k,)h^^\k,) / dfc ^^(p^ + PB + Pc) , (3.73) 
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where 
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and the "reduced" fragmentation vertex 3^^ is 



(3.74) 



By deriving with respect to Ins, performing the integrals and inserting the resuhs in 
Eq. we obtain p5 
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(3.75) 



the small circle o denoting composition of operators, i.e., a ^-integration. 



3.5.2 Next-to-leading BFKL kernel 



The left over elastic contributions (the first and the third diagram in the 2 — >■ 2, 2-loops 



group of Fig. p^ ) contribute to the cross section with 
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(2) 
2^2 
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At least we are at the final step: by adding Eqs. (P . 5 7| , P . 6 2| , |3 . 6 3|j3l68|J3 . 75| J3l7^ ) and 



identifying the corresponding terms in the expansion ( 3.54 ), we are able to single out the 
NLx BFKL kernel at energy-scale Sq = kik2 p 



= 2Q^'\kl)6ik, - fc2) + + K^^ + K^f + 
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where in the second line we have explicitly shown the scale-invariance violating term 
and K^^^ defines the scale-invariant part of the NLx kernel. Note that the coefficient of 
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the In/i^ term is precisely the leading kernel with the /3-function coefficient Bq given in 
Eq. (|1.2| ). Therefore it can be interpreted as a running factor and we can express the total 
Lx + NLx kernel in the form 



)C = as{jj, 



(l-6oas(/x')ln4) 



(3.78a) 

,{kl)[K('^+a.{^,')K^l]+0{a!) (3.78b) 

Factorizing the running coupling at the scale k^ is an asymmetric procedure, but is con- 
venient for the discussion of the non-scale-invariant BFKL equation. Using a different 
scale (e.g., as{k'i_)) implies changing K^^l so as to leave the total NL kernel invariant (see 
below) . 



3.5.3 Eigenvalues of the NLx kernel 

In order to investigate the physical features emerging from the NLx kernel, we determine 
the eigenvalues of the scale-invariant part by applying K, to the (spherically symmetric) 
eigenfunctions of the leading kernel: 

a,{k') I dk'{K('\k,k')+a,K!^l{k,k')){k"y-' = a,{k') [x^'\j) + a,x^'\j)] {k'^-' . 

(3.79) 

The leading eigenvalue x*'^^ is the same of Eq. ( p.60|) with m = 0: 

X(°H7) = 2^(1) - ^(7) - ^(1 - 7) , (3.80) 

while the eigenvalue function x^^^ of the next-to-leading kernel K^^^ (= K''^^ from now 
on) is the sum of a gluonic part 
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and a quark — A^j- dependent — one 
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We want to stress that the particular form of the NLx kernel, and hence the NLx 
eigenvalue, depends on the choice of both the scale of as and the scale of the energy sq. 
For instance, if we adopt as the argument of ds, in Eq. ( ^.771 ) we should replace 

Mn4i^(°) = -Mn4^(°) -6efc,.,lnMirW , (3.84) 

so that the NLx kernel undergoes a shift 



Air(i) = -60,, In = -6777^^ In ^ (3.85) 

k( |fci — fc2r fcf 

and, correspondingly, the NLx eigenvalue function acquires the additional term 

Ax^'H7) = -&^'(l-7) (3.86) 

which symmetrizes completely itself. 

In DIS the energy dependence of the SF is commonly expressed by means of the 
Bjorken variable x ~ Q"^/ s, s being the photon-proton CM energy squared. In the limit of 

^ fcp, where fep is the typical transverse momentum scale of the proton's constituents, 
the correct energy-scale is Sq = fc^ ~ Q^. The shift in the NLx kernel corresponding to 
the change Sq = kk^ k"^ can be accomplished by a simple trick. If we denote with -ff (7) 
the eigenvalue function of the impact kernel H of Eq. ( |3.51| ) 



2^^ r(7)n! ; (7 + n)2 27^^"^^^ 



(3.87) 



the GGF (actually, its scale-invariant part) has the representation (cfr. Eq. ( p.61| ) with 
m = 0) 

dcu d7 / s \- , , /PV+^" 1 ^ ^ ^ 

where 



l + as[ij^(7) + g(l-7)] . , , ^ ^ 
^.. 7 = r (ou ^ ^ - (Df M ^0 = • 3.88c 

Therefore the energy-scale change kko k'^ is equivalent to replacing 7 with 7 + ^u. By 
shifting 7 in g^^ and then expanding in u up to NLx level, the GGF Mellin transform at 
the "upper" energy-scale sq = k'^ reads 

Mr ^ - l + as[g(7)-k^°)X7)+^(l-7)] . _ oQ^ 

^^^^^ ^-«s[x(°H7) + «.x(^H7)-^«.x(°)x(°)'] ^ ^ ^ 
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This means that the shift in the kernel's eigenvalue function is 
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The shift ( |3.90|) cancels the cubic singularity (and the ({3) term) of 
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(3.91) 



providing the right collinear behaviour for k"^ ^ /cq, as explained in the next section. 
Furthermore, by Eq. ( ^.871 ), the numerator in Eq. ( |3.89| ) becomes regular for 7 = too, 
where it takes the value 1 + asH{l) = 1 — ag'ip'i^), which renormalizes the impact factors. 
This confirms the importance of factorizing H in the GGF in Eq. (|3.53|) . 



3.6 NLx resummed anomalous dimension 

In order to determine the resummed anomalous dimension in NLx approximation, we 
have to check whether the NLx BFKL equation is consistent with the RG equation. The 
running of the coupling spoils the arguments of Sec. |2.3.2| and a new analysis is needed. 

The choice of k"^ as the scale of as is suitable for that analysis. In fact, the inhomoge- 
neous BFKL equation for the unintegrated gluon density (|2.62| ) reads 

uj[j^^{k)-J^i^\k)]=-^Jdk'K{k,k')J^^{k') , b:=^b, (3.92) 

where we have explicitly written the perturbative expression of the running coupling as{k'^) 
and the scale- invariant Lx + NLx kernel K = K^^^+a^K^^K In 7-space, Eq. ( p.92|) assumes 
a simpler form: with the definition ( |U.19c ) and by observing that 61nA;^/A^ —>■ —bd^ and 



K ^ X = X*'"'* + Q^sX''^'*) we obtain 

-6c.[J-^(7)-^(°H7)r = x(7)-^.(7) , r = (3.93) 
which is a first order linear differential equation for J-'^^ admitting the formal solution 0] 

^^(7) = r dy e^ J-(°)'(7') , (3.94a) 

J 00 



d7 f k 



.2 \ 7 
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where X is any primitive of x- 



It has been argued [^,^ that, in the hmit A;^ ^ /cq, the GGF assume the factorized 
form 



6„(fe,/!„) = F„(fc)F„(fc„) X 



(3.95) 



where F is purely perturbative while F is determined by the properties of non perturbative 
physics. After integration with the bare fco-dependent gluon distribution (Eq. (|2.28| )), the 
unintegrated gluon density J-'uj{k) inherits the factorized form 



J^uj{k) = F^{k) f^ + higher twist . 



(3.96) 



A comparison with Eqs. ( 3.94| ) led us to identify the non perturbative factor 



+ 



oo— le ^ oo+i£ 



dy e^4°)'(y) , 



(3.97) 



where we have set the upper integration limit 7 ~ (see below). The perturbative factor, 
on the other hand, carries the whole fc-dependence of and is given by 



Fik)=^i^ r^Ve-^ 

' nk^ J 2m UV 



(3.98) 



In the anomalous dimension limit k'^ ^ A^, Eq. ( p.98| ) can be evaluated by the saddle 
point method, which yields |^ 



(3.99) 



where 7 = 7(^(t) = dt[^t — X{^)/huj\ is the position of the saddle point determined by the 
condition 



(3.100) 



The representation ( ^.99 ) is valid only in the anomalous dimension regime 



hut > Xm:= Min y(7) • 

76]0,1[ 



(3.101) 



If the above limitation were not satisfied, the representation ( |3.99| ) would breaks down, 
due to large 7-fluctuations. 
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Finally, we consider the singlet^ density ( p.27| ) 



2m VA^ J 7 



exp J %{t) dr {tQ := In . (3.102) 



7 V 27r|x'(7)| 

The coefficient function in front of the exponential term is simply a constant in the ^ 
(7 — > 0) limit and the singlet density satisfies 

;fL^\Q')=UQ')fL''\Q'), (3.103) 



dlnQ2' 

consistently with the RG equation (|1.36|) . Therefore, we interpret the function implic- 
itly defined by ( p.lOO| ) as the (highest eigenvalue of the) resummed anomalous dimension 
in the singlet sector. 

The structure of the NLx eigenvalue function at coupling-scale and energy-scale k"^ is 
(cfr. Eqs. ( F^ICTQ ) 

X^'\l) = ^ + — + A, + 0{^) (3.104) 
This imply that, up to 3-loops, the NLx anomalous dimension has the expansion 

7..(as) = 7l(^) + «s7nl(^) (^7nl = (3-105) 



= - + a, f + A2- + As^] + Oiat) (3.106) 

which can be checked to be consistent with the known expression up to 0(o?^). It is 
apparent that, in order to recover the expansion (|3.106|) , the double pole of x'^^^ and 
hence the shift to scale sq = fc^, plays a crucial role: if a cubic singularity were been 
present, the anomalous dimension would assume the non realistic form 7 ~ c^J"^ juj^l'^ . 



3.7 Pomeron: perturbative versus non perturbative 
features 

The study of the cu-singularities of the GGF in the NLa; approximation shows a completely 

different scenario with respect to the Lx one. A first point concerns the running of the 

^In the NLcc approximation the solution of the BFKL equation, taking contributions also from quark 
intermediate states, has to be identified with the eigenvector of the singlet sector. In fact the anomalous 



dimension in Eq. (3.106) corresponds to the larger eigenvalue of the singlet anomalous dimension matrix. 
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coupling, which changes drastically the spectral properties of the GGF. In particular, the 
regularization of a^i^Q"^) to values of around and below the Landau pole = 
affects the position and the strength of the singularities of . 



The second point, to which we would address the reader's attention, concerns the 
nature of those singularities. Note that both the perturbative and the non perturbative 
factors in Eq. ( |3.95| ) carry an tu-dependence and could contribute to the singular behaviour 
of the GGF. 

The questions then arise: 

• are both and singular in u;? 

• if so, which of them carries the rightmost cu-singularity? 

The structure of the cu-singular points and their dependence on the particular regulariza- 
tion of the IR coupling has been studied in some models whose evolution equation 
resembles in their main features the BFKL running coupling one. It turns out that the 
rightmost singularity of the GGF (the pomeron Uf) is a really non perturbative phe- 
nomenon, which is very much dependent on the behaviour of the strong coupling in the 
soft region fc2 ~ A2. For that reason, it should be related to the universal singurarity of 
soft physics, also. 

0.5 0.3 " ' 




10^ 
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10^ 



10" 
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Figure 3.13: The BFKL hard pomeron a;p(as) in leading and next-to-leading approximation. 



As far as the perturbative factor is concerned, the integral representation ( |3.98| ) has an 
essential singularity at u = 0. However, the anomalous dimension representation ( p. 99 ), 
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depending on uj through 7^^, is singular in correspondence of the branch cut u;p(as) of 
u ^ (as) given by 

^p(as) ■=asXrn ^ as[x^°^ (^) + "sX^^^ (^) 

=4^"^^ (as) [1 - 6.47as] • (3. 107) 

Note that the value of uOf^{as) is not affected by the choice of the energy-scale within the 
factorized class Sq = k^k^'^, since the shift Ax^^Hl/2) is proportional to x^°^'(l/2) = 



(cfr. Eq. (^m )- 



Due to the formal analogy of the definitions ( |3.107| ) and ( p.73|) , ujf>{as) is often called 



the hard pomeron singularity. To be fair, however, only the hard nature of ujp{as) is 
a fact. In fact we cannot guarantee it to be the rightmost singularity of the GGF. In 
addition, there may be sound doubts it to be a true singularity. In fact, uj]i>{as) signals 
the breakdown of the anomalous dimension representation (|3.99|) where the saddle point 
estimates does no more work. 

Nevertheless, even if Ljp{as) doesn't correspond to a singular behaviour of Q^^, but 
rather to a singularity of the anomalous dimension, it may be related to a power-like 
behaviour in an intermediate small-x moderate-Q^ region. This assertion is confirmed by 
studies of some simplified models [^,^ one of which we shall deal with in Chap. |^. 



The last term in the RHS of Eq. ( |3.107|) emphasizes the large and negative corrections 
stemming from the NLx BFKL kernel. In Fig. p.l3| you can see that the maximum value of 
ujp{as) — 0.11 is reached for as — 0.08 and in the range of HERA data ~ 1^10^ GeV^, 
i.e., as — 0.1 -7- 0.4, the NLx hard pomeron is very small and even negative, in deep 
disagreement with the data (see Figs. |0| ). 

The situation is even worse if we look for the solution of Eq. (|3.107|) for finite (non 
vanishing) values of as- 

Mas) = Min a,[x^°H7) + "sX^'Ht)] • (3.108) 

7e]o,i[ 



The plot of the eigenvalue function in Fig. |3.14| shows that a stable minimum for x exists 



only for very small values of as < 0.04, whereas for as > 0.04 the shape of x is the 
opposite of the Lx eigenvalue. Among the physical consequences of this pathological 
behaviour there is the loss of positivity for the total cross section, as it was pointed out 
in Refs. |3,|S|. 



Does a reliable phenomenological analysis require the high energy In 1/z resummation 
in NNLx approximation? A simple arithmetic-geometric extrapolation based on the Lx 
and NLx computation-time would suggest something like 20 -7- 100 years of work. But it 
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Figure 3.14: The BFKL eigenvalue function for different values of as and Nj = 0. It is apparent 
the change in shape and the disappearance of the minimum for Os > 0.04. 

seems that the whole In 1/x BFKL hierarchy has slow or even bad convergence problems. 
Rather than a 3-loops calculation it would be better to add some physical information, 
some consistency constraints with well established physics... 



Chapter 4 



Improvement of the small-x equation 
by RG analysis 



After many years of calculations and expectations, the NLx corrections to the kernel of 
the BFKL equation turned out to be so large as to question the very meaning of the high 
energy expansion ( p.9| ), rising the compelling question of how to improve it. 

Exact higher order calculations would be for sure a formidable task, both for their 



size and because they mix with unitarity effects and thus cannot be described within 
the BFKL equation alone. Furthermore, there seems to be serious problems with the 
convergence of the kernel, as we will show in the next section. 

The main issue is to understand the reasons of the pathological behaviour of the 
NLx corrections. Just like for the problems concerning the subtraction of the leading 
terms from the reducible contributions of the NLx kernel, the discussion of the large NLx 
correction is fundamentally related to the issue of the choice of scales. 



4.1 Origin of the double logarithms 

In a high energy scattering process of two objects with transverse scales and fcg, the 
Lx BFKL equation resums all the leading logarithmic terms of the cross section of the 
form 

asK^°Un—) , n^O (4.1) 
So J 

and we have already remarked that, in Lx approximation, the choice of sq is immaterial. 

In the NLx approximation, the energy-scale is intimately connected to the form of 
the NLx kernel (Sec. |3.5.2|) and also to the structure of the GGF (Sec. p^) . However, 
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there is not a particular scale which seems to be preferable, rather there are various 
candidates according to the regime being investigated. We have already discussed the 
collinear-safe scale sm and the factorized Regge-motivated ones kko, k"^. Nevertheless, in 
the situation where k"^ ^ k^, the DGLAP approach (Chap. |I]) tells us that each power of 
(is is accompanied by a large Iu/c^/Zcq and that the correct scaling variable, carrying the 
whole s-dependence, is given by a; ~ k'^/s. In practice, the cross section has terms like 

1 / P s \ 

as In 72 In — {k^ > k^) . (4.2) 



1.2 \ ^ h.2 h.2 



Now, if we rewrite the general term ( [4.2|) by adopting kko as the scale of the energy, 
double logarithms (In^) of the transverse momenta appear in the perturbative series: 



For instance, the O^a"^) cross section at Lx level reads (see Eq. ( |3.54| )) 



- a'^h^^\k)h^^\ko) \n' -K^'^'ik, fco) . (4.4) 



dfcdfco s a V . b V 

In the collinear limit fc ^ feg we have Sq = k'^ and 

1 1 



xW2(fe,fco) ^ J dk' 



1 



TT{k - k')^ n{k' - koY 



k^ 



(4.5) 



which shows the single logarithmic term and, in Mellin 7-space, gives rise to a double 
pole as/7^. 

By shifting the energy-scale sq = fc^ — > kko, we get 

^ - ^ F '"^ F - F '4 AS - F F + "'^^^ ■ 

rv rv rvQ rv rv rvQ rvrvQ rv rvQ rvrvQ 

The Lx term, i.e., the coefficient of In^ s, remains unchanged whilst the NLx term acquires 
the double logarithmic contribution K^'^'^ ln^(fc^//cQ) whose eigenvalue gives rise to a cubic 
pole as/7^. 

In general, the appropriate DGLAP single logs corresponds to and «"/(! — 7)'^ 

poles with k ^ n + 1. On the contrary, the presence of double logs involves poles with 
n + l<k^2n+l. 



We conclude with Salam [£0| that the Os-expansion of the kernel in conjunction with 
a non collinear-safe energy scale, generate double logs, i.e., strong collinear singularities, 
in the lower order kernels. 
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4.2 A toy kernel 

Now, with the aid of a toy model we want to show that those strong colhnear 
singularities are responsible of the instability of the finite order as In 1 /a; truncation. We 



can get an idea of that mechanism by considering the Lund model consisting in a 
small-x evolution equation for the unintegrated gluon density — derived in the context 
of the linked- dip ole- chain (LCD) model — similar to the BFKL one but differing in the 
colhnear region. The relevant equation can be written 



= + /C^ , (4.7) 



where u is the variable conjugate to x = k^/s — so that we are considering sq = fc^ — 
and the kernel JC^^ admits the usual power-like eigenfunctions (|2.59|) and its eigenvalue is 



asX(7, ^) = as[2^(l) - ^(7) - V^(l - 7 + ^)] ■ (4-8) 

In complete analogy with the method of solution explained in the context of the BFKL 
equation, the anomalous dimension 'juj^ds) is determined by the equation 

'^ = asXil,^) (4.9) 

or, equivalently, 

= ttsXeff(7,as) , (4.10) 

where XefT(7, «s) is implicitly defined by the solution of Eq. ([4.9|) once 7 and ds have been 
fixed. 

In the DGLAP limit, corresponding to the region close to 7 = 0, the effective eigenvalue 
function has only simple poles I/7 to all orders in ds'- 

Xefr(7, «s) = - [1 + 0(«s)] + 0(7°) (so = k') . (4.11) 

7 

In the opposite DGLAP limit, k^ ^ k"^, the relevant region of 7 is close to 7 = 1. Making 
the transformation to the relevant Bjorken scale Sq = k^, using 7 — 7 + C(j (cfr. Sees. p. 5. 3 
and ^73|) , there is only the pole 1/(1 — 7): 

Xesil, «s) = [1 + ^("s)] + 0{{l- 7)°) (^0 = kl) . (4.12) 

So one is completely free of double transverse logarithms in both DGLAP limits. 
In terms of the symmetric energy-scale Sq = kk^, Eq. ( |4.8| ) becomes (7—^7 + ^1^) 

X(7, uj) = 2^(1) - ^(7 + i^) - ^(1 - 7 + iw) . (4.13) 
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The as-expansion of the ensuing effective eigenvalue function 

oo 

Xefr(7,«s) = 5^«rxi?(7), (4.14a) 



n=0 

xS(7) = 2^(1) - ^(7) - ^(1 - 7) , (4.14b) 
xS(7) = -^xlJ(7)[^'(7) + V^'(l -7)] , (4.14c) 
xS(7) = -^xS(7)[^'(7) + V^'(l - 7)] - ^XS(7)^K(7) + ^"(1 - 7)] , (4.14d) 
signals the presence of double logarithms from NLx level on, as one can see from 

^S(7) = + ^(^) + {7 1 - 7} , (4.15a) 

xS(7) = +^ + <^{^s) + {7 - 1 - 7} . (4.15b) 

Note that the most divergent part of xilr in the regions of 7 close to and close to 
1 is exactly the same of the NLx BFKL eigenvalue (at scale sq = kko)\ Therefore, as 
previously remarked, the cubic poles of Eq. ( p. 791 ) are a consequence or the fact that the 
high-energy-factorization scale kko is not collinear safe, i.e., is different from the Bjorken 
scale k'i. 



If we look at the pomeron singularity ujp{as) = ttsXcfr(l/2) of the toy-kernel in Fig. 
we observe an alternate behaviour of the perturbative estimates, whose convergence is 
rather slow and limited to very low values of in the lower approximations. For realistic 
values of as, say 0.2, the fixed-order expansion are quite unreliable, and going from NLx 
to NNLx does not bring one any closer to the all-order result except for < 0.1. 

To summarize, what one learns from this toy kernel is that the double transverse 
logarithms lead to a very poor convergence of the BFKL kernel at subleading orders 
and in particular the features of the analytic structure of the kernel and of the physi- 
cal observables after the resummation are quite unrelated to the NLx ones. Therefore, 
essential ingredients of any improvement of the BFKL approach should be the correct 
treatment of the collinear behaviour, as predicted by the RG, and the resummation of the 
corresponding collinear singularities to all orders. 



4.3 Improved small- cc kernel 

In small-x physics, Regge theory and the RG come to a sort of clash, which provides non 
trivial consistency requirements. The general form of the RG improved small-x kernel 
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is constrained hj (i) the exact Lx and NLx calculation and (ii) the coUinear singulatity 
structure which we are going to specify. The latter requires to study the kernel in the 
coUinear limits k ^ kg and k <^ k^ where the relevant energy-scales are k"^ and k^ 
respectively. It is also useful to consider the symmetric high-energy-factorization scale 
kkf). 

It is easy to see that, under energy-scale transformations of that kind, the GGF in 



Eq. (|3.9|) transforms according to 



(4.16) 



where the symmetric "s", upper "u" and lower "/" suffixes refer to the energy-scales kk^, 
k^ and k^. The similarity transformations ( 4.16| ) apply straightforwardly to the small-x 
kernel defined by 



g. = [uj - ic^]-' 

which we should consider tu-dependent [ p^ : 



^] L{k,ko) 



(4.17) 



(4.18) 
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Figure 4.1: The toy-kernel pomeron singlularity as a function of Os, at leading, NLx, NNLx 
N^Lx and all orders. 
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Note that in the definition ( [4.17D of our improved small-x kernel, the GGF does not 
contain any impact kerneL 

We have aheady expressed the opinion that cxg is not a suitable expansion parameter 
for the kernel. The introduction of the w-dependence of /C^ will reveal crucial in order to 
derive a novel small-x expansion which automatically resums the coUinear singularities. 

From a perturbative point of view, however, we observe that, in general, the renor- 
malized kernel /Ct^(fc, k'; /i^, Q;s(yU^)) for non vanishing values of its arguments k and k' , is 
IR finite (see Sees. p.3.1| and |3.5.3|) and RG invariant, so that it can be expanded as a 



power series in as{k ) with scale-invariant coefficient kernels 



/C^(fc,fc';/i^a,(/x2)) = J2 [<ys{k')Y^'Ki-\k,k') . (4.19) 

n=0 



In other words, the only non scale-invariant source is the running coupling, which in 
Eq. ( |4.19| ) has been evaluated, for simplicity, at scale k^. Different coupling-scales corre- 



spond to different, but always scale- invariant, coefficient kernels. 



4.4 Form of the collinear singularities 



We have shown in Sec. p.6| (see also Refs. |^,^) that the BFKL equation satisfies RG 
factorization in an asymptotic way. The asymptotic form of the GGF for t := In fc^/A^ ^ 
to := Iu/cq/A^ is given by 

k''gUk,ko) ~ a(as(t)) exp / 7^(«s(r)) dr FM , (4-20) 

where is the larger eigenvalue of the singlet anomalous dimension matrix, defined by 
the saddle point condition (|3.100|) . 

Therefore, the RG invariant kernel in the LHS of Eq. ( [4.19| ) acquires collinear singu- 
rarities for k'/k —>■ {k/k' 0), which corresponds to strong ordering of the transverse 
momenta in the direction of the "upper" Bjorken scale k"^ ("lower" scale k^). Such sin- 
gularities are due to the non singular part 7 of the singlet anomalous dimension (p.l06|) 
which, neglecting the (small) qq contribution, is 

70. = 7f - - = as^i(^) + asM^) + " • • , (4-21) 

UJ 
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the singular part being taken into account by the BFKL iteration itself. It follows |£2 
that, for k ^ k', 

L{k,k') ^ ^ exp I'^Mt)) dr (t » t') (4.22) 

= ^(l-^«s(t)ln^ 

Expanding Eq. ( |4.22| ) in as{t) and comparing with the general definition ( [4.19| ), leads 
to the identification of the kernels K}^ in the collinear limit, whose eigenvalue functions 
turn out to have the singularities 

+ + (4.23) 

which correspond to single logarithmic scaling violations for k ^ k^. A similar reasoning 

1 

yields the collinear behavior of /C^^ in the opposite strong ordering region k' ^ k 

L(k, k') ~ exp !' iM.{t)) dT («' » i) (4.24) 



and to the singularities 

X^'W^ ^-'^'^^;),!^--'"" (1-7«1). (4.25) 

u 1 u 

However, the similarity relations ( 4.18| ) connect the kernels /C and /C. Therefore /C has 

1 

the singularities (f4.25|) shifted at 7 = l+u also, and similarly /C has the singularities ( |4.23| ) 

s 

shifted at 7 = —uj. As a consequence, the symmetric kernel /C — for the energy-scale 
So = kk^ — has both kinds of singularities shifted by ±u;/2: 

Xl Kl) ^ (^^i^)n+i ^ -2^) ' (4-26a) 

1- (Ai-b)(Ai-2b)---(Ai-nb) , , , . ,^ 

Note the 6-dependent asymmetry of the singularities in Eq. ( f4.26| ) under the 7^1 — 7 
transformation. It is due to the fact that the expansion ( ^4.19|) involves a^ii) (and not 
as(t')). Of course, the kernel K.^ itself must be symmetric under t ^ i! exchange, so that 
expressing aj(t!) in terms of a.J(C) 



«s(0 = ^^^^^ (4.27) 

l-6as(t) In — 
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leads to the symmetry constraints (in the following we omit the suffix "s") 



"^(^) = E Q i-bd.r-'xtHl - 7) • (4.28) 



It is straightforward to check by the binomial identity 

I \ " / \ / 
r + n \ I r \ I n 



( TTj ) ( TTj ) (4.29) 



n / — ' \ m / \ m 

' m=0 ^ ' ^ 

that the symmetry constraints (^4.28|) are indeed satisfied by Eq. (|4.26|) . In particular we 
must have 

xL^Hl - 7) = XL°H7) , xl;ni-7) = xL^n7) + &xL°^'(7), (4.30) 

showing that the antisymmetric part of xi^''(7) is ici). 

The coefficient kernels -ft'd"^ take up coUinear singularities not only from the non singu- 
lar part of the gluon anomalous dimension 7^^, but also from qq states which are coupled 
to it in the one- loop gluon/quark-sea anomalous dimension matrix 

7.:=7.--fc °) =«sA(u;). (4.31) 

Although the inclusion of the two-channel evolution (|4.31|) involves the collinear problem 
conceptually, the numerical effect of the quark-sea contribution to the gluon anomalous 
dimension is pretty small p4| , |53| . Therefore, in the following, we will restrict our analysis 
to the gluon sector only. 



4.4.1 Form of the leading coefficient kernel 

Having determined the RG constraints, we are going to develop a procedure for defining 
the improved small-x kernel at NLx accuracy. The improved coefficient kernels k]^^ are 
constructed by requiring that 

(1) the Green's function reproduce the known NLx calculations; 

(2) the collinear singularities be as in Eq. ( [4.26|) . 

In order to implement condition (Ijwe have first to relate the w-dependent formulation 
of in Eq.( [4.17l ) to the customary expression of the BFKL kernel at NLx level 



(4.32) 
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The cu-dependent formulation of Eq. ( [4.19| ) yields instead the NLx expansion 



(4.33) 



which is actually more general than Eq. ( [4.32|) because the a^oJ term, coming from the 
cu-expansion , is a possible NLx contribution too. 

Now it turns out that, at NLx level, the formulation ( |4.33| ) reduces to the one in ( |4.32| ), 
provided the impact kernels and of Eq. (|3.53|) are taken into account |Q. In fact, 
by using the expansion ( 4.33 ) and simple operator identities, we can write 

[c^ - ^J-i = (1 - _ ( -^^(0) ^ -2^(1) + -1 (1 _ ^^K^o,i)yh (4^34) 

provided we set 



5 



(4.35a) 
(4.35b) 



Eqs. ( [4.34]) and (|3.53|) show that the two formulations above differ by just a redefinition of 
the impact kernels, while Eq. (|4.35b| ) means that k'^^ is given by K^^\ after subtraction 
of the term already accounted for in the w-dependence of K''^'^ . 

The simplest way to define a Lx improved coefficient kernel fulfilling Eqs. (^4.35a|) and 
( [4.26|) |yt=o is to displace the poles of xi?"* in Eq. (p.80|) by shifting the arguments of the 
^/'-f unctions, as suggested in Ref. [pO|] : 



Xi°n7) = 2^(1) - ^(7 + \u:) - ^(1 - 7 + \uj) 

,2 



(4.36) 



1 vr' 
-oj- 



+ 



sm 7r7 



The kernel K^"* , corresponding to Eq. ( |4.36| ) is that occurring in the Lund model [^] and 
is given by 



(4.37) 



where /c> := Max(A;, A;') and := Min(A;, A;'). It is thus related to the customary lead- 
ing kernel K^^^ by the "threshold factor" (k^/k^)^. This means that the s-dependence 
provided by its inverse Mellin transform is [0] 



K^'>\s;k,k') 



du / s 
27d \kk' 



^ 1 



Kj:^\k,k')=K^'^\k,k')Q{s-k 



(4.38) 
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Can one justify the form of the kernel (|4.37| ) "a priori"? From the point of view of the 
RG improved equation, any kernel which (i) reduces to K^^'^ in the cu ^ limit and (ii) 
has the leading simple poles of Eq. ( [4.26D for n = 0, is an acceptable starting point. An 
alternative choice of this kind will differ from k!;^^ by a NLx kernel without 7 = or 7 = 1 
singularities. The resulting ambiguity can thus be reabsorbed by a proper subtraction in 
the NLx coefficient kernel. 

Nevertheless, the threshold interpretation of Eqs. (|4.37| ) and (|4.38|) is appealing. For 
instance, the first iteration of such a kernel provides the expression 

2 



-A'<»> 



(4.39) 



: j dk i^(°)(fci> k) (^In ^ - r]{h, k) - r]{k2, k)^ K^'\k, k. 



where 



cosh?7(A;j, k) = 



k^ + kl 



The threshold condition implied by Eq. (|4.39|) 

——— = cosh?] > cosh (rjiki, k) 
2kik2 



ri{k, k2)) 



(4.40) 



(4.41) 



is reminiscent of phase space in Toller variables and may be regarded as an alternative 
way of stating coherence effects [3C], as implied in the original version of the Lund model 
itself. 



Whether or not such hints will eventually provide a more direct justification of 



(0) 



the fact remains that Eq. ( [4. 36] ) resums the cu-dependence of the 7-singularities, and thus 
provides the correct singularities of the scale-dependent terms of the NL kernel. Therefore, 
it is a good starting point, yielding NL contributions which are smoother than those in 
the as(t)-expansion, as we now discuss. 



4.4.2 Form of the next-to-leading contribution 

The NLx improved coefficient kernel is constrained by Eq. ( |4.26| )|ra=i and by Eq. ( |4.35b[ ) 
which yields the u = limit of the eigenvalue function 

X?io(7) = X^'Hl) + y^Hl)-^ ■ (4.42) 
2 sm 7r7 



The subtraction term so obtained is important because it has cubic poles at 7 = 0, 1 
which cancel the corresponding ones occurring in x^^Hl)- Furthermore, the impact kernels 
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of Eq. ( |4.34| ) have quadratic poles which similarly account for the ones occurring in Hj^ 
and Hr- This means that the remaining contributions are, in both cases, much smoother 
in the cj-dependent formulation. 

In order to implement condition (2) on x2'\ we note that the uj = limit ( [4.42| ) 
still contains double and single poles at 7 = 0, 1, which should be shifted according to 
Eq. ([4.26|) . By neglecting the (small) qq contributions, the explicit form of Eq. (|4.42| ), 
following from Eq. ( 3.81|) for the energy-scale Sq = kki 



0, IS 



xi=o(7) = - 




-ix^'^"{i) + ^x^'\i)- 



TT 



sm 7r7 



+ 



TT 



cos7r7 
sin 7r7 J 3(1 — 27) 

2^ 3 



11 + 



7(1-7) 



12 



x^'\i) + ^m + 



;i + 27)(3 - 27)^ 
$(7) , 



+ 



TT 



-w., (4.43) 

4 sm TT'f 

where $(7) has been defined in Eq. ( |3.82| ). Here we have singled out some singular terms 
which have a natural physical interpretation, namely the running coupling terms (in round 
brackets), the energy-scale-dependent terms (in square brackets) and the coUinear terms 
(in curly brackets). 

The running coupling terms have a double pole at 7 = 1 only, and account for the 
asymmetric part of x^^^ (given by —^X^^^') which provides the 6-dependent double pole 
on xi^"* in Eq. ([4.26|) . The coUinear terms have symmetric double poles with residue 
Ai{lj = 0), in accordance with Eq. ( [4.26| ) also. Both types of singularities can be shifted 
by adding a NNLx term, vanishing in the u = limit, which we take to be 

A,{u;)^'{^ + I) - Ai(0)^'(7) + {A^iuj) - b)^'{l - 7 + f ) - (^i(O) - b)^'{l - 7) • 

(4.44) 

This term incorporates the cu-dependence of the one-loop anomalous dimension ( |4.21|) too. 

The energy-scale-dependent term in square brackets contains the subtraction ( [4.42| ) 
and has, therefore, simple poles at 7 = 0, 1 only, which we can shift by adding the 
contribution^ 



TC 



Y(x^(7)-Xo(7)). 



(4.45) 



By then collecting Eqs. ( [4.42|) , ([4.44|) and ( [4.45| ) we obtain the final eigenvalue function 



xr(7) = Xi(7) + ^i(^)^'(7 + 1^) + (^i(^) - &)^'(1 - 7 + 1^) + yXo (7) , (4-46) 



^Of course, such simple poles, which are dependent on the choice ( [4.36| ) of xS\ do not occur — by 
construction — in the NLx eigenvalue function x^^Ht)- They are just part of the NNLx ambiguity of 
our resummation scheme, whose size is evaluated in Sec. 4.10.3. 
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where 



I sm 7r7 



TT 



(4.47) 



is a symmetric function without 7 = or 7 = 1 singularities at all. The expression ( [4.46 ) 
satisfies in addition the symmetry constraints ( [4.30|) , having antisymmetric part — Ixi^"*'. 



Of course, there is some ambiguity involved in the choice of the subtraction terms 
(|4.44| , |4.45| ), which boils down to the possibility of adding to (|4.47| ) a term, vanishing in 
the = limit, and having only higher twist 7-singularities, around 7 = — 1, —2, ■ ■ ■ and 
7 = 2, 3, ■ ■ ■ . We will discuss that ambiguity in Sec. |4.10.3. 



4.4.3 Numerical importance of collinear effects at NLO 




-15 



0.2 



0.4 



0.6 



0.8 y 



Figure 4.2: A comparison of the coUinearly-enhanced (double and triple poles only) part of the 
NLx corrections with the full NLx corrections; Nf = 0. 



Above we have given the general form for the collinear singularities of the kernel at 
all orders. It is of interest to consider at NLx level just how much of the full corrections 
come from these collinearly enhanced terms. Accordingly we look at the part of the NLx 
corrections which contains just double and triple poles, Xcoii- 

1 



^(1) 

A coll 



Ai Ai-b 1 
(1-7)2 ~2y 



2(1-7)= 



(4.48) 
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This is compared with the full x^^^ Fig- |4.2| , where we have plotted their ratios to 
X^'^^ ■ The remarkable observation is that over a range of 7, the collinear approximation 
reproduces the true corrections to within 7%. It is obviously impossible to say whether 
this is true at higher orders as well. However the fact that the study of collinear terms 
has such predictive power at NLx is a non-trivial point in favour of our resummation 
approach. 



4.5 Factorization of non-pert urbative effects 

The next step is to obtain the improved GGF ( |4.17|) , i.e., to solve 



ujQUk, ko) = 6\k -ko) + j dk' /C^(fc, k')g^{k', ko) . (4.49) 



where an extention of the representation ( [4.19|) in the region around the Landau pole 
k"^ ~ (t = 0) is needed. However, for perturbation theory to be applicable, the non 
perturbative effects of such region should be factored out, as is predicted by the RG, and 



has been argued for at Lx and NLx level 



We assume that, by a suitable regularization of asit) around the Landau pole, /C^ can 
be defined as a hermitian operator bounded from above in an Hilbert space, with a 
continuum (or possibly discrete) spectrum Sp(/C^) C ]— 00, /ip(a;)]. Typical regularizations 
of this kind may 

(a) cut-off as{t) below some value t = t > 0: 

«s(t) = ^e(t-t); (4.50) 
bt 



(b) freeze it in the form 



«.(t) = ^e(t -t) + i0(t - t) , (4.51) 



possibly with some smoothing out around the cusp. 

In such a framework, a formal solution for the Green's function Q^^ is given by the spectral 
representation 

J-oc TT UJ-fl 
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in terms of a complete and orthonormal set of (real) eigenf unctions : /i G Sp(/C^) 

/C^F^^ = f^F^ . (4.53) 



Let's refer for definiteness to the frozen-as(t) regularization ( [4.51| ), which allows a 



simple classification of the eigenfunctions F^{k) of Eq. (^33|), according to their behavior 
for t — i> — oo. In fact, since as{t) is fixed for t < t, in the t — oo limit the eigenfunctions 
must behave just like the scale-invariant ones 

(A;2)^M-i = le''^^^^' , 7 = 1 + iz/ . (4.54) 

K A 

Because of the symmetry property of the kernel with respect to the exchange of its argu- 
ments k ^ fc', to each eigenvalue [i there corresponds two opposite values ±z^(/u) and we 
can decompose the eigenfunctions F^{k) as 

F^{k) = j^lE^^^^k) - E-'''^^'\k)) (4.55) 



_ 1- (c(z/)e''^(^)* - c* Me-''^"^') 



for suitable functions E^{\z) having a plane- wave asymptotic behavior for large and neg- 
ative t (the 10 index has been dropped). 

The precise superposition coefficient civ) of left- and right-moving waves occurring in 
Eq. ( [4.55| ) is determined by the condition that F^ik) be regular for t —>■ +oo, i.e., be 
vanishing at least as rapidly as 1/k, so as to allow an normalization. 

If we assume E''^^^ and E~'^^^^ to be boundary values of an imaginary analytic function 
F^{k) of /i, whose branch cut lies along the spectrum (see Fig. |4.3| ), we can rewrite the 
spectral representation ( [4. 521 ) as a contour integral 



6„(,.,„)./ ^^i£»M. (4.56) 

encircling the spectrum Sp(/Ctj). By distorting the yU-contour (because is well behaved, 
for 3fte(iz/) > 0) and by picking up the residue at the fi = u pole, we end up with the 
factorized expression 



g^k, ko) = F^{k)F^{ko) {e > kl) (4.57) 

where F^{k) : = F^{k) and the same for F. This procedure can be carried through 
provided k/k^ is large enough for the decrease of F^ to compensate the increase of E'^^'^\ 
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Figure 4.3: The function is defined as that analytic function of fi whose hmit on its branch 
cut — corresponding to the spectrum of IC — is given by E'^^^^ (^E^'^^'^^ ) over (under) the cut. 



In this way the real integration in Eq. (4.52 ) can be replaced by a complex integration along the 
contour C{lo) of Eq. ([4.5^ . 



The plausibility argument above is further supported by the explicit model we will 
give in Chap. ^ for arbitrary values of t and to, and hints at the general validity of 
Eq. ( [4.57] ). Therefore, for k ^ ko, the Green's function is asymptotically proportional 
— up to higher twist contributions — to the regular (for t oo) solution F^{k) of the 
homogeneous BFKL equation ( [4.53| ) with eigenvalue fj, = u: 



IC F 



(4.58) 



which becomes the basic quantity to be found. 



4.6 The small-u; expansion 

Having determined the factorization property of the GGF in Eq. (|4.57| ), we are now in 
a position to determine its perturbative (large-A;) behaviour determined by the small-x 
equation (^4.58|) whose kernel we know at NLx level. In this section the scale of the energy 
So is arbitrary, but one have to keep in mind that the resulting quantities depend on the 
scale chosen. 

As usual it is better to work in 7-space, where we represent the eigenfunctions F^} in 
a formQ hinted by Eq. { ^9^ : 

^The canonical dimension of the eigenfunctions is [F] = — 1. 
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where X^(7,/i) is to be found by solving Eq. ( [4.53[ ). 

By using the above representation for and keeping the general expression (|4.19|) 
for the kernel, the eigenvalue equation ( [4.53D becomes 



= [/C^-/i]F^'(fc) 



27ri 



n=0 



We assume that in the small-//, large-t regime 



> - > - » 1 



(4.60) 



(4.61) 



— resembling the anomalous dimension regime ( p.lOlj) — the integral ( [4.60|) is dominated 
by a stable saddle point 7^(as(^)) determined by the condition 



5,|(7-i)t-^X^(7,/i)} =0 ^ h^t = d,X^M 



5^X^(7,/i) <0. 
which need to be checked "a posteriori" . 



(4.62a) 
(4.62b) 



By expanding in 7 the exponent and the eigenvalue functions xi""* around 7 = 7 and 



by keeping the ensuing fluctuations up to the relevant order, as explained in Ref. [53 



Eq. ( |4.60 ) provides relations between the eigenvalues of the improved coefficient kernels 
and the 7-derivatives of X^(7, /i) which allows us to write the latter in terms of the former. 
More precisely, in the regime ([4.61|) , we can express 



(4.63) 



in) 

by an expansion in /i whose coefficients are the eigenvalues Xw and their 7-derivatives: 



X. (7, /i) = (7) + /i r^i'^ (7) + r/(2^ (7) + • • • 



(4.64) 



where the //-expansion coefficient r^^ could be derived, in principle, to all orders, and up 
to the third order (NNNLx) reads |53| 



(0) 




+ 6r/(^)'-2r/(^(^) 



(4.65a) 
(4.65b) 

(4.65c) 
(4.65d) 
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In the 11 = oj case, Eq. ( |4.64| ) becomes an cu-expansion [^] for 



+ NNLx , 



(4.66) 



which provides, through Eqs. (|4.63| ), ( [4.59[ ) and (|4.57| ), the perturbative expression for 
the GGF we were looking for. 



The analogy of Eqs. ([4.62 



I ^l=L^J 



( glOOD , (|39D ^ (|3:98|) and the Lx truncation of 
Eq. ( [4.66|) leads to the identification of Xw(7) with the effective eigenvalue function of the 
improved kernel /C^. 




1 y 



Figure 4.4: The NLx truncation ( 4.6(\ ) of the effective eigenvalue function Xcuij); = 0- 



Before entering the details of the determination of the improved anomalous dimension 
and of the hard pomeron, we would like to spend a few words about the major features of 
the improved eigenvalue Xuj{'j)- In Fig. [4.4| you can see the effective eigenvalue function 
for various values of the new expansion parameter uj. The stability of the shape, i.e., 
the presence of a minimum even for sizeable values of u is apparent, at variance with 
Fig. |3.14| . The shift of the poles is evident too, providing the consistency with the RG. 
The asymmetric shape is due to our (asymmetric) choice of the coupling- scale k"^. 
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4.7 Improved anomalous dimension 

In order to study the anomalous dimension limit ^ of the integrated gluon density, 
let's set the energy scale for the GGF to Sq = fc^. Due to the validity of RG factorization 
in the large-t limit of Eq. ( |4.57| ), we can state that the gluon density f!f^^\Q'^) in the 
hadron H has a universal t := InQ^/A^ dependence 



where the if-dependent coefficient arise from Eq. (|2.28|) 



(4.67) 



(4.68) 



and is in general non perturbative, while the perturbative one is independent on the probe 
and is determined by 



d7 1 
27ri 7 



dk F^{k) = I 7 exp <{ - -^X^i^j) 



(4.69) 



where we have specified the function X^^ at the "upper" energy-scale k"^. 

The asymptotic behaviour of Eq. ( 4.69| ) in the RG regime can be found from the saddle 
point ( [4.62|) which, thanks to the identity 



iJ - X^il) = / 7a;(as(r)) dr + const 



(4.70) 



(check it by differentiating with respect to t) yields the result 



bio 



%y 27r\xUl)\ 



exp 



7c<;(«s(t")) dr . 



(4.71) 



The coefficient in front, coming from the saddle point fluctuations, has been evaluated at 
NL level only. If we work at NL level, the saddle point approximation ( [4.71|) is enough, 
and provides the effective anomalous dimension 



7f (as(t)) := 37 ^^Q^i^) = % - -\n 



dt 



dt 



huj 



^(7 



1 ^ 



X'^{%) V^- 2x^,(7^) 
where in the last equality use have been made of the relation 



+ 



(4.72) 



= xL(7^)^7c.(as(t)) • 



(4.73) 
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Figure 4.5: The anomalous dimension in various approximations. 



Fig. [4.5| shows the (purely gluonic) anomalous dimension as a function of uj for as = 0.2. 

'^{uj/as) and has the familiar branch-cut 



X 



(0)- 



The Lx anomalous dimension is just 7l 
at UJ = 41n2a;s- The NLx anomalous dimension ( p.l05| ) has a divergent structure around 
the same point as 7l. The resummed result, defined in Eqs. ( [1.69| ) and ( |4.72| ), shows a 
divergence at a much lower u, defined by ujc{t), corresponding to the vanishing of the 
gluon density 



9.At)it) = 0, (4.74) 

thus extending the resummed anomalous dimension beyond the spurious singularity of the 
BFKL truncation. What is particularly remarkable is the similarity to the DGLAP result 
until very close to the divergence. The momentum sum rule is automatically conserved: 
for C(j = 1 we have 7^^^^ = 0. Note however that, since ujcit) comes from a zero of guit), the 
u = ujc(t) singularity does not transfer to g^it) itself. Furthermore, the singularity of the 
anomalous dimension representation (|4.71| ), coming from the failure of the saddle point 
expansion at the value uj = Ugit) such that 

XL(7..W) = 0, (4.75) 



is different from uJc{t). 
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4.8 Improved hard pomeron 
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Figure 4.6: lOc sjid uJs as a function of Og for the BFKL kernel with Nj = 0. 



The values of the exponents and as a function of as (and Q"^), are shown in 
Fig. |4.6| and compared with the Lx and pure NLx results. It is apparent that the im- 
proved equation provides sensible predictions even for sizeable values of ctg. The above 
difference should not be too confusing. The exponent ujs{t) signals the breakdown of 
the formal small-x expansion of the anomalous dimension of Eq. ( [4.72|) , due to infinite 
saddle-point fluctuations, while ujcit) tells us the position of the singularity of the re- 
summed anomalous dimension. Their difference arises from their different definitions, not 
from some instability of our approach. 

What are the roles of the two hard pomeron estimates Uc and a;<j? What, if any, their 
relation to the pomeron singularity cjp, the leading singularity of the GGF? We are not in 
a position to say a lot about those questions yet. Probably there is an intermediate small- 
X moderate-Q^ regime where the perturbative physics dominates, where Uc or could be 
the effective small-x growth exponent. The non perturbative soft pomeron, determined 
by the spectral condition = /ip(u;p), should however govern the very low- a; region. 



In order to gain some hints on these issues, we have devised a simplified model |^ 
which allows one to treat in a semi-analytic way the whole t-space dependence and thus 
to estimete the (in) dependence of some physical quantities on the strong coupling region 
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and also to clarify possible transition mechanisms between various x-regimes. But this is 
the subject of the last chapter. 



4.9 Estimate of the NLcc truncation error 

What is the error that we make in the NL truncation of the RG improved equation? 
Our claim is that, in the improved formulation, based on the cj-expansion ([4.66|) , this 
error is smaller than in the formal NL expansion in as{t). Let us in fact estimate the 
remaining terms in Eq. ( [4.64|) . According to Eq. ([4.26|) further subleading eigenvalue 



(2) (3) 

functions contain at least higher order coUinear poles which contribute to ritj , rju and 
so on. A first observation is that, even if xL"^ has (n + l)-th order poles, the r^i"'* : ^ 
have at most simple poles, due to the powers of x^^^ i^i the denominator, roughly due to 
the replacement as{t) ~ uj/x'^^ ■ Therefore, their contribution cannot be too big, even for 
small values of 7 = O^u). 

Furthermore, one can check that, if qq contributions are neglected, the leading collinear 
poles actually cancel out in the expansions (|4.65| ) of ri!^\ r]u \ ■ ■ ■ around both 7 = and 
7=1. The mechanism of this cancellation is explained in Ref. [Q. 

From a physical point of view, it is not possible for simple poles to survive in r^i"^ : 
n ^ 2, because, when replaced in the saddle point condition ( |4.73| ), they would provide 

: 72 ^ 2 corrections to the 1-loop anomalous dimensions which cannot possibly be there. 
In fact, the full anomalous dimension is accounted for by Eqs. ( [4.66| , ^!62l) as follows 



bu;tc^- + ^ ^^ = aJ- + A,{u;)] , (4.76) 

where we have taken the small-7 limit of the collinear safe eigenvalue function Xuj{l)- 

We therefore conclude that, in the purely gluonic case, the NL tu-expansion ( |4.62| ) 
takes into account the collinear behavior to all-orders, and that no further resummation 
is needed. 



4.10 Stability 

The effective eigenvalue Xuj{'j) reveals its stability in shape providing a stable minimum 
and thus a reliable framework for small-x physics. Nevertheless, there are some am- 
biguities, one of which is proper of the present resummation method and concerns the 
pole-shifting procedure of Sec. W^, others are typical of the perturbative treatments and 



concern the renormalization scheme prescription and the renormalization scale choice. 



110 



The original Lx+NLx formalism suffered from considerable instabilities under renor- 
malization group scale and scheme changes |5^. An important characteristic of any re- 



summed approach is that it should be relatively insensitive to such changes, and generally 
stable. 



4.10.1 Renormalization scale 

Note first that in our approach the renormalization scale only enters through the RG 
invariant A parameter of t = Infc^/A^ (Eqs. ([4.19|) and ( [4.69|) ). It is then easy to see that 



the physical results are A-independent. A redefinition of A is essentially a shift in t, say 
by an amount At. There is a corresponding modification of x'^\ x'^\ • • • by the amounts 

Xi^^ - XL^^ + &Atxi?) , xi^)-xL^) + 26Atxi^)+6^(At)^xL°^ , .... (4.77) 



In the 7" representation ( [4.59|) , this corresponds to a modification of X^^ by an amount 
bfi'-fAt. In fact the transformation ( |4.77| ) changes the coefficient rj^^ only, the remaining 
ones ri!^\ ri^\ . . . being left invariant. This change exactly cancels the modification of t 
itself: 

exp - ^X,(7,/u)| ^ exp S^-f{t + At) - i-(X^(7,^) + 6^7 At) | , (4.78) 
thus implying that the physical results are independent of the A-parameter choice . 



This automatic resummation of the renormalization scale alleviates the need for techniques 



such as BLM resummation |5g], advocated for example in |5^, which show a strong 
renormalization scheme dependence. 

4.10.2 Renormalization scale 

The issue of renormalization scheme dependence is in fact closely related. Consider a 
scheme S related to the MS scheme by 

af) = ap) + ra^ (4.79) 

with an appropriate modification of x^^^ ■ Except for terms of 0{a^) and higher, this is 
identical to a renormalization scale change. Indeed if one defines the scheme change by a 
modification of A then renormalization scheme changes behave exactly as renormalization 
scale changes, and so have no effect on the answer. Using instead |4.79| there is some 



residual dependence on the scheme at 0{al), but as one can see in Fig. for the T 



scheme, which has T ^ 1.17 (for = 0), the effect of the change of scheme is small. 
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Figure 4.7: Renormalization scheme uncertainty of the two exponents; MS scheme and T 
scheme; Os is always shown in the MS scheme, and is connected to the T scheme value of ag 



via (4.75). 



4.10.3 Resummation scheme 

In resumming the double transverse logarithms (energy-scale terms), there is some free- 
dom in one's choice of how to shift the poles around 7 = and 7 = 1. In a similar manner 
to what was done in 0] we consider two choices. The one explicitly discussed in this 
chapter can be summarized as 

with an equivalent procedure around 7 = 1. We refer to this as resummation scheme (a). 
An alternative possibility is (scheme (b)) 



1 



1 



7" (7 + icj)*^ ■ 

For instance, in this second scheme, the Lx improved eigenfunction would be 

1 



7 1 - 7 7 + 



7 + 5^ 



(4.80) 



A comparison of these two resummation schemes pBI is given in Fig. |4.(j| and the difference 
between them is again reasonably small. 
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Figure 4.8: Resummation scheme uncertainty of the two critical exponents ujs and ujc- 



Aside from the explicit renormalization-scale independence, the stabihty of our ap- 
proach is connected with the resummation of the colhnear poles, for both the double-log, 
energy-scale dependent terms (the 1/7^ and 1/(1 — 7)^ poles at NLx order) and for the 
single-log ones of Eq. (|4.76| ) . 



Chapter 5 

The collinear model 



Let's briefly review the results of the previous chapters. The study of perturbative high 
energy processes requires resummation techniques in order to correctly take into account 
the energy dependence. The BFKL approach provides such a framework, but exhibits an 
unstable behaviour as soon as the coupling exceeds very small values. It has been shown 
that the implementation of RG constrained collinear physics cures the BFKL hierarchy 
pathologies and the use of oj instead of as as expansion parameter furnishes a stable 
framework where to study small-x physics. 

However, the improved formulation of Chap. ^ still leaves some open questions, con- 
cerning in particular the relative importance of perturbative versus non perturbative ef- 
fects and the identification of the physical observables. Of course a lot of work has to be 
done yet in checking the validity of the factorization properties ([4.57|) of the GGF and on 
the novel small-x equation ( |4.6CI| ) determining the perturbative solution. 

For those purposes, we have devised a simple, but powerful tool for studying the 
problem of small-x physics that we call the collinear model — namely a model where 
all and only the coUinearly-enhanced physics is correctly included, in particular the full 
dependence on the one-loop running coupling, the sphtting functions, and the so-called 
energy-scale terms (the ones in square brackets in Eq. ( |4.43| ). The model has the prop- 
erties of correctly reproducing the one-loop renormalization group results and of being 
symmetric (given a problem with two transverse scales, its results do not depend on which 
of the scales is larger). These are properties desired from the resummation of the NLL 
corrections in the case of the full BFKL kernel. 

While it does not correctly resum the series of leading and subleading logarithms of 
s (i.e. the non-coUinear part of the problem), in that region it does have a structure 
which qualitatively is very similar to that of the BFKL equation (cfr. Fig. |5.1| ) and can 
usefully serve as a model. In contrast to the BFKL equation, it is very easily soluble, as 
a Schrodinger-like problem. 
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5.1 Definition of the model 



Our model |^6[ is defined in such a way to reproduce the DGLAP hmits for branchings of 
gluons with ordered transverse momenta, and the anti-DGLAP hmit for branchings with 
anti-ordered transverse momenta. In practice the kernel of this model extends the coUinear 
asymptotic form ( [4.22|) of the small-x kernel down tot:= lnA;^/A^ ^ t':= Infc'^/A^ while 
in the opposite region t ^ t' it is defined by Eq. ([4.24|) . It is clear that the two expressions 
refer to the "upper" k'^ and "lower" energy-scales of the GGF respectively. By taking 
into account the energy-scale transformation relations (|4.18|) the symmetric expression of 
the collinear kernel reads 



X^{t,t'):= kk'}C^{k,k') 



(5.1) 



exp - t') + A,{uj) j^' a,{T) drj e{t - t') + (t ^ t') 



Aiiuj) being the non singular part of the gluon anomalous dimension (|4.21|) at 1-loop. 

Therefore the kernel for a single small-x branching actually resums many branchings, 
of which the last (and only the last) is governed by the l/uj part of the splitting function 
and is taken into account in the iteration of the kernel giving the GGF 



Qu) = — 



(5.2) 



of the model. 

The collinear properties of the kernel (|5.1|) can also be seen in 7-space, by the expansion 



/C^(fc,A;') = X^[a,(P)r^KW(A;,A:') 



n=0 



where the scale-invariant kernels K^df^ have eigenvalue functions x)j^\'^) given by 
(„) 1 ■ A,{Ai + h)---{A, + {n- 1)6) , I- {A,- h){A, - 2b) ■ ■ ■ {A, - nh) 



in). 



Xr(7) 



(7 + k^r^' 



The leading order eigenvalue (cfr. Eq. (|4.35a| )) 



(5.3) 



(0) 



7 



7 



(5.4) 



differs from BFKL quite significantly numerically, but retains a very similar structure — 
a saddle point at 7 = 1/2 (see Fig. |5.1| a), implying a power growth of the cross section, 
and diffusion. The NLx eigenvalue (cfr. Eq. ( 4.35b|) ) 
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Figure 5.1: Comparison of the effective eigenvalue functions of the coUinear model and of the 
full small-x kernel of Chap. ^. In (a) the leading (uj = Oj ones and in (b) the Lx+NLx ones for 
uj = 0.3. 

(5.5) 

turns out to be very close, even numerically, to the full BFKL NLO kernel, reproducing 
it to better than 7% over the whole range of 7 from to 1 (cfr. Fig. |4.2| ). This confirms 
that collinearly enhanced effects dominate the NLx kernel. 



5.2 Differential equation formulations 



5.2.1 First order formulation 



The main advantage of our collinear kernel from the point of view of this article is its 
relative simplicity. Specifically it can be written in factorized form: 



Xuit.t') = U^it) V^it') Q{t - t') + U^it') VL(t) Q{t' - t) , 
where (the a;-dependence will be understood) 

U{t) = as{t)e^^[-^—^t + Ai{uj) [ as(r)dr 



V{t) =exp{^-^t- Ai{uj) [ as(r)dr 



(5.6) 

(5.7a) 
(5.7b) 
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This allows us to recast the homogeneous BFKL equation 

/t poo 
dt' V{t') J{t') + V{t) / dt' U{t') 'J{t') 
oo Jt 

as a differential equation. Dividing 5" into two parts, 

Ait) = U{t) f dt' V{t') -Jit') 

/oo 
dt' U{t') '^{t') 

u-Jit) = A{t) + Bit) 

and taking the derivative leads to a pair of coupled differential equations: 

d^ U' 



dt U 



dt V 



For the specific kernel ( p.l| ) that we consider, we have 

dA _ f l + uj 
dF ~ 

dB _ fl+uj 



2 Os 



— AiUs I B — as?" 



(5.8) 

(5.9a) 

(5.9b) 
(5.9c) 

(5.10a) 
(5.10b) 

(5.11a) 
(5.11b) 



Since we have two coupled equations, there are two independent solutions. Examining the 
equation for large and positive t, where as is small, one sees that they can be classified 
as a regular solution 



\Jr ~ exp I — t 



which is dominated by A, and an irregular solution 



J- 1 ~ exp I — - — t 



(5.12) 



(5.13) 



dominated by B. 



5.2.2 Second order formulation 



The coupled set of differential equations ( |5.10| ) can be recast in the form of a simple 
second order equation for 5". In fact, by using ( ^.9|) , we can first rewrite (|5.10|) in the form 



UVJ. 



(5.14) 
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Then, in order to eliminate the resolvents appearing in ( [5.14| ) we introduce the operator 



v,.,a,..^-^A(a.-^)Ja.^'^-^)(a.-Y;) (5.15) 



U W J y U ) \ V W J \ V 

^2 w' ^ U'V" - U"V' 
d + 



w w 

and the wronskian 

w{t) = W[U, V] = UV - U'V = a,{t) (l+uj- 2Ai«s - |^ ) • (5.16) 



By applying to ( |5.14| ) we finally obtain |^ 



uDt^ = ( — - — j UV3' = -w;(t)?, (5.17) 

which is the second order formulation that we were looking for. 

Eq. ( |5.17[ ) can be recast in normal form by the similarity transformation 

g^(t) = const ■ V^t)/(t) , (5.18) 

which leads to a Schrodinger type equation, 

{-d^ + V,s)f = (5.19a) 

^ f^'V ^ f^'\ U"V'-V"U' 
A \ w J 2 \ w J UJ w 

Note that the above derivation is valid for any form of the running coupling as{t) which 
extrapolates the perturbative form {bt)~^ into the strong coupling region around the 
Landau pole t = 0. 

In what follows we consider various regularizations of the Landau pole, in particular 
the cut-off and the frozen cases defined in Eqs. ([4.501 ) and (|4.51|) . It is obvious that 
such different regularizations will change the form of the potential (Fig. p.2| ) and thus 
the boundary conditions on 3^ (or /) coming from the strong coupling region. A similar 
Schrodinger formulation was found |^ for the (Airy) diffusion model @,^] with running 



coupling, with a potential which roughly corresponds to the bottom of the well in Fig. 

We consider in particular the solution 5'j?(t) (3'L(t)) of the homogeneous equation 
( [5.17|) which is regular for t — > +00 {t —00). If both conditions are satisfied, = "Jl 



is an eigenfunction of the BFKL equation. The pomeron singularity = cup is the 
maximum value of uj for which this occurs (ground state). 

\{ uo > C(jp, 7^ and the two solutions have rather different features. Due to the 
locality of the differential equation, 5"/? is independent of the regularization in the region 
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Figure 5.2: Qualitative t-dependence of the effective potential for the regularizations of type 
(a) and (b) of the coupling strength. 

t > t. On the other hand, will be dependent on the behaviour of as{t) for t < i. For 
instance, in the case of as{t) being frozen (4.51), has the exponential behaviour 



9^z.(t)~e'^S {t<t) 
where k = ^/Vcfi{t < t) is found from ( |5.19b| ) to be 



-{1 + u) - Aiasit) 



2(1 + -) 



Aia,{t) 



2as(t) 



(5.20) 



(5.21) 



5.2.3 Factorization of non-perturbative effects 

The basic tool for describing BFKL evolution is the Green's function Sujit, to) = kkoQ^{k, ko) 
which satisfies the inhomogeneous small-x equation 



toSUt, to) = 5{t - to) + X^0 to) 



(5.22) 



and is supposed to be well-behaved for t,to ±00. The problem of factorization is the 
question of the (in) dependence on the non-perturbative strong-coupling region. 

For t 7^ Sui satisfies the same differential equation as 3" and is thus a superposition 
of two independent solutions. The large-t behaviour implies a regularity condition and 
suggests the expression 



9Ut, to) = 3'R{t)3^Lito)eit - to) + 3'Lit)3'R{to)Qito -t) , (t^to 



(5.23) 



Actually (|5.23|) is a rigorous consequence of the second-order formulation. In fact, Sa 
satisfies the differential equation 



ht + -) = -Dt5it - to) . 



(5.24) 
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With a little thought, one can realize that must contain a delta function term in the 
form 

S^(t, to) = -6{t - to) + to) (5.25) 

UJ 

and hence 

^t + -)L = -^S{t-to) (5.26) 

showing that § is continuous function at t = to with discontinuous derivative. Eq. ( ^.26| ) 
is an inhomogeneous Schrodinger type equation with a delta source, and its solution is 
just like the RHS of Eq. (|5.23|) also for t = to- We conclude that 

g(t, to) = -S{t - to) + 3^nmL{to)Q{t - to) + Mt)3'R{to)e{to - t) (5.27) 

00 

with the normalization 

9^R,L(t) = ^^/ij,L(t) , W[fRjL]:=fRf'L-fnfL = l. (5.28) 

UJ 

The main consequence of ( |5.23| ) is that the regularization dependence is factorized away 
in 3^L, whenever t or to are large enough. This happens in particular in the coUinear limit 
t — to ^ 1 relevant for structure functions. 



5.3 Solutions: analytical features 

The collinear model just defined can be solved in principle as a Schrodinger problem by 
known analytical and numerical techniques and for both types of solutions occurring in 
the Green's function (|5.23| ) (i.e., the left-regular and the right-regular ones). 

The right-regular solution 5"^ is, for large t, perturbative, i.e., independent of the 
potential in the strong coupling region t ^ t , while the left-regular one is dependent 
on the strong coupling. The left-regular solution 3"^ can be written as a superposition of 
the right-regular solution S^ji and an other generic (linearly independent on 5"/?) solution 
3^ J (which is necessarily right-irregular, i.e., 3"/ — > oo for t +oo)|^. By normalizing 3^/ 
in such a way that 

W[fnJi] = l , Ut)=--^^^m, (5.29) 

UJ 

^The right-irregular solution ^F/ has to be specified by suitable boundary conditions for t — > +oo, thus 
it is independent of the strong coupling region. This is the reason why 3^i is useful for parametrizing the 
left-regular function J'l. 
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we have 



(5.30) 



The coefficient S{uj) is formally found from the matching conditions of both sides of 
Eq. ( p.30| ) and their first derivatives at t = t: 



3^R{i+) ■ Cl{uj) - Cniuj) 



C,:-- 



Hi' 



(5.31) 



The coefficient £l is provided by the strong coupling boundary conditions. In the two 
cases of regularization that we are considering we have 



-'^ I freezing — 

~ '1+iU 



C 

'^-f' I cutoff 



as{i) 
uj{l + uj) 



Kit) 

iu(l + LU) 



1 + 



u;(l + cj) 



(5.32) 



(5.33) 



where as{t) is the t ^ limit of as{t). 

Let us remark that, while 3"/ and S^r are expected to have an essential singularity 
at u; = only, (cfr. Sec. |5?7|) , the coefficient S{u) displays the leading part of the u- 
spectrum in ^euj > 0, and in particular the leading Pomeron singularity cjp. In fact, if 
u —>■ Uq & Sp(/Ctj), thenQ ~^ const x 3^ji and the normalization condition ( |5.29D can be 
satisfied only if S{uj) — > oo. 

The soft pomeron up is the greatest cj-singularity of S{u!) and in the frozen-Og regu- 
larization it corresponds to the endpoint of the continuous spectrum. In the cut-off case, 
we expect also a discrete component in the spectrum (since there is a potential well, see 
Fig. p. 21 ) with poles at u;- values such that Cl{^^) = 'Cr{lj), both being analytic functions in 
3?ea; > 0. Mixed cases are also possible, depending on the detailed shape of the potential 
in the strong coupling region. 



Eq. ( 5.30 ) allows us to rewrite the Green's function for t > and lo > up in terms of 
S{u!) in Eq. ( ^.31|) as follows: 



(5.34) 



If both t, to 3> 1, but t — to = 0{1), Eq. (|5.34|) is dominated by the ffist term, the second 
being suppressed exponentially in because the regular solution is much smaller than the 

^This argument apply in the case of luq in the discrete spectrum of ICuj. If luq belongs to the continuous 
spectrum, then S{uj) is discontinuous at cj = o^o, the singularity being a branch cut. 
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irregular one for to ^ 1- This term is, on the other hand, defined by boundary conditions 
for t, to — > +00 only, and is therefore independent of the strong coupling region. Its 
analytical and numerical form will be discussed in more detail in the following. 

The second term in Eq. ( ^.34| ) carries the regularization dependence and contains the 
leading cu-singularities, in particular (Ref. [0). In Y := \og{s/kko) space, the sum in 
Eq. ( |5.34| ) defines two asymptotic regimes, as we shall see. 



5.4 Strong-coupling features 

Within the collinear model we can study the spectrum of /C^, which provides the uj- 
singularities of the Green's function embodied in S{uj), and in turn determine the high- 
energy behaviour of the cross section. 

The leading cj-singularity is the Pomeron i^p, i.e., the maximum uj value for which 
a zero energy bound state for the potential ( p.l9b| ) is present. The pomeron properties 
depends on the coupling regularization. For smooth and non negative couplings the 
effective potential is a regular function bounded from below and increasing with uj so 
that there exists a value of = cjp > such that for uj > uj-p there are no more zero 
energy solutions. This demonstrates that the spectrum of the collinear kernel is bounded 
from above (Sp(/C(^) C ] — 00,0;?]) for any regularization (tup depends however on the 
regularization!). 

In the following, we consider the two cases of freezing ( [4.51|) and cutting-off ( [4.501 ) 



as{t). The first case is fairly simple, because for t <i the effective potential is constant 
and at its minimum value. Therefore wp, the upper limit of the continuous spectrum, 
given by V^^{t < t) — 0, is determined by the condition 

1 + cjp — la^\t)Ax\yjTp) 



This is to be compared with the saddle point definition of the hard Pomeron ( [4.75| ) co's(t), 
which is obtained by minimizing the solution (as a function of 7) of 

uj, = aAxt{l)^^s%\ . (5.36) 
V Xo / 

Taking the 6 = form for the minimum is at 7 = 1/2, and is given by 
which is identical to the true u)w>. 
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The second critical exponent ujcit) introduced in Sec. ^?7t corresponds to the rightmost 
singularity of the anomalous dimension. In the colhnear model the integrated gluon 
distribution can be defined as 

g^tj ■■= P dt exp l^^t + A, J'" as(r) drj 3^^{t) = ^ (5.38) 

because gio{t) can be shown to satisfy, in the collinear limit, the usual DGLAP equation 
with anomalous dimension [l/u + Ai)as- The singularity of the anomalous dimension (at 
energy scale k"^) 

7f(«s(t)):=^ln<7.(t) (5.39) 
is at the point u = ujdt) where guj{t) goes to zero, i.e., where A{t) goes to zero. 

It is easier to discuss ujc{t) in the cut-off case, in which as{t) = for t < i. In this 
case ( |5.38D allows us to write 

u-J{t) = A{t) + B{t) = as{t)g{t) + B{t) . (5.40) 

By substituting (|5.4CI|) in the first order differential equation ( ^.11|) it is not difficult to 
show that both g{t) and B{t) are continuous at t = t and exponentially behaved for t < t: 
g B exp (i(l + uj)t). It follows that 

For LiJ = ujp, 3^ujp{t) corresponds to the ground state of the system and have no zeroes: 



^ :r^,(t+) = — (^(t+) + Bit)) = — 



a;p(l + 



B{t) , (5.42) 



hence B{t) ^ and also g^jfit) 7^ 0. In order to have A{t~^) = one has to increase the 
depth of the potential well, i.e., to lower u < u^. This means that ujcii) < up. 

The relationships just found {u^""^'''^ = ujs(t), (t)) represent two extreme 

cases of the boundary condition dependence of up. If the coupling agit) > is positive 
but has intermediate size and shape for t < i, we expect in general that 

Uc{t) < tup < Us{t) . (5.43) 
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5.5 Perturbative regime: u;-expansion and WKB limit 



Approximate homogeneous solutions in the large-t region can be found, as in the full 
small-x equation of Chap. H, by the method of the 7-representation ([4.59|) and cj-expansion 
( [4.66| ). The regular solution is approximated, for ^ 1, by the expression: 



(7 



2^ buj 



with 



xr(7) 



where, for the coUinear model, 

Xo(7) 



1 



1 + 



+ 



1 



+ 



7 + t^ 
Ai-b 



(5.44) 

(5.45) 

(5.46a) 
(5.46b) 



(^+1^)2 ■ (l-^+i,;) 

In Chap. ^ this representation was extensively studied, and shown to be a solution of the 
problem up to next-to-leading order in uj and to all orders for the coUinear structure. But 
this was only guaranteed to work true for small values of u (whereas for the continuation 
with the DGLAP anomalous dimensions it is useful to be able to access high u as well). 
Also there was no way of determining the coefficient of any higher-order error introduced 
by the procedure. A comparison of this representation with the exact solution, as is 
possible in the coUinear model, is therefore important (cfr. Sec. |5.6.1j ). 

The expressions ( p.44| - 3^ ) can be further specialized in the large-t limit, where 
is dominated by a saddle point at 7 = 7a;(t) : but = ^(7, uj) and yields 



buo 



2vr|xL(7)l 



7TT Gxp / 7(r) dr . 



(5.47) 



We can obtain the large-t behaviour of the right-regular solution by means of the WKB ap- 
proximation for solving the differential equation ( ^.19| ). In fact, on the basis of Eqs. ( p.l9| ) 
and ( ^.28|) , one can prove the asymptotic expansion 



Mil 

2 Kit) 



exp 



nir) dr 



X 



(5.48) 



where K{t) is defined in terms of the effective potential (|5.19b| ) as 
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and is the zero of the WKB momentum: ^{1^) = 0. 

Eq. (|5.47| ) and ( ^.491) are asymptotic representations of the same function and are 
related by 

7 = ^1 + ^) + ■ (5.50) 

It is interesting to note that in the colhnear model, because of the differential equation 



( 5.19 ), the present method yields the irregular solution also, which is obtained by just 
changing the sign of the WKB momentum /t(t), i.e.. 



g^,(t) = -./-!^exp(+ f K{T)dT\ . (5.51) 



UJ \j 2K{t) 

This is useful for evaluating the Green's function according to Eq. (|5.34]) . 



5.6 Test of the u?-expansion 
5.6.1 Right-regular solution 

The possibility of evaluating (at least numerically) the GGF of the colhnear model to 
arbitrary precision, allows us to test the validity of the cu-expansion approximation method 
introduced in Sec. by comparing the exact right-regular solution of the differential 
equation ( ^.11| ), or equivalently (|5.19|) , with the 7-representation ( |5.44|) where ^^(7) is 



determined by the NLx truncation of the ti;-expansion ( ^.45] ). Fig. shows the exact 
and approximated function 3^R(t) for u = 0.15. One sees that for larger values of t, the 
a;-expansion is in good agreement with the exact solution, while for smaller t, where the 
solutions oscillate the results from the tu-expansion are slightly out of phase with the exact 
solution. In general we are interested in the behaviour to the right of the rightmost zero. 



5.6.2 Critical exponent 

One quantitative test of the w-expansion concerns the critical exponent Uc, i.e., the value of 
UJ at which the effective anomalous dimension of Eq. ( [5. 391 ) diverges. We recall that this is 



connected with the position of the rightmost zero of the regular solution: 3^(^^(t) = 0. The 
a;-expansion determination of the position of this zero, i.e., the zero of 3^R{t) determined 
by means of the 7-representation ( 5.44|) and the cj-expansion ( ^.45 ), involves a small 
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Figure 5.3: Right-regular solution from the explicit solution of the differential equation and 
from the cu-expansion in arbitrary normalization; shown for iv = 0.15. 

error which we call 6uJc- Fig. ^]4| shows 5uJc/uJc as a function of Uc-, and we see that the 
relative error on ojc goes roughly as cj^, or equivalently as . This corresponds to a NNLx 
correction and is beyond our level of approximation. We note also that even for relatively 
large values of ~ 0.3, the relative correction remains of the order of 5% which is quite 
acceptable. In other words the NNLx correction that arises is not accompanied by a large 
coefficient. 

5.6.3 Anomalous dimension 

Our second quantitative test of the cj-expansion concerns the anomalous dimension 7(^, as 
defined in Eq. ( p.39| ). The error in the cu-expansion result for the anomalous dimension, ^7, 
is plotted in Fig. |5.5| as a function of as for two values of uj. Let us first concentrate on the 
region for < 0.01. We see that the error is roughly independent of and proportional 
to al- In other words the difference between the exact result and the w-expansion is a term 
of 0{a1). We recall that the terms that we wish to include properly are the leading terms 
(as/cu)"', the NL terms as(«s/t^)" and the collinear terms {as/^^)oj^- A first correction of 
0{a1) is consistent with all these terms having been correctly included. 

For (Ss > 0.01 we see that the error in the anomalous dimension has a more complicated 
behaviour: it changes sign (the dip) and then diverges, at the same point as the divergence 
in the anomalous dimension itself (solid curve): this is just a consequence of the position 
of the divergence of the anomalous dimension from the cj-expansion being slightly different 
from the true position 
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Figure 5.4: The error, 6uJc in the determination of uJc within the uj-expansion. Shown for both 
the unintegrated and integrated solutions. 




Figure 5.5: Error in the effective anomalous dimension of Eq. ( }5.3£ } from the co-expansion 
plotted against as{t). 



127 



5.7 High energy behaviour 

It is widely believed that a two-scale process — described by a small-x equation of BFKL 
type — is perturbative for large enough t and to, while it becomes a strong coupling process 
if the energy is so large as to allow diffusion to small values of t ^ (fc^ ~ A^) 



In the collinear model the Green's function has the explicit expression (|5.34]), in which 



the strong-coupling information is clearly embodied in the reflection coefficient S{uj). 
Therefore it allows a direct study of the relative importance of the "perturbative" part 
® 3^/ and of the "strong-coupling" part S 3^r^ 3^ji, induced by diffusion through the 
boundary conditions at t = i. 

For large t and to; the factorized form ( |5.23|) of can be approximated by the WKB 



expression of Eq. (|5.48| , |5.5l|) , that we study in the special case t = Iq, so that the GGF 



in Y := In s/kk^ space is 



g(F;t,t)^ / — e^ 



(5.52) 



for large enough t. The (perturbative) first term is characterized by a saddle point at 
uj = Lj(t) > ujsit) > placed between the vertical asymptote of 1/Kuj(t) at ~ i-i-'s(^) 
on the left and the fast rising exponent e^^ on the right. Since around uJs we have 
fi;^(t) ~ const (t) [uj — ujs{t)], the saddle point condition reads 



^ (5.53) 



2(lj — u 



hence 



(^(F,t)^o;,(t) + ^. (5.54) 



The saddle point ( ^.54|) is self consistent provided the WKB approximation is valid, i.e. 



the variation of the WKB momentum K{t) is small on the typical wave-length of the 
solution: 

2 

— /— <1 ^ 2r< - uJs{t)-^ ^ asity^ 5.55a 

K I K \0 J 

and provided Y is not too small so that the saddle point is not too shallow: 

^ < Us{t) ~ a^{t) . (5.55b) 
The saddle point estimate is finally performed by means of the integral 
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rather than by naive quadratic fluctuations around x = 1/2 and yields 



g(P-*)(r;t,t) ^ ^=i=^e-W^ X [1 + OialY')] , x'L{t) ■= Min x(7, ^s(t)) • 

(5.57) 



The contribution of the (non perturbative) second term of Eq. ( |5.52| ) is evaluated with 
the residue method 



poles} 



(5.58) 



Retaining only the leading pole and observing that HujpiT) — (1 + '-^p)/2 + C(l/r), we 
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Figure 5.6: The contour of integration for the Green's function Saj(^; i, t), displaying the pole 
contribution of the non perturbative part and the background integral parallel to the imaginary 
axis relevant for the perturbative part of the GGF. 



end up with the result |^ 
9{Y;t,t)c. ^ 



-(l+ajp)i 



(5.59) 



showing that in the well-deflned intermediate regime ( |5.55| ) the exponent ujs{t) appears 
to be an observable quantity. Increasing Y the GGF enters a second, non perturbative 
regime where the second term of Eq. (|5.59| ) dominates on the flrst one, i.e.. 



(5.60) 
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Two remarks are in order. 
It is possible to calculate the higher order corrections to Eq. ( p. 57] )) as shown in Ref. 



[l + 0(a,V')] ^ EM^^ =exp r^-^-7'^ . (5.61) 

n=0 ■ \ / \ / 

It turns out that these corrections coincide with the "non-Regge correction" to the "Regge 
exponent" uJs{t)Y found by other authors 5£] in different but related context . We 
notice, however, that the true Regge contribution is the second term in Eq. ( |5.59D , which 
is of strong-coupling type, with a t-independent and eventually leading exponent up. 
The perturbative part, which dominates in the large-t limit, comes from the background 
integral and has no reason to be Regge behaved. 

• The quantity uJc{t) — the formal anomalous dimension singularity in the definition 
( |5.39|) — does not directly appear in the Y dependence, because the oscillating behaviour 
of the S^'s, which is relevant for increasing values of Y, is masked by the onset of up 
dominance for Y > Yt. 

We conclude that the two-scale Green's function shows a perturbative (non-Regge) 
regime where the exponent uJs{t)Y shows up with calculable corrections (Eq. ( [5.61| )), 
provided the parameter a^Y^ is small. Even before the latter gets large, at Y > Yt the 
Pomeron-dominated regime takes over, characterized by the regular solution, which is 
confined to the strong-coupling region of small fs. 



5.8 The tunneling mechanism 

The transition mechanism between the perturbative and non perturbative regimes de- 
scribed in the previous section reveals completely new features with respect to the well- 
known diffusion one [0], as we will see soon. We are going to present some numerical 
results, the analytical description being rather hard. 

Note first that several aspects of the two regimes are due to just running ctg effects. 
The latter enters at two levels: first, because of the "attraction" in the potential well 



of Fig. ^.2| , it modifies the traditional fixed-ag kind of diffusion, weighting it towards 
lower transverse scales. Secondly, because of the strong-coupling boundary condition, it 
also introduces a qualitatively new kind of diffusion, perhaps more properly referred to 
as 'tunneling': namely there is a certain Y = Yt, defined in Eq. ( ^.60|) , at which the 
non-perturbative pomeron suddenly takes over and beyond which the Green's function is 
dominated by the regular solution, and thus confined to small t-values. 

This 'tunneling' phenomenon is qualitatively different from diffusion, in so far as there 
is not a gradual decrease in the relevant scale for the evolution, until non-perturbative 
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scales become important, but rather there is a point beyond which low scales suddenly 
dominate, without intermediate scales ever having been relevant. 

To study these effects at a qualitative level it suffices to consider a very simplified 
version of the coUinear model: one which retains only the running of the coupling, but 



not the C(j-shifts of the 7 = 0, 1 poles, nor the Ai component of the NL corrections 
We then examine the solution to 



S{Y,t,to) = 6{t-to)+ [ dyK^Siy) 

Jo 



for this simplified kernel and we consider the effective exponent of the evolution 



(5.62) 



(5.63) 



as a function of Y, where S(Y,t,to) is simply S(y, t,to) with 6(t — to) subtracted, in 
analogy with Eq. ( ^.25| ). Fig. ^.7| illustrates the basic behaviour of uj^s for rather extreme 



^ 0.4 




Figure 5.7: The effective intercept ujes as a function ofY for t = 13 and two different values of 
the infra-red cutoff i. 

kinematics — not intended to be phenomeno logically relevant, but rather to show clearly 
the relevant features. Two values of the infra-red cutoff are considered. What is seen is 
that the exponent at first increases slowly and smoothly, and then at a certain threshold 
Y increases rapidly towards cup. The saturation of ujf occurs later in Y for increased i 
(i.e., decreased tup) as expected from ( |5.60|) . Traditional smooth diffusion into the infra- 
red would have led to the opposite behaviour, namely the higher i case (lower cjp) being 
saturated first. 
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In order to have a closer look to the underlying physical mechanism, we analyse the 
F-dependence of the so called Bartels' cigar consisting in the contour plot of the function 

S(Y-y;t,t'my;t',to) 

b{y,t,to) 

Once one has fixed Y, t and to, the y-section of Eq. ( p.64| ) as a function of t' describes the 
relative importance of the intermediate states of different transverse momentum k' = e*'/^ 
in the evolution of S from to Y . 

Let us start with a not too large rapidity interval Y and consider for simplicity to = 
t ^ 1. The first plot of Fig. ^.8| shows that the domain of the dominant paths in transverse 
momentum-rapidity space, represented by the light coloured region, is restricted around 
the value of t of the initial and final state, with a small spread. Increasing the value Y, 
the cigar gets thicker and thicker. In normal diffusion, the reach and the broadening of 
the cigar in the non perturbative region is gradual, according to its width ~ ■>/ as{t)Y. 

In the coUinear model instead, we observe (Fig. |5.8| ) a jump (tunneling) in the IR 
region of an increasing fraction of paths as Y gets larger, until all the relevant intermediate 
states belong to the strong coupling domain, thus providing the non perturbative regime 
mentioned in Sec. 

An impressing picture of the tunneling phenomenon can be obtained by plotting the 
distribution of transverse momenta in the middle of the evolution interval Y/2: Fig. ^ 
displays the contour plot of the function 

in t' and Y/2. In the left region (small evolution) the transverse momenta are spread 
around t in the perturbative region and the small-x effective growth exponent 

dln9(y;t,t) 

uJes{Y, t) := — (5.66) 

is given by the saddle point one Usit) plus the (y-derivative of the) corrections (|5.61| ). 
For Y > Yt tunneling is at work and the small-x exponent is given by Uf. 

To conclude this chapter, we stress that all the estimates based on the collinear model 
are by no means intended to be quantitative results. This model reveals itself useful in 
order to justify some of the assumptions adopted in the treatment of the complete small-x 
kernel such as the factorization properties and the cj-expansion, and to clarify the role of 
the critical exponents ujs(t), udt) and up. 

In particular, for two-scales processes, we have elucidated a completely new transition 
mechanism into the non perturbative domain which, for large enough energies Y >t{l + 
u!p)/{uf> — ujs(t)) and for sizeable values of — ujs{t) may be spectacular. 
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Figure 5.8: Domains of the dominant paths in transverse momentum t' and rapidity y space 
for the evolution of the GGF in Y between two hard and equal scales t = t^. In the first plots 
the relevant paths stays in the perturbative region, their spread increasing with Y. The last 
plots show the sharp transition of the paths in the non perturbative region corresponding to the 
pomeron dominated cross section. 
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Figure 5.9: Middle section of the Bartels' cigars of Fig. 5^ for continuous value ofY and the 
corresponding small-x effective growth exponent 



Chapter 6 



Conclusive remarks 



This last and short chapter summarizes the basic results of the work presented in this 
thesis and indicates the main open questions and directions for future studies as well. 

We have investigated the properties of QCD in the kinematic region of small-x with 
particular attention to deeply inelastic scattering of electrons on hadrons and in particular 
to the physics at HERA. The most precise experimental results in the small-x domain 
concerns the measurements of the F2 structure function which clearly exhibits a marked 
rise (Fig. towards low values of a; in a wide range of Q^. Even though the standard 
DGLAP evolution equations are able to describe the whole set of data, nevertheless they 



have to be supplied with an initial condition (see Eq. (|3.8| )) for the partonic distribution 
functions which cannot be explained within the DGLAP approach itself. Furthermore, 
there is a variety of two-scale processes, e.g. 7*7* scattering and forward-jets, whose 
x-behaviour is determined independently of Q^-evolution which is suppressed. 

The natural framework where to study the small-x behaviour of those inclusive pro- 
cesses is the BFKL one, where the large logarithmic coefficients of the perturbative series 
have been resummed at leading order and, recently, also at next-to-leading order. Even 
though the BFKL equation provides the right power-like shape for the small-x growth 
of structure functions, nevertheless the quantitative estimate of the growth exponent is 
absolutely not compatible with its measurements (Fig. This is true both at leading 
order — which overestimates the small- a; rise — and in the NLx approximation, where the 
corrections are so large and negative that, for the accessible values of Q^, the x-behaviour 
changes its slope (Fig. |3.13| ). 

The analysis performed in Sees. and have shown that the responsibles of such 
large corrections and of the ensuing instability of the BFKL hierarchy are the big collinear 
contributions to the coefficient kernels in the BFKL expansion. More precisely, the co- 
efficient kernels K^'^\k, k') of Eq. ( p. 11 ) are affected by single and double logarithms of 
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k"^ /k'"^ in the collinear limits k"^ ^ k'"^ and k"^ -C k'"^. The single logarithms are essentially 
those predicted by the renormalization group equation which are responsibles for the log- 
arithmic scaling violations in for the partonic densities. The double logarithms, on the 
other hand, are due to the mismatch between the high energy factorization scale Sq = kk^ 
— which is supposed to be the correct scale which allows the factorization formula ( p.9|) — 
and the Bjorken scale sb = Max(/c^, /cq) providing the relevant scaling variable xb — sb/s 
in the collinear limits k^ ^ k^ and k^ k^. 

Those large logarithms are unavoidable in the ctg-expansion of the BFKL kernel ( p.ll| ) 
and the corresponding series is oscillatory and unstable for not very small values of Og- In 
addition, the total cross section stemming from the naive NLx truncation of that series 
may turn out to be negative in the very large s-limit! A resummation of the large collinear 
contributions mentioned before is therefore mandatory. 

By a RG analysis of the collinear singularities of the GGF (|4.20| ) and hence of the kernel 
( [4.22| , ^^2^ it is possible to determine the coefficients of single and double logarithms of the 
transverse momenta k, feg to all orders in In 1/x. The resummation of those logarithms 
can be accomplished by a suitable shift of the 7-singularities of the eigenvalue functions 
xi"^(7) of the improved coefficient kernels of Eq. ( [4. 191 ) as indicated in Eqs. ( |4.26| ). By then 
incorporating the constraints of the known Lx and NLx BFKL eigenvalue functions, we 
are able to construct — up to NNLx corrections — the resummed Lx and NLx eigenvalues, 
given in Eqs. ( ^.361 ) and ( [4.4(j| ). 

Having determined the coUinear-resummed Lx and NLx kernels, it remains the ques- 
tion of determining the improved GGF, i.e., the solution of Eq. ( ^4.49|) . The technical 
difficulty in solving such an equation is mainly concerned with the running of as(^^) 
and in particular with the singularity of its perturbative expression at the Landau point 
k"^ = A^. In practice, a suitable regularization of as(fc^) has to be adopted. 

At this point we proceed in two directions: the first one consists in considering the large 
t = Infc^/A^ regime (t ^ to) where perturbative (regularization independent) features and 
strong coupling (regularization dependent) ones are argued to factorize and therefore can 
be considered independently. The second one is to define a simplified small-x model which 
can be exactly solved (at least numerically) where the qualitative features of the real case 
can be studied as well as the issue of factorization. 

The complete small-x equation has been studied in Chap. ^ where, in the regime 
k^ ^ kQ, the GGF is argued to assume the factorized form ([4.57|) in terms of a perturbative 
factor F^{k) which obeys the homogeneous equation (|4.58|) . The solution of the latter 
has been given in terms of the 7-representation ( [4. 591 ) and of the tu-expansion ( [4.661 ). 

The merit of the cu-expansion is to provide an effective eigenvalue function Xuj{'y) 
whose shape is stable even for sizeable values of the expansion parameter uj (Fig. 
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and whose error stemming from neglecting NNLx contributions is much smaller than the 
corresponding NNLx truncation of the BFKL kernel in the ttg-expansion. As far as the 
stability of the shape is concerned, this means that a minimum in the function 7 t— > Xuj{'y) 
is always present, even for large values of u. Because of that, a stable saddle point exists, 
preventing us from negative cross sections and providing a small-x growth exponent uJs{t) 
— defined in Eq. ( [4 .751 ) and plotted in Fig. — which should be observable in some 
intermediate-x regime (see Sec. p.7| ). In the range — 4 40 GeV^, i.e., — 0.2 0.3, 
the critical exponent ujs{t) ^ 0.27 ^ 0.32 is not far from the corresponding values of the 
HERA data. 

It is to be noted that, while uJs{t) decreases for large t, the trend of the experimental 
growth exponent is to increas with Q^. This is due to the scaling violations, which the 
DGLAP approach takes into account by resumming the In contributions which are not 
directly incorporated in uJs{t). A more direct observation of the critical exponent ujs{t) is 
expected from two-scale processes where the Q^-evolution is less important (cfr. beginning 
of Chap. I). 

The error coming from the NL truncation in the cj-expansion is uniformly (i.e., inde- 
pendently of 7) (9(ct;^), the neglected coefficients having no 7 = nor 7 = 1 singularities 
at all. This error is therefore of the same size as the ambiguity in the definition of Xf. 
The corresponding error in the saddle point condition ( [4.62| ) is a roughly 7-independent 
change of scale A{bt) = 0{uj), affecting the anomalous dimension with a relative uncer- 
tainty Aos/as ~ o^s^ while for the critical exponent Aujs{t) ~ uJs{t)At. The important 
thing to note is that all these errors are NNLx and that their coefficients do not show a 
strong t-dependence. 

In connection with the ambiguities of the finite order truncation, one can consider 
also the renormalization scale and scheme dependence as well as the resummation scheme 
dependence. In Sec. |4.10| it is shown that such dependences are very small and do not 
affect considerably the critical tu-exponents Us and Uc- 

The collinear model presented in Chap. ^ allows us to investigate in more detail the 
physical consequences of the coUinear-resummed kernel. First of all, the model offers a 
method of analysis of non perturbative effects on the observable quantities, in particular 



on the very small- a; growth exponent (Sec. |5.4| ). The main results of that analysis is 
that, in two-scale processes, for large enough values of the external scales fc^ and fcg and 
for not too large values of the energy s, there exists a regime governed by perturbative 
physics where the effective small-x growth exponent is given by ujs{t). A preliminary 
phenomenological analysis seems to indicate that the x-growth in 7*7* and forward-jet 
processes are similar, and consistent with the value quoted before. For increasing values 
of the energy, there is a transition point beyond which most of the exchanged gluons lie in 
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the strong coupling region and the growth exponent is given by the "pomeron" which is 
of non perturbative nature. This transition mechanism does not take place continuously 
hke "diffusion" . For this reason, we refer to it as "tunneling" . 

The study of the RG improved small-x equation has just begun and leaves many open 
questions yet, which hopefully will find an answer after further investigation. Just to 
mention some of them, there is to clarify the role of the critical exponent C0c{t) — the true 
singularity of the resummed anomalous dimension. A second point concerns the evaluation 
of higher twist contributions stemming from poles in the effective eigenvalue functions at 
values of 7 < —uj/2 and j > 1 + uj/2. A preliminary study of such corrections can be 
done in the coUinear model by suitably modifying the pole-structure of the eigenvalue. 
A major task concerns the implementation of unitarity which however should involve the 
study of strong coupling effects. 

To conclude, small-x physics emerges from our work as a rather unique framework 
where the perturbative and strong coupling features play a complementary role. I think 
we have shown that by combining Inl/x and InQ^ evolutions we can progress a great 
deal towards understanding the perturbative aspects of the problem. Further progress 
will require a deeper understanding of non perturbative physics. 



Appendix A 



The Structure Functions in DIS 



According to the definition of tlie kinematics variables of DIS in Sec. [1.4| (see Fig. |l.l|) ,the 
transition matrix element is given, in the lowest order of the EM interaction, by 

{I' X\T\IH) =uAv'i)'^e^''u„{Vi)^{X\{-e)J,mV2A) , (A.l) 

where a', a and A are the spin components of the scattered electron, initial electron and 
target proton respectively, e is the electronic EM charge and J^{x) is the quark part of 
the EM current. 

The inclusive unpolarized differential cross section (with respect to the outgoing elec- 
tron) is obtained by multiplying ( [A.l|) by the appropriate phase space factor and reads 

where we have defined the leptonic tensor 

cr.cr' 

=2 (Kp- + KK^ ) + ? (A.3b) 

and the hadronic tensor 

W,AP2, q) 5^(2vr)^5^(p2 + q-Vx)\ $^(^2, A| J^(0)|X) (X| 7,(0)^2, A) (A.4a) 

X A 

=i- j d'ze^'^^\Y.^p,A\[J,{z).UmP2A) ■ (A.4b) 
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Because of current conservation at the hadronic vertex, (9^J^(z)|X) = and hence 
Q^W^iv = = q^W^y. The most general form of the hadronic tensor fulfilhng the require- 
ments of Lorentz covariance, space-inversion and time-reversal conservation and current 
conservation is (p29 = P2 ' q) 



W^"" = Fi 



P2q V r 



P2-^<l'') , (A.5) 



where the Lorentz invariant coefficient Fi{x^Q'^ 



: I 



1,2 are called structure functions 



of the proton. It is worthwile to stress that the structure functions depends only on the 
Lorentz invariants of the lower vertex (the blob) in Fig. |1 . 1| , i.e., on x and Q^. 



Keeping p2 and q fixed, let's introduce the longitudinal polarization unit vector 



Q_ 

P2q 



{P2qy 



1 + 



P2 - 



q - €1 = , ei- eL 



(A.6) 



and a couple of transverse unit vectors ei and €2 such that 

q-ei = , eL-ei = , ■ = -1 , ei ■ €2 = , (« = 1, 2) . 

The sets of vectors {£1,62} and {£^,£1,62} form pseudo-orthonormal basis for the sub- 
spaces < P2, q >"*" and < q >-*- respectively. Defining the projectors 



n 



(A.7) 
(A.8) 

(A.9) 



in the longitudinal, transverse and orthogonal (to q^) space respectively, we can decompose 
the hadronic tensor as 



-F,n- + F2||(i + ^)nr 



1 

2x 



,/2 



1 + 1 F2 n^" 



(A.lOa) 
(A.lOb) 
(A.lOc) 



where the longitudinal structure function is given by 



FLix,Q''):= [l + %]F2~2xF, 



(A.ll) 
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If we consider virtual photon absorption by a protonQ, the ensuing cross section involves 
just the hadronic tensor ( [A .41) and the former assumes the form 

a{^*H)^>e''W,,, (A.12) 

being the polarization vector of the virtual photon. It is easy then to extract any of 
the structure functions simply by choosing the appropriate vector e^. If we parametrize 
the generic polarization vector = Ae^, + tet where is a unitary transverse vector, 
then *e^^e^ W^y gives 

• Fi for |A| = and |r| = 1; 

• Fl for |A| = 1 and |r| = 0; 

• F2 for |A| = |r| = 2a;/(l + g). 



A.l M,m ^ limit 

For high values of the energy a/s and of the momentum transfer Q with respect to the 
typical hadronic masses, the latter can safely be neglected, together with the electronic 
mass which is even smaller. In this case z/~^ = and for the relevant kinematical variables 
holds 

= xys . (A.13) 

In terms of these variables and of the azimuthal angle (f) in the CM frame, the phase space 
element in Eq. (|A.2| ) becomes 

,Q d 2sx^ d 

so that we can write the differential cross section with respect to x, and at fixed s 
as 

= ^,{xy^F, + {l-y)F^] (A.14a) 

' {[l + {l-yf]F^-y^FL] (A.14b) 



a 



xQ^ 

where a: = e'^/An is the EM fine structure constant. Note that in this massless case 
Fl = F2-2xFi. 



^Dealing with an off-shell incoming particle, a particular flux convention has to be adopted. 



Appendix B 



Dimensional regularization 



The dimensional regularization technique consists in performing the loop and phase space 
(divergent) integrals of Feynman diagrams and squared matrix elements in a modified 
space-time of generic dimension D ^ 4, where such integrals exist. The ensuing expres- 
sions are then analytically continued for complex values of D. 

For our purposes, the main modifications with respect to the D = 4 case concern the 
phase space measure d^*^"^ and the expression of the coupling Og. 

The generally adopted convention for momentum integration is to replace the volume 
element as 



(27r)4 {2tt)D 

Introducing the (half) extra-dimension parameter e : D =: 4+2e, the phase space measure 

n+1 , , ^ , n 1 n+1 



2+n) 



]_r[^r[ , ^ 1 TT TT (^9) 



(2.)2 (47r)- l\ f \ (27r)2 (2.)^ (47r)" l\ fj^ (27r)^ 

i.e., only the (27r)^ denominator under dkj is modified. 

In addition, in going from Lx to NLx approximation, also the longitudinal measure is 
modified. In particular the two-body phase space reads 

and the three-body one is 

d0(^) = — ^^^^^^ . (B.4) 

^ 167rs2zi(l-2i) (27r)4+4= ^ ' 
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In 4 + 2£ dimensions the QCD coupling constant g acquire a mass-dimension [g] — —e, 
therefore a massive parameter /i has to be introduced, so to consider g/i^ as the effective 
dimensionless couphng. The convention generally adopted in the papers reporting the 
calculations of the NLx BFKL kernel is to use the MS renormalization scheme and to 
define the dimensionless colour strength as 



which reduces to the customary expression in the e — ^ limit. Accordingly, we have 
adopted the following definition for the Born impact factor: 




(B.5) 



1 




As transverse measure we keep the natural one 



dfc:=d2+^^fc 



5{k): 




Appendix C 



Integral representations for functions 
in transverse momentum space 

C.l The leading BFKL eigenvalue function 

Given a set Q, of complex functions defined on a vectorial space V, an operator /C acting 
on Q is said to be rotationally invariant if it commutes with all the rotation operators 
on V. If /C(fc, k') is the kernel of the hnear integral operator /C, then JC is rotationally 
invariant if and only if 

)C{Rk,Rk') =IC{k,k') , k,k' eV (C.l) 

where R is the generic rotation operator in V. 

Let Sx denote the scaling operator on fl defined by 

[Sxfm:=f{Xk) , fen. (C.2) 

An operator K, is said to be scale invariant if it commutes with Sx for all A > 0. In a 
d-dimensional vector space V, the integral operator K, is scale invariant if and only if 

IC{Xk, Xk') = A-'^/C(fe, k') . (C.3) 

Under rather general conditions, if the operator /C is rotationally invariant, then it 
admits a complete set of eigenfunction {fj : j & J} of the form 

fj{k) = ^,{k')YM (C.4) 
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where (j) denotes the set of angular variables and Yj is a spherical harmonic in d dimensions. 
Furthermore, if /C is also scale-invariant, the radial eigenf unctions \l/ can only be powers 
of the variable . 

The BFKL operator /C acts on functions defined in the d = 2 dimensional transverse 
space V\ the only angular variable is the azimuthal angle 0. The action of /C on a generic 
function f eVL can be explicitly written as 



(C.5) 



the first and second term corresponding to the real and virtual part of the kernel re- 
spectively. It is easy to check that, for fixed coupling ots-, the BFKL kernel defines a 
rotationally and scale invariant operator. 

Since in two dimensions spherical harmonics corresponds to complex exponentials 



the eigenfunctions of /C are of the form 

f'y,m{k) = k 



7-1' 



and a complete set is obtained by allowing m G Z and 7 = 1/2 + iz/ : u G 
eigenfunctions ( p.7|) are normalized according to 



5\k -ko) = J2 ^ f,Ak)fU{ko) 



(C.6) 

(C.7) 
The 

(C.8) 



meZ 2 

and satisfy the orthonormality conditions 

j dfc f*^rn{k)f^',m'{k) = 271 8 mm' " 



-i=\ + iy , -i' = \ + w' . (C.9) 



The IR finiteness of the kernel can be checked directly on the eigenfunctions. For 
simplicity we restrict ourselves to the m = case. By using the dimensional regularization 
m d = 2 + 2e and the formula 

d^-"^^fc 1 _ T{l-a-P- e)T{a + e)T{P + e) . 

T{a + [5 + 2e)T{l-a)T{l- p) ' 



we get: 

d^+^'k' 



(C.IO) 



}C^''\k,k'){k'y-^ 



- + V(l)-V'(7)-^(l-7) + (^(^) 



(C.ll) 
(C.12) 
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where the real and virtual BFKL kernel are defined in Eqs. ( |2.54| ) and ip : = F'/F is the 
logarithmic derivative of the gamma function. It is evident the singular behaviour of both 
the real and virtual terms m. d = 2 dimension corresponding to the poles in e. However, 
when adding the two contributions, the poles cancel and in the e — limit we recover 
the m = eigenvalue function xil) '■= x(7)0) = 2^/'(l) — ipi^j) — ip{l — 7). When acting 
on the general eigenfunction the kernel yields 

J dk }C{k, k') f^^Uk') = «sX(7, m) (C.13) 

771 771 

X(7, m) = 2ij{l) - V.(7 + y ) - ^(1 - 7 + y ) ■ (C.14) 



C.2 Mellin representation in transverse momentum 
space 



The BFKL kernel can be expressed in terms of its eigenfunctions ( C7 ) and its eigenvalues 
( |C.13|) by means of the spectral representation 

/-i+ioo 1 

/c(fc,fc') = E / ^ A,4fc)«sx(7,m)/;^(fe) (C.15) 

In particular, the azimuthal average of the kernel, defined by 

— /c(fc,fc') , e = kk' (C.16) 

2n 

has (A;^)"^"^ as eigenfunctions and xil) '■ = x(7)0) as eigenvalue function. Thanks to 
Eq. (|C.15|) , it can be given the integral representation 

;C(fc^fc-) = 4, / ^ «sX(7), (C.17) 

nk^ yi_ioo 27ri V 

a.x(7) = Jdk (^^^ \{k\k'') , (C.18) 

showing that x{l) is nothing but the Mellin transform of the azimuthally averaged kernel 
/C. This is due to the fact that the radial eigenfunctions of K, are powers of k^. 

The Mellin transformation allows to diagonalize convolutions of function in both lon- 
gitudinal and transverse space, such as the high energy factorization formula (|2.20|) . By 
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defining 



F2,M--= I ^ f'-x-F,{x,Q^) (C.19a) 



g2 \A2j X 



(ruj{V-= ~ P ^(^'Z') ' ^ = T ' P = — (C.19b) 



-7 ^1 



1 dk (C.19c) 



the mentioned factorization formula can be written in tlie simple algebraic form ( p.22 ) or 
as a fc-integral if the Mellin transform is carried out only in the longitudinal variables x 
and z. 

The Mellin transforms ( |C.19| ) can be inverted by using 

F2{x,Q')= r^^P: (^Y I Px-F,M (C.20a) 



27ri V A V J 27ri 

where the integration in the complex tu-plane is parallel to the imaginary axis and to the 
right of all the singularities of the integrating function. 
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